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I.  PROGRAM  SLMMARY  1969-73 


A.  Most  Significant  Accomplishments 

1.  The  establishment  of  a  totally  interdisciplinary  effort  in  which  the 
areas  of  physiology,  psychology,  and  engineering  interface  to  provide 
an  optimal  research  capability.  This  interdisciplinary  effort  has  also 
provided  a  major  impact  on  the  research  and  teaching  capabilities  of 
this  University's  Biomedical  Engineering  Program. 

2.  Development  of  a  center  of  excellence  as  a  resource  for  support  of 
the  DOD  mission. 

3.  Establislunent  of  effective  communication  links  and  collaborative 
efforts  among  other  Air  Force  Contractors  and  Air  Force  Research 
Laboratories  to  provide  continuing  two  way  flow  of  information. 

4.  Organization  and  presentation  of  the  first  national  symposium  on 
Chronically  Implanted  Cardiovascular  Instrumentation,  and  publi¬ 
cation  of  the  results  by  Academic  Press. 

5.  Development  of  a  Standard  Experimental  Preparation  (SEP)  for  eval¬ 
uation  of  the  Physiological  and  Biomechanical  responses  to  environ¬ 
mental  stresses.  The  SEP  includes  awake  and  anesthetized  animals 
with  acutely  and  chronically  placed  instrumentation  providing 

a  broad  range  of  variables  appropriate  for  analysis  of  the  effects 
of  single  and  multiple  stressful  stimuli. 

6.  Implementation  of  the  SEP  for  study  of  awake,  chronically  instru¬ 
mented  animals  exposed  to  vibration  stress,  the  first  series  of  in¬ 
vestigations  of  this  type  to  be  reported. 

7.  Combining  the  capabilities  supported  by  this  contract  with  additional 
support  from  the  Vibration  Branch,  Wr ight-Patterson  AFB,  to  measure 
multiple  physiological  variables  in  a  monkey  performing  a  tracking 


_> 
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task.  Preliminary  results  were  described  in  a  Technical  Report 
published  by  the  Vibrations  Branch  and  a  second  Technical  Report 
is  in  press. 

8.  Determination  of  the  facts  that  (1)  a  major  factor  underlying  the 
cardiovascular  effects  of  vibration  in  doga  is  peak  net  transmitted 
force;  (2)  maximum  force  transmission  in  dogs  and  monkeys  is  in  the 
3  to  6  Hz  range  for  acceleration  levels  of  1G  and  above;  (3)  use  of 
frequency,  displacement,  or  acceleration  alone  does  not  give  the 
consistent  response  pattern  observed  with  peak  net  force  as  the  para¬ 
meter;  and  (4)  with  repeated  exposure,  force  transmission  and  heart 
rate  extremes  decrease,  just  as  in  human  subjects. 

9.  Development  of  a  closed-loop,  distributed  parameter  computer  model 
which  accounts  for  some  of  the  wave  characteristics  of  transient 
flow  and  for  the  motion  induced  in  the  circulatory  system  by  vibra¬ 
tion. 

10.  Development  of  a  novel  approach  for  quantifying  the  coefficients  of 

a  two-mass,  single-degree-of-f reedom, bio-vibrational  model  and  demon¬ 
stration  of  its  validity  for  predicting  the  mechanical  impedance  of 
the  sitting  primate  during  sinuaiodal  vibration. 

11.  Demonstration  in  the  intact  left  ventricle  that  the  elastic  modulus 
(E)  is  a  linear  function  of  the  mean  tangential  stress  (o)  through¬ 
out  iaovoiumetric.  systole,  and  the  slope  (K,  modulus  of  stiffness) 
of  the  E  vs.  o  curve  is  18.8,  with  one  standard  deviation  of  0.9. 

This  is  the  first  known  evaluation  of  the  relationship  between  the 
elastic  modulus  and  the  mean  tangential  stress  of  canine  cardiac 
muscle. 

12.  Demonstration  of  the  relationship  between  a  decrease  in  the  amplitude 
of  the  somatosensory  cortical  evoked  potential  and  the  development 

of  irreversible  hemorrhagic  shock. 
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C.  Summary  of  Experimental  Program 

Of  the  various  stresses  imposed  upon  subjects  performing  critical  tasks, 
time  dependent  accelerations  (vibration)  can  produce  some  of  the  greatest  dis¬ 
comfort  and  decrements  in  performance.  The  effects  of  this  stress  on  the  car¬ 
diovascular  system  are  being  investigated  at  this  laboratory.  Instrumented 
animals  are  used  to  evaluate  the  cardiovascular  changes  produced  by  whole  body, 
sirusoidal  vibration.  Experiments  are  designed  to  evaLuate  the  role  of  tbe  five 
major  mechanisms  which  current  evidence  indicates  are  Che  main  factors  produc¬ 
ing  these  changes.  The  mechanisms  important  for  vibrational  stre  s  of  a  given 
frequency,  displacement,  acceleration,  transmitted  force,  axis,  and  duration  are: 

1.  Reaction  of  the  fluid  and  vessel  system  as  determined  by  its  mass, 
elastic,  and  damping  characteristics. 

2.  Reaction  of  large  body  organ  systems  and  the  musculoskeletal  sys¬ 
tem  as  determined  by  their  mass,  elastic  and  damping  characteristics. 

3.  Reaction  of  the  mechanoreceptors  of  the  cardiovascular  regulatory 
system  as  determined  by  the  relationships  to  the  fluid  and  vessel 
system,  large  body  organ  systems,  and  the  musculoskeletal  Bystem 
and  its  mechanoreceptors. 

4.  Reaction  of  the  endocrine,  metabolic  -  hematological  system  as  deter¬ 
mined  by  the  relationships  to  the  musculoskeletal  and  large  body 
organ  systems. 

5.  Response  mediation  of  the  central  nervous  system  and  the  psychophvsio- 
loglcal  system  of  conscious  subjects  acting  through  the  above  pathways 
to  produce  further  alterations  in  metabolic  and  cardiovascular  function. 

-Cardiovascular,  hematological,  and  mechanical  variables  are  used  to  evalu¬ 
ate  the  overall  stress  level  and  the  role  of  the  above  mechanisms.  Anesthe¬ 
tized  and  av.ke  animals  are  used  for  both  short  and  long  term  (30  sec.  to  6  hr.) 
vibration  stress  exposures. 
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General  Experimental  Protocol:  Responses  of  acutely  and  chronically  instru¬ 


mented  dogs,  primates  and  pig9  have  been  studied.  Special  emphasis  has  been 
placed  on  the  chronically  implanted  preparation  which  permits  repeated  awake 
and  anesthetized  tests  on  the  same  animal.  Additional  advantages  include 
the  closed  chest  protocol  and  more  stable  transducer  characteristics,  espe¬ 
cially  important  for  the  vibration  environment.  Animals  are  restrained  on 
the  hydraulic  table  with  the  spine  vertical,  and  vibrated  along  that  axis. 

Short  term  tests  consist  of  sequential  exposures  to  sinusoidal  it  nation  of 
2-12  Hz  and  acceleration  amplitudes  ranging  from  0.5  G  to  3.0  G  for  10  sec., 

2  min.,  and  5  min.  intervals.  Long  term  exposures  include  durations  vary¬ 
ing  from  30  min.  to  6  hrs.  at  a  selected  vibration  frequency  and  constant 
acceleration  amplitude.  The  highest  g-level  for  the  awake  animal  is  1.5  g, 
and  for  the  anesthetized  animal  is  3.0  g.  The  cardiovascular  variables  in¬ 
clude  heart  rate,  flow  velocity  from  the  aorta  and  peripheral  arteries,  ar¬ 
terial  pressures  from  acceleration- insensitive  intraaortic  and  intraventri¬ 
cular  transducers,  0^  consumption  and  body  temp,  and  blood  chemistries  such 
as  blood  gases,  cortisol,  free  fatty  acids,  and  glucose,  measured  on  samples 
withdrawn  from  indwelling  atrial  and  arterial  catheterB.  The  design  of  the 
vibration  exciters  permits  recording  of  force  transmitted  between  subject  and 
table,  and  also  the  vibration  exciter  velocity.  From  these  variables  and  their 
phase  relationship,  complex  mechanical  velocity  impedance  can  be  determined. 

In  addition,  accelerometers  implanted  on  large  body  organs  provide  Information 
on  the  displacement  of  these  organs  re'ative  to  the  animal  torso  and  the  input 
vibration.  A  Raytheon  704  data  acquisition  system  is  being  programed  for  off¬ 
line  and  on-line  analysis  of  the  massive  amount  of  cardiovascular  and  mechanical 
data.  Additional  aspects  of  the  research  program  include  analog  and  digital 
computer  simulation  of  the  effects  of  sinusoidal  vibration  on  the  hydraulic 
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and  mechanical  aspects  of  the  cardiovascular  system.  The  results  from  the  com¬ 
puter  study  guide  the  planning  of  animal  experiments  and  aid  in  the  analysis  of 
experimental  data. 

Results 

A.  Short  term  exposure:  In  an  effort  to  identify  the  parameters  of  vibration 
which  may  serve  best  as  an  index  of  stress  to  the  cardiovascular  system,  special 
emphasis  has  been  given  to  the  role  of  the  peak  net  force  (PNF)  transmitted  be¬ 
tween  the  animal  and  vibration  exciter.  A  whole  body  resonance  at  3-4  Hz  was 
found  for  seven  awake  dogs  tested  at  1G  acceleration  amplitude.  At  this  frequency 
(f ) ,  peak  net  force  was  about  1.7  times  the  body  weight.  For  frequencies  up  to 
and  including  6  Hz,  the  animal  still  transmitted  a  p^.ak  force  greater  than  body 
weight  tBW) .  Beyond  6  Hz  this  force  continually  decreased  until  at  30  Hz  a 
peak  force  of  only  about  15%  of  body  weight  resulted.  The  effects  of  anesthesia 
as  reflected  in  the  force  transmission  data  show,  in  general,  that  the  anes¬ 
thetized  animal  transmits  approximately  30%  more  force  at  resonance  and  has  a 
resonant  frequency  1  to  2  Hz  lower  than  the  awake  animal.  At  higher  frequencies 
there  have  been  no  consistent  force  transmission  differences  between  awake  and 
anesthetized  dogs.  When  the  percent  change  in  mean  aortic  flow  (MAF)  for  the  awake 
dogs  was  plotted  versus  PNF/BW,  a  linear  relationship  wa6  found.  Changes  in  MAF 
ranged  from  -30  to  +150%.  At  this  time,  it  appears  PNF/BW  describes  best  the 
changes  in  MAF,  with  only  secondary  correlations  existing  for  frequency,  acceleration 
and  displacement.  The  mechanisms  responsible  for  these  changes  in  dogs  (7  dogs, 

10  experiments)  can  be  provisionally  catalogued  as  follows: 

1.  Changes  in  mean  heart  rate  (MHR)  and  relatively  small  changes  in 

stroke  volume:  Data  from  three  tests  showed  a  linear  relationship 
between  %  change  in  MHR  and  peak  net  force,  with  stroke  volume 
relatively  unchanged. 
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Changes  in  stroke  volume  (SV)  and  relatively  small  changes  in  MHR: 

Data  from  four  tests  showed  a  linear  relationship  between  7.  change 
in  stroke  volume  and  peak  net  force,  with  X  change  in  MHR  relatively 
unchanged. 

3.  Equal  changes  in  MHR  and  stroke  volume:  Data  from  three  tests  showed 
linear  relationship  between  both  %  change  in  MHR,  and  X  change  in 
stroke  volume  and  peak  net  force. 

Furthermore,  the  type  of  changes  in  MAF  could  be  correlated  with  pre-vibra¬ 
tion  values  of  MHR  and  SV.  For  example,  if  the  previbration  MHR  was  low,  in¬ 
creases  in  MAF  occurred  principally  by  an  increase  in  HR.  (Category  1).  If, 
however,  previbration  MHR  was  high,  increases  in  MAF  resulted  from  an  increase 
in  SV.  (Category  2) . 

Similar  observations  were  made  for  the  data  of  the  anesthetized  dogs  except 
that  Category  3  was  not  as  evident  and  that  the  previbration  anesthesia  profile 
also  influenced  '■he  dura.  The  slope  of  the  line  of  MAF  vs  PNF/BW  (best  fit)  for 
the  data  from  10  tests  on  7  awake  dogs  and  7  anesthetized  dogs  was  analyzed.  For 
those  experiments  in  which  increases  in  MAF  were  produced  principally  by  increases 
in  SV,  no  significant  differences  were  found  between  awake  and  anesthetized  dogs  run 
on  the  same  vibration  t3sts.  When  an  increase  in  MHR  was  used  to  increase  MAF, 
the  slope  of  the  line  ot  MAF  vs  (PNF/BW)  for  the  awake  animals  was  67%  higher  than 
that  of  the  anesthetized  animals  data.  Changes  in  MAF  were  found  also  to  vary 
linearly  with  the  log  (MHR/vlb .  freq.). 

In  addition,  transmissibility  of  the  root  of  the  aorta  was  found  to  be  lin¬ 
early  related  to  PNF/BW,  extremely  dependent  on  vibration  frequency,  and  to  have 
a  resonance  at  approximately  A  Hz,  at  which  frequency  the  applied  vibration  ampli¬ 
tude  (acceleration)  was  amplified  approximately  2.4  times. 
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B.  Long  term  exposure:  For  90  m'n.  exposures  at  a  constant  frequency  and  G 
level,  no  consistent  changes  in  MAF  were  observed.  The  fluctuations  which 
occurred  over  time  were  more  numerous  at  4  Hz  (Q.75G;  than  at  6  or  10  Hz.  Of 
the  limited  tests  performed,  one  out  of  5  awake  dogs  changed  MAF  through  signif¬ 
icant  changes  in  MHR  and  lesser  changes  in  SV.  Two  out  of  five  dogs  changed 
MAF  through  significant  changes  ir.  SV  and  lesser  changes  in  MHR.  No  consistant 
changes  in  PNF/BW  vs  time  were  observed.  Similarly,  no  definite  trends  were 
seen  in  tracking  monkeys  repeatedly  exposed  to  multiple  frequencies  up  to  3 
hours  duration. 

C.  Computer  Model:  A  lumped  parameter  model  was  used  to  aid  in  determining 

the  passive  mechanical  and  hydraulic  contribution  to  the  total  physiological  re¬ 
sponse  of  the  awake  dogs  tj  vibration.  The  effects  of  vibration  are  modeled  by 
pressure  source  terms  generated  from  values  of  the  vibration  g-level  or  result¬ 
ing  transmitted  force.  These  pressure  disturbances  propagate  and  reflect,  adding 
and  subtracting  with  the  pressure  generated  by  the  action  of  the  heart,  and  as 

a  result  produce  changes  in  flow  as  well.  The  value  of  PNF/BW  vs  f  from  the 
animal  experiments  was  used  to  establish  the  values  of  the  source  terms  for  the 
different  simulated  vibration  frequencies.  Heart  rate  was  varied  in  the  model 
as  observed  in  the  animal  experiment.  The  computer  response  of  MAF  vs  PNF/BW 
was  similar  to  those  of  the  animal  experiments  with  the  actual  %  change  values 
close  to  those  found  fer  the  anesthetized  dogs,  in  which  SV  remained  relatively 
constant  and  HR  changed,  e.g.  maximum  %  change  in  MAF  was  approximately  25%. 
Discussion: 

From  a  preliminary  analysis  of  the  data,  it  appears  that  of  the  changes  in 
MAF  (via  changes  in  SV)  of  awake  and  anesthetized  dogs  subjected  to  vibration, 
approximately  60  to  70%  is  due  to  one  or  a  combination  of  the  mechanisms  3,  4 
and  5  (i.e.  reaction  of  mechanoreceptors  of  the  cardiovascular  regulatory  sys- 
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remaining  30  to  40%  due  to  mechanisms  i  and  2  (i.e.  reaction  of  the  fluid  vessel- 
system,  and  the  large  body  organ  and  musculoskeletal  system) .  A  comparison  of 
the  computer  model  studies  with  the  anesthetized  animal  data  in  which  changes  in 
MAF  were  produced  principally  by  changes  in  MHR,  indicates  that  this  anesthetized 
preparation,  to  the  first  approximation,  responds  as  a  passive  "mechanical"  sys¬ 
tem.  Furthermore,  the  correlation  of  MAF  with  the  ratio  of  heart  frequency/table 
frequency  suggests  the  influence  of  the  phase  relationship  between  the  heart  cycle 
and  the  vibration  cycle  (through  the  pathways  of  mechanism  1)  in  producing  changes 
in  MAF. 

Analysis  of  these  data  indicates  that  in  addition  to  the  passive  mechanical 
aspects,  substantial  changes  are  produced  by  other  regulatory  feedback  mechanisms 
sensitive  to  vibrations.  These  vibration  receptors  may  lead  to  sympathetic  stimu¬ 
lation  and/or  possible  heterometric  regulation  in  the  heart.  The  relative  role 
of  these  various  active  mechanisms  in  producing  changes  in  MAF  is  certainly  not 
clear  at  this  time.  However,  similar  interactions  between  HR  and  SV  in  generat¬ 
ing  MAF  as  observed  for  these  vibration  studies  have  also  been  described  for  ex¬ 
ercise.  It  will  be  of  value  to  investigate  this  point  further.  Nevertheless, 
the  use  of  awake  animals  for  assessing  the  effects  of  vibrational  stress  on  the 
cardiovascular  system  has  been  found  most  valuable.  In  addition,  PNF/BW  is  an 
especially  important  parameter  in  assessing  changes  in  MAF  during  vibrational 
stress.  It  is  recommended  that  future  studies  should  investigate  the  neural  and 
mechanical  mediation  involved  in  the  response  of  the  cardiovascular  system  to  vi¬ 
bration. 
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II.  PROGRESS  REPORT  1972-73 


A.  Introduction 

Of  the  various  stresses  imposed  upon  subjects  performing  critical  tasks, 
time  dependent  accelerations  (vibration)  can  produce  some  of  the  greatest  dis¬ 
comfort  and  decrements  in  performance.  The  effects  of  this  stress  on  the 
cardiovascui.tr  system  are  being  Investigated  and  are  the  subject  of  this  re¬ 
port.  Instrumented  animals  are  used  to  evaluate  the  cardiovascular  changes  pro¬ 
duced  by  whole  body,  sinusoidal  vibration.  Experiments  are  designed  to  eval¬ 
uate  the  role  of  the  five  major  mechanisms  which  current  evidence  indicates  are 
the  main  factors  producing  these  changes.  The  mechanisms  important  for  vibra¬ 
tional  stress  of  a  given  frequency,  displacement,  acceleration,  transmitted 
force,  axis,  and  duration  are: 

1.  Reaction  of  the  fluid  and  vessel  system  as  determined  by  its  mass, 
elastic,  ar.d  damping  characteristics. 

2.  Reaction  of  large  body  organ  systems  and  the  musculoskeletal  sys¬ 
tem  as  determined  by  their  mass,  elastic  and  damping  characteristics. 

3.  Reaction  of  the  mechanoreceptors  of  the  cardiovascular  regulatory 
system  as  determined  by  the  relationships  to  the  fluid  and  vessel 
system  large  body  organ  systems,  and  the  musculoskeletal  system 
and  its  mechanoreceptors. 

4.  Reaction  of  the  hormonal,  metabolic  -  hematological  systems  as  deter¬ 
mined  by  the  relationships  to  the  musculoskeletal  and  large  body  organ 
systems . 

5.  Response  mediation  of  the  central  nervous  system  and  the  psychophysio- 
logical  system  of  conscious  subjects  acting  through  the  above  pathways 
to  produce  further  alterations  in  metabolic  and  cardiovascular  function. 

A  proposed  model  of  the  interaction  of  these  five  mechanisms  is  illustrated 
in  the  block  diagram  of  Figure  1.  Cardiovascular,  metabolic,  hematological,  and 
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mechanical  variables  are  used  to  evaluate  the  overall  stress  level  and  the 


role  of  the  above  mechanisms.  The  following  is  a  review  of  the  available 
information  on  the  relationship  between  the  mechanical  and  physiological 
responses  in  animals  and  man  exposed  to  vibrational  stress. 

Vibration  Parameters:  Any  given  vibration  is  characterized  by  its  frequency, 
displacement,  velocity,  acceleration  and  waveform.  While  these  parameters  may 
define  the  vibration  applied  to  a  subject,  they  are  not  necessarily  transmitted 
equally  to  all  parts  of  the  subject.  For  example  the  head  might  well  vibrate 
with  a  displacement,  velocity,  acceleration,  and  waveform  vastly  different  from 
that  applied  at  the  subject-vibration  exciter  interface.  Thus,  while  certain 
parameters  define  the  applied  vibration,  other  parameters,  dependent  upon  the 
properties  of  the  vibrated  subject,  are  necessary  to  characterize  the  response 
of  the  subject  to  the  applied  vibration.  Such  parameters  include  the  trans¬ 
mitted  force  at  the  subject-vibration  exciter  interface,  resonant  frequencies 
(those  frequencies  corresponding  to  the  maxima  of  total,  integrated  whole-bocy 
displacement  with  respect  to  the  applied  displacement) ,  and  the  transmissibllity 
of  a  particular  body  segment  (transmissibllity  is  defined  simply  as  the  ratio 
of  the  vibration  amplitude  of  any  given  segment,  e.g.  the  head,  heart,  shoulder, 
etc.,  to  the  applied  vibration  amplitude,  measured  at  vibration  exciter). 

When  attempting  to  determine  the  vibration  response  of  any  system  the  ini¬ 
tial  question  arises  as  to  what  scheme  of  applied  vibration  parameters  should 
be  used?  Consider  the  case  of  sinusoidal  vibration  where  the  displacement  (x) 

is  x  ■  x  sin  wt,  the  velocity  (v)  is  v  ■  x  w  cos  wt  and  the  acceleration  (a) 
o  o 

2 

is  a  *  -x  w  sin  wt.  The  peak  acceleration  amplitude  relative  to  the  gravita- 

0  2 

X  w 

tional  constant  g  is  often  times  written  as  a  ■  o  and  referred  to  as  "G 

~g~ 

level" . 

Some  investigators  favor  applying  various  vibration  frequencies  each  hav¬ 
ing  the  same  displacement  amplitude,  x  ;  however,  for  this  approach  as  frequency 
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increases  the  peak,  acceleration,  a,  increases  with  frequency  squared.  There¬ 
fore  going  from  1  to  10  Hz  requires  a  peak  acceleration  100  times  greater  at 
10  Hz  than  at  1  Hz.  Thus,  this  approach  necessitates  a  very  low  acceleration 
amplitude  at  the  lowest  frequency  in  order  to  avoid  lethal  accelerations  at 
the  higher  frequencies.  A  similar,  but  less  extreme,  condition  results  by 
applying  different  frequencies  each  with  the  same  peak  velocity,  xqw.  For 
this  case  as  ^he  frequency  is  increased  the  displacement  amplitude  must  be 
reduced  inversely  with  the  frequency  and  the  acceleration  amplitude  will  in¬ 
crease  linearly  with  frequency.  The  most  connr'n  method  employed  by  biovibra- 

tory  investigators  has  been  to  apply  various  frequencies  each  having  the  same 

2 

acceleration  amplitude,  xow  .  For  increasing  frequencies,  this  approach  re- 

g 

quires  decreasing  the  displacement  amplitude  by  the  inverse  frequency  squared. 
For  any  given  peak  acceleration  the  largest  displacement  amplitude  is  therefore 
required  at  the  lowest  frequency,  and  decreases  rapidly  with  increasing  fre¬ 
quency.  From  a  mechanical  viewpoint  this  constant  peak  acceleration  approach 
has  two  important  advantages;  namely  (1)  mechanical  measurements  can  be  com¬ 
pared  directly  to  the  response  of  an  "inert"  mass  of  weight  equal  to  that  of 
the  subject,  since  in  this  case  the  peak  force  transmitted  between  the  vibra¬ 
tion  exciter  and  the  imaginary  inert  mass  is  independent  of  frequency,  and 
(2)  from  a  practical  view  point,  the  displacement,  velocity,  and  acceleration 
limits  of  the  vibration  exciter  are  not  as  easily  exceeded  using  this  approach 
as  compared  to  applying  various  frequencies  at  either  a  constant  displacement 
or  velocity  amplitude.  Further  information  pertinent  to  this  discussion  is 
given  by  Lange  and  Edwards  (1970). 

A  unique  feature  incorporated  with  the  Wenner-Gren  Research  Laboratory's 
vibration  exciters  is  a  force  measurement  unit  (force  cell)  which,  during  vi¬ 
bration,  provides  an  analog  output  voltage  proportional  to  the  net  force  trans- 
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quires  decreasing  the  displacement  amplitude  by  the  inverse  frequency  squared. 
For  any  given  peak  acceleration  the  largest  displacement  amplitude  is  therefore 
required  at  the  lowest  frequency,  and  decreases  r:?idly  with  increasing  fre¬ 
quency.  From  a  mechanical  viewpoint  this  constant  peak  acceleration  approach 
has  two  important  advantages;  namely  (1)  mechanical  measurements  can  be  com¬ 
pared  directly  to  the  response  of  an  "inert"  mass  of  weight  equal  to  that  of 
the  subject,  since  in  this  case  the  peak  force  transmitted  between  the  vibra¬ 
tion  exciter  and  the  imaginary  inert  mass  is  independent  of  frequency,  and 
(2)  from  a  practical  view  point ,  the  displacement,  velocity,  and  acceleration 
limits  of  the  vibration  exciter  are  not  as  easily  exceeded  using  this  approach 
as  compared  to  applying  various  frequencies  at  either  a  constant  displacement 
or  velocity  amplitude.  Further  information  pertinent  to  th^s  discussion  is 
given  by  Lange  and  Edwards  (1970) . 

A  unique  feature  incorporated  with  the  Wenner-Gren  Research  Laboratory's 
vibration  exciters  is  a  force  measurement  unit  (force  cell)  which,  during  vi¬ 
bration,  provides  an  analog  output  voltage  proportional  to  the  net  force  trans¬ 
mission  between  subject  and  restraint,  i.c.  the  "live-load"  force  transmission 
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occurring  at  the  animal-seat  interface,  Sharp  (1963).  The  force  transmitted 
between  animal  and  support  (restraint)  during  vibration  reflects  the  relative 
movement  of  internal  body  organ  systems  and  parts.  Correlations  of  large 
transmitted  force  values  to  corresponding  changes  in  circulatory  parameters 
can  indicate  that  movement  of  internal  systems  caused  the  observed  circula¬ 
tory  changes  via  vessel  occlusion,  mechanical  stimulation  of  the  heart  and/ 
or  vessels,  reciprocating  internal  pressures,  mechanical  stimulation  of  regu¬ 
latory  receptors,  etc. 

Literature  Review:  Numerous  invectig;  lions  have  been  concerned  with  establish¬ 
ing  the  mechanical  response  of  man  and  animals  to  whole-body  mechanical  vibra¬ 
tion.  Coermann  (1962)  applied  the  theory  of  mechanical  impedance  of  systems 
with  one  or  more  degrees  of  freedom  to  sitting  and  standing  man  vibrated  at 
frequencies  from  1  to  20  Hz.  he  was  one  of  the  first  investigators  to  demon¬ 
strate  that  the  force  transmission  response  of  man  to  whole-body  vibration 
could  be  approximated  by  that  of  a  simple  mechanical  system  consisting  of  inter¬ 
connected  linear  inertial,  elastic,  and  damping  components.  Experimental  data 
confirmed  results  predicted  by  his  model,  i.e.,  that  man  had  a  primary  whole- 
body  resonance  in  the  4  to  8  Hz  frequency  range.  Thus,  it  was  quantitatively 
established  that  vibration  frequencies  between  4  and  8  Hz  induced  the  greatest 
internal  body  organ  movement  relative  to  that  input  to  the  body  by  the  vibra¬ 
tion  exciter.  White,  et  al.  (1962)  reported  the  effects  of  mechanically  induced 
vibration  upon  the  human  abdomen  by  measuring  the  pressure  at  the  rectal  end 
of  the  colon  sigmoideum.  For  the  semi -supine  posture  a  mean  "peak-to-peak" 
vibration  induced  fluctuation  of  57  mmHg  at  4  and  4.5  Hz  was  recorded.  The 
pri.  saure  resonances  were  correlated  to  resonances  of  body  organs  in  the  thoraco¬ 
abdominal  region.  Clark  et  al.  (1963)  determined  the  effects  of  forced  vibra¬ 
tion  by  measuring  from  semi-supine  man  the  circumferential  body  strain  of 
chest,  abdomen,  pelvis,  and  thigh.  All  subjects  showed  maximum  body  strain 
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near  7  Hz.  Zechman  et  al.  (1965)  reported  maxima  of  seated  man's  transpulmonary 
pressure  during  5  Hz  vibration  for  0.5  G  acceleration  amplitude.  Lange  and 
Edwards  (1970)  vibrated  human  subjects  in  the  supine  attitude  from  2  to  14  H2 
and  0.2  to  0.5  G.  A  principal  body  resonance  was  reported  between  5  and  7.5 
Hz,  where  also  the  maximum  body  strain  was  recorded.  Using  skeletal  muscle 
tonus  as  a  parameter  they  found  that  muscle  tightening  increased  the  whole-body 
transmitted  force  by  as  much  as  50%  during  vibration  at  the  resonant  frequency. 
Edwards  and  Knapp  (3972)  recorded  the  force  transmitted  during  repeated  vibra¬ 
tion  of  awake  dogs.  Results  indicated  reductions  as  great  as  36%  in  the  trans¬ 
mitted  force  recorded  from  the  initial  to  the  7th  exposure  of  one  animal  during 
4  Hz  vibration  at  0.7g.  Evces  and  McElhaney  (1971)  studied  the  effects  of 
drugs  on  the  mechanical  response  of  dogs  to  vibration  and  reported  that  drugs 
produced  measurable  changes. 

The  nonlinear  behavior  of  biomechanical  systems  has  been  the  subject  of 
several  investigations.  Krause  and  Lange  (1963)  described  such  behavior  during 
vibration  at  various  acceleration  amplitudes  over  the  2  to  15  Hz  frequency  range. 
Wittman  and  Phillips  (1969)  also  presented  data  that  indicated  the  recognition 
of  nonlinearities  as  seen  in  the  measured  impedance  response  of  human  subjects 
seated  erect  in  several  impact  and  vibrational  environments.  Vykukal  (1968) 
combined  a  sustained  acceleration  (via  a  centrifuge)  with  vibration  to  study 
man's  nonlinear  mechanical  response,  and  Vogt  et  al.  (1973)  used  a  similar  ex¬ 
perimental  technique  to  conclude  that  the  human  body  stiffens  with  increasing 
sustained  accelerations,  and  that  the  resonant  frequency  increases  from  6  Hz 
under  normal  gravity  to  3,  11,  13,  and  15  Hz  with  2  Gx,  3  Gx,  4  Gx,  and  5  Gx, 
respectively. 

Efforts  have  also  been  directed  toward  establishing  tolerance  criteria  for 
man  exposed  to  vibration.  Magid  and  Coermann  (1960)  vibrated  seated  human,  sub¬ 
jects  at  various  frequencies,  each  with  increasing  amplitude  at  0.75  mm  per  sec¬ 
ond  until  the  subject  reported  that  he  felt  actual  bodily  harm  would  result  (this 


17 


limit  being  beyond  discomfort) .  Results  indicated  minimal  short  time  tolerance 
between  4  and  8  He  at  accelerations  between  1.5  and  2.0  G.  Ziegenruecker  and 
Magid  (1959)  in  a  similar  investigation  presented  a  "short  time  tolerance"  curve 
to  sinusoidal  vibration  by  gathering  subjective  data  from  man  under  vibration  at 
various  frequencies  and  amplitudes,  i.e. ,  at  various  accelerations. 

Von  Gierke  (1971)  chaired  a  symposium  on  "Biodynamic  Models  and  their 
Application"  for  establishment  of  environmental  exposure  limits,  for  intrepre- 
tation  of  animal,  dummy,  and  operational  experiments,  mechanical  characteri¬ 
sation  of  living  tissue  and  isolated  organs,  models  to  describe  man's  response 
to  impact,  blast,  and  acoustic  energy,  and  performance  in  biodynamic  environ¬ 
ments.  Because  of  man's  obvious  limitations  as  a  subject  in  biodynamic  experi¬ 
ments,  animal  surrogates  have  been  used  in  many  studies.  Bantle  (1971)  used 
three  vibration  conditions  to  determine  which  was  the  most  damaging  to  4  1/2 
and  7  day  mouse  embryos.  It  was  found  that  a  20  Hz  (visceral  resonant  frequency) 
vibration  was  the  most  destructive,  and  the  41/4  day  embryo  was  the  most  labile. 
He  reported  that  vibration  typically  inhibited  embryonic  growth.  Ashe  et  al. 
(1961)  recorded  body  weight  and  camperature  of  rats  subjected  to  whole  body 
vibration  and  showed  that  observed  and  measured  differences  to  changes  in  fre¬ 
quency  and  amplitude  indicated  that  both  factors  played  a  role,  but  calculations 
indicated  that  the  common  denominator  was  not  simply  acceleration. 

In  animals  subjected  to  very  high  levels  of  acceleration  (0.6  to  50  g)  and 
various  frequencies,  tissue  damage  and  survival  correlated  with  mechanical  re- 
ponse  (Pape,  et  al. ,  1963;  Boorstin,  et  al.,  1966).  Clearly,  there  are  signifi¬ 
cant  body  resonances,  particularly  in  the  range  below  12  Hz.  Organ  system  re¬ 
sonances  account,  in  part,  for  the  subjective  discomforts  experienced  by  humans 
(Pape,  et  al.,  1963).  The  type  of  symptoms  occurring  in  human  subjects  appear 
related  to  both  displacement  and  acceleration  (Linder,  1962).  The  symptoms 
appeared  somewhat  dependent  upon  the  acceleration  level  reached  in  the  studies 
of  Temple,  et  al.  (1964)  also. 
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Several  studies  have  dealt  with  analytical  models  which  help  to  explain  the 
mechanical  response  of  subjects  to  vibration.  Broderson  and  Von  Gierke  (1971)  in¬ 
vestigated  the  biomechanical  parameters  of  the  sitting  rhesus  monkey  and  their  tem¬ 
poral  changes  during  sinusoidal  vibration  at  frequencies  from  6  to  30  Hz.  Results 
showed  impedance  magnitude  and  phase  angle  to  decrease  with  time.  The  authors  also 
derived  a  simple  model  for  the  impedance  response  of  the  sitting  rhesus  monkey. 
Payne  and  Band  (1971)  developed  a  "four-degree-of-f reedom"  lumped  parameter  model 
for  seated  man  to  study  the  problem  of  his  response  to  vertical  accelerations. 

Broderson  (1972)  conducted  a  study  investigating  the  biothermal  response  of 
rhesus  monkeys  to  mechanical  vibration.  He  concluded  that  core  temperature  in- 
cr'j.',»‘d  particularly  at  resonance  and  varied  with  frequency.  Prolonged  l.Og  vibra- 
••  30  Hz  were  shown  to  increase  core  temperature  by  about  1/2  C,  and  the 
■AL-y  mechanism  was  speculated  to  be  frictional  dissipation.  Lafferty  et  al. 
(1973)  and  Edwards  and  Lafferty  (1973)  also  conducted  studies  concerned  with  the 
response  of  the  rhesus  monkey  to  mechanical  vibration.  The  mechanical  impedance 
response,  its  variation  with  time  during  vibration  at  a  given  frequency  and  ac  el- 
eration  amplitude,  and  a  mechanical  model  to  predict  the  impedance  response  at 
different  frequencies  and  accelerations  up  to  1.0  G  are  described  by  these  investi¬ 
gators.  The  effect  of  vibration  on  a  performance  task,  and  related  changes  in 
physiological  parameters  were  also  discussed. 

The  effects  of  vibration  on  performance  has  been  the  subject  of  many  investi¬ 
gations.  Lange  and  Coermann  (1963)  measured  the  visual  acuity  of  human  subjects 
during  whole-body  vibration  at  various  low  frequencies.  They  reported  that  maxi¬ 
mum  decrements  occurred  at  those  frequencies  where  resonance  of  whole  body  and  or¬ 
gan  complexes  had  previously  been  determined,  i.e.,  4-8  Hz.  O'Braint  and  Ohl- 
baum  (1970)  presented  a  method  to  describe  the  visual  efficiency  of  vibrated  man 
in  terms  of  a  Performance  Index. 

Studies  of  human  tracking  performance  showed  inconsistent  effects  when 
different  frequencies  and  accelerations  were  the  parameters  (Frazer,  et  al., 
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1961).  In  seated  subjects  performing  a  light-pattern  psychomotot  task  using 
a  control  stick,  Fraser,  et  al.  (1961)  found  displacement  to  be  the  predo¬ 
minant  factor  in  determining  performance  decrement.  An  increase  in  either 
displacement  or  frequency  was  significantly  related  to  performance  decrement, 
but  a  combination  of  any  frequency  with  the  lowest  displacement  did  not  pro¬ 
duce  a  performance  decrement.  Increase  in  displacement  at  any  frequency  did 
produce  a  decrement  which  could  be  intensified  by  change:,  in  frequency.  The 
frequencies  used  were  2,  4,  7,  and  12  Hz  at  displacements  o'  0.1'>5,  0  25,  0.375, 
and  0.5  in. double  amplitudes.  Maximum  peak  acceleration  was  3.67  g.  Regression 
analysis  showed  mean  error  was  most  closely  correlated  to  the  product  of  the 
displacement  and  the  square  root  of  the  frequency.  Interestingly,  the  frictional 
resistance  component  of  surface  impedance  is  also  proportional  to  the  square 
root  of  frequency.  A  function  of  both  displacement  and  impedance  was  suggested 
as  the  most  satisfactory  parameter  (Fraser,  et  al.,  1961). 

In  related  studies,  Catterson,  et  al.  (1962)  extended  the  analysis  of  the 
role  of  displacement  and  frequency  in  human  tracking  performance.  Exposure 
duration  23  min  at  6,  8,  11,  and  15  Hz  and  displacements  of  0.13  and  0.26  in. 
double  amplitude.  A  decrease  in  error  occurred  at  all  frequencies  for  0.13  in. 
displacement,  and  an  increase  in  error  occurred  at  all  frequencies  for  0.26  in. 
displacement.  As  pointed  out  by  Catterson,  et  al.,  acceleration  is  proportional 
to  displacement  times  frequency  squared,  and  functions  determined  by  acceleration 
must  reflect  changes  in  frequency  to  a  far  greater  extent  than  amplitude. 

Grether  et  al.  (1972)  extended  an  earlier  study  in  which  they  reported  that 
a  combination  of  heat,  nois®,  and  vibration  stress  had  no  greater,  and  for  some 
measures  slightly  less,  effect  on  physiological  and  performance  functions  than 
did  the  same  levels  of  heat  or  vibration  alone.  The  study  searched  for  possible 
explanatory  mechanisms  for  the  apparently  antagonistic  stress  interaction.  It 
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was  postulated  that  the  heat  stress  night  in  some  way  have  modified  the  bio- 
dynamic  response  of  seated  man  so  as  to  reduce  the  vibration  transmitted  to 
the  subject's  body.  However,  transmissibility  data  for  the  shoulder  shewed 
no  evidence  of  any  such  effect.  It  was  postulated  that  increased  motivation 
and  application  of  effort  during  the  combined  stress  exposures  might  be  the 
responsible  mechanism.  Sommer  and  Harris  (1972)  exposed  10  human  subjects 
to  vibration  (5  Hz,  0.25  g)  and  noise  (110  dB)  while  asking  them  to  perform 
a  mental  arithemetric  task.  Results  suggested  that  the  phase  of  the  cir¬ 
cadian  cycle  could  be  a  variable  to  be  considered  in  studies  on  the  effects 
of  stress  on  human  performance. 

The  effects  of  vibration  on  behavior  has  also  been  investigated.  Wike 
and  Wike  (1972)  reported  results  from  seven  experiments  on  low-frequency, 
whole-body  vibration  of  rats  as  related  to  their  escape  conditioning.  They 
concluded  chat  vibration  of  sufficient  amplitude  is  aversive  for  rats  and  that 
its  brief  termination  is  reinforcing.  A  direct  relationship  was  evident  be¬ 
tween  amplitude  and  the  number  u£  escape  responses  when  frequency  was  held 
constant  and  amplitude  was  varied. 

The  area  receiving  the  most  attention  is  the  effect  of  vibration  on  the 
cardiovascular  system.  Roberts  and  Dines  (Oct.,  1966)  administered  propra¬ 
nolol  and  atropine  to  anesthetized  vibrated  dogs  and  measured  heart  rate, 
cardiac  output,  total  peripheral  resistance,  and  (d’.Vdt)^^  le£t  ventri¬ 
cular  pressure.  They  concluded  that  tachycardia  during  vibration  resulted 
from  a  decrease  in  vagal  efferent  activity  to  the  heart.  In  vibration  tests 
■using  human  subjects  Roberts  et  al .  (1969)  concluded:  "(1)  clinically  use¬ 
ful  ECG's  can  be  obtained  from  human  subjects  while  they  are  bein^  vibrated, 
even  at  relatively  severe  intensities,  if  our  recommendations  to  accomplish 
this  are  followed,  (2)  short  time  vibration  in  the  practical  range  uBed  were 
only  moderately  stressing  to  healthy  subjects,  (3)  vibration  can  elicit  pre- 
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mature  extrasystoles  in  man,  (4)  changes  in  configuration  of  ECG's  observed 
during  vibration  were  transient,  terminating  with  termination  of  vibration, 
and  (5)  patterns  of  pulse  rate  changes  during  vibration  are  variable."  Ed¬ 
wards  (1970)  and  Edwards,  McCutcheon  and  Knapp  (1972)  performed  tests  in  which 
dogs  were  vibrated  at  1-3G  from  2  to  12  Hz,  during  which  chronically  implanted 
transducers  enabled  recording  cardiovascular  data,  and  showed  the  greatest  changes 
occurred  in  the  3  to  9  Hz  frequency  range  for  any  given  acceleration  amplitude. 

3G  vibration  to  4  Hz  produced  (1)  a  maximum  aortic  peak  flow  rate  of  more  than 
twice  the  maximum  recorded  from  the  resting  animal,  and  (2)  a  minimum  aortic 
peak  flow  rate  of  only  10%  of  th>  minimum  observed  in  the  resting  animal.  The 
results  also  indicated  a  correlation  between  relative  motion  of  internal  masses 
and  large  changes  in  the  circulatory  system  during  the  vibration  exposure. 

Hooks,  Nerem,  and  Benson  (1969)  performed  a  theoretical  analysis  of  pul¬ 
satile  flow  in  a  rigid  tube  with  longitudinal  vibration  by  Introducing  as  a 
boundary  condition  a  sinusoidal  wall  velocity.  Using  a  sinusoidal  pressure 
gradient  of  constant  amplitude  they  predicted  flow  rate  increases  of  400%  or 
reductions  to  zero  relative  to  the  flow  rate  which  occurred  without  vibration. 
Catnil,  Knapp,  and  Collins  (1971)  used  the  analog  computer  to  model  the  response 
of  the  human  cardiovascular  system  to  1-10  Hz  sinusoidal  vibration.  Their 
study  showed  the  most  pronounced  effects  at  3-4  Hz  for  all  acceleration  ampli¬ 
tudes  investigated.  For  an  acceleration  amplitude  of  3G  they  reported  the 
maximum  difference  in  peak  aortic  pressure  during  one  heart  cycle  tripled  and 
the  maximum  peek  aortic  flow  doubled,  relative  to  the  control  values. 

Arntzenius,  Koopa,  and  Hugenholtz  (1969),  Blok  et  al .  (1969),  Verdouw, 
Arntzenius,  and  Noordergraaf  (1969),  and  Jackson  (1971)  reported  a  technique 
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employing  Body  Acceleration  Synchronous  with  the  Heartbeat  (BASH)  in  which 
a  particular  type  of  whole-body  vibration  applied  to  both  animals  and  man 
was  reported  to  enhance  cardiovascular  performance.  Ihe  use  of  vibration 
as  a  noninvasive  cardiac  assist  device  was  thus  reported. 

Any  effect  primarily  dependent  on  the  degree  of  mechanical  stimulation 
that  increases  with  increasing  amplitude  of  movement  should  increase 
with  decreasing  frequency  if  peak  acceleration  ia  held  constant.  If  ampli¬ 
tude  is  held  constant,  the  effect  will  increase  with  increasing  peak  accele¬ 
ration  and  consequently  increase  with  increasing  frequency.  This  was  found 
to  be  true  for  the  responses  of  hyperventilation  (Dixon,  et  al.,  1961)  and 
peripheral  vasodilation  (Liedtke  and  Schmid,  1969).  Gaeuman,  et  al .  (1962) 
cited  the  previous  reports  of  Fraser,  et  al.  (1961)  and  Catterson,  et  al. 

(1962)  to  substantiate  their  experimental  design  for  the  measurement  of  oxy¬ 
gen  consumption.  The  subjects  were  exposed  to  a  constant  peak-to-peak  dis¬ 
placement  of  0.264  in. double  amplitude,  and  frequencies  of  2,  6,  8,  11,  and 
15  Hz.  Since  displacement  amplitude  was  constant,  acceleration  should  increase 
with  frequency.  The  values  given  for  peak  acceleration  at  each  frequency 
were  0.05,  0.45,  0.82,  1.55,  and  2.38  g  respectively.  A  clear  frequency- 
displacement  dependence  of  oxygen  consumption  was  demonstrated,  suggesting 
constant  displacement  was  a  useful  parameter.  Hoover  and  Ashe  (1962)  used 
displacements  of  0.124  and  0.250  in. double  amplitude  at  frequencies  from  2 
to  15  Hz.  The  peak  acceleration  range  was  0.03  to  2.88  g.  The  magnitude  of 
the  resulting  hyperventilation  showed  a  mixed  frequency-displacement  depen¬ 
dence.  Respiratory  rate  changes  were  predominant  at  a  displacement  of  0.124  in., 
while  tidal  volume  changes  were  predominant  at  a  displacement  of  0.250  in.  They 
concluded  the  respiratory  response  tended  to  be  related  primarily  to  table 


displacement . 


An  initial  tachycardia,  hypotension,  increased  cardiac  output  and  de¬ 
creased  total  peripheral  resistance  was  the  characteristic  change  in  the 
anesthetized  dogs  studied  by  Dines,  et  al.  (1965).  The  parameters  were 
frequencies  of  4,  7,  and  11  Hz  at  a  displacement  of  0.5  in.  double  ampli¬ 
tude.  They  averaged  the  responses  at  all  three  frequencies,  and  did  state 
the  blood  pressure  dip  was  not  evident  at  11  Hz. 

In  contrast,  a  large  body  of  data  supports  the  primacy  of  a  frequency- 
acceleration  dependence  of  the  above  responses  (Guignard,  I960;  Goldman  &  von 
Gierke,  I960;  Zechman,  et  al.,  1965;  Cavagna,  1970;  Edwards,  et  al.,  1972). 

Acceleration  dependence  of  cardiovascular  changes  has  also  been  demon¬ 
strated  in  animals  (Hood  and  Higgins,  1965;  Liedtke  and  Schmid,  1969).  More 
marked  physiological  effects  were  found  at  1.2  g  and  8  and  10  Hz  than  at 
0.6  g  and  frequencies  on  either  side  of  that  range  (Hood,  et  al.,  1966). 
Clark,  et  al.  (1967),  in  studies  of  anesthetized  dogs,  found  significantly 
greater  change  in  multiple  systemic  cardiovascular  variables  at  2  than  at 
1  g  for  frequencies  of  10,  12,  and  14  Hz.  (Apparently  2  g  runs  were  not 
made  at  6  and  8  Hz). 

The  changes  in  oxygen  consumption  observed  by  Duffner,  et  al.  (1962) 
were  maximal  at  the  lovest  frequencies  and  greater  at  0.35  than  at  0.15  g 
acceleration . 

Another  type  of  physiological  interaction  of  frequency  and  displace¬ 
ment  is  illustrated  by  the  machanoreceptor  responses.  Although  the  exact 
route  by  which  muscle  stimulation  is  connected  to  caruiovascular  and  res¬ 
piratory  responses  is  unclear,  increases  in  systemic  ventilation  and  blood 
pressure  are  produced  by  muscle  stimulation  such  as  rhythmic  manual  squeez¬ 
ing  at  2  to  3  Hz,  and  the  response  i  <;  abolished  by  cutting  the  nerve  (Guy¬ 
ton,  et  al.,  1972,  and  McCloskey,  er  nl.,  1972).  Of  course  the  carotid 
sinu9  also  has  distinctive  frequency-amplitude  characteristics,  but  its 


participation  in  the  responses  is  likely  to  be  minimal  since  its  sensitivity 
is  probably  reduced  (Bristow,  et  al.,  1971). 

Acceleration  was  a  much  more  important  factor  than  duration  in  cats  ex¬ 
posed  to  severe  vibration  (10  to  15  g;  Pape,  et  al.,  1963).  It  would  thus 
appear  that  dogmatic  statements  about  the  interaction  of  displacement,  fre¬ 
quency,  and  acceleration  are  inappropriate,  and  each  must  be  considered  in 
evaluating  the  effects  of  vibration  stimulus. 

Previous  studies  of  the  integrated  response  of  the  cardiovascular  sys¬ 
tem  to  vibration  have  also  produced  variable  and  often  conflicting  results. 
Components  of  cardiovascular  function  including  heart  rate,  blood  pressure, 
and  cardiac  output  have  been  reported  to  decrease,  increase :  o’"  remain  un¬ 
changed  du;ing  vibration.  Hood  and  Higgins  (1965)  studied  anesthetized, 
supine  dogs  restrained  In  a  form-fitting  metal  frame  and  vibrated  in  the  x- 
axis  at  selected  frequency  and  acceleration  levels.  Cardiac  output  increased 
as  a  result  of  increased  heart  rate  at  the  higher  acceleration  levels;  i.e. 
stroke  volume  remained  constant.  At  6  Hz,  mean  arterial  pressure  increased 
slightly;  at  10  Hz,  it  decreased.  Total  peripheral  resistance  was  not  signif¬ 
icantly  altered  at  6  Hz  but  increased  at  10  Hz.  Overall,  they  were  unable 
to  identify  any  main  effects  of  frequency  or  duration.  Main  effects  of  accele¬ 
ration  were  present  for  most  of  the  variables.  There  were  interactions  be¬ 
tween  acceleration  and  frequency  for  heart  rate  and  mean  arterial  pressure. 

In  other  studies,  tachycardia,  decreased  arterial  pressure,  increased 
cardiac  output,  and  decreased  total  peripheral  resistance  were  indicated 
(Clark,  et  al.,  1967;  Dines,  et  al. ,  1965;  Hoover,  et  al.,  1961).  Clark, 
et  al.,  1967  found  mean  arterial  blood  pressure  decreased  initially  and  then 
returned  to  control  levels;  stroke  index  remained  constant.  The  effects 
were  greater  at  higher  G  levels. 


In  the  presence  of  pharmacological  blockage  of  skeletal  muscle  contraction, 
dog  forelimb  resistance  was  decreased  by  vibration  (Liedtke,  and  Schmid,  1969). 

The  magnitude  of  the  decrease  was  reduced  with  limb  devervation  but  not  eliminated. 

There  are  no  previous  reports  of  the  detailed  cardiovascular  changes  pro¬ 
duced  by  vibration  of  unanesthetized  animals. 

In  the  human  subjects  studied  by  Lamb  and  Tenney  (Lamb,  and  Tenney ,1966) , 
heart  rate  and  blood  pressure  were  unchanged.  Hood,  et  al.,  1966,  described  in¬ 
creases  in  mean  arterial  blood  pressure,  heart  rate,  and  cardiac  output  in  tightly 
restrained,  semisupine  humans  exposed  to  7  minutes  of  wholebody,  x-axis  vibration. 
In  the  tightly  restrained,  upright  human  subjects  studied  by  Clark,  et  al.,  1967, 
mean  arterial  blood  pressure  and  heart  rate  increased  slightly,  and  cardiac  index 
approximately  doubled. 

The  metabolic-hormonal  effects  of  vibration  have  been  similarly  variable. 

Oxygen  consumption  has  been  found  to  increase,  decrease,  or  remain  unchanged 
in  both  anesthetized  dogs  and  humans  (Dixon,  et  al.,  1961;  Lamb  and  Tenney, 

1966-  Hoover  and  Ashe,  1962;  Hutt,  Horvath,  and  Spurr,  1958;  Duffner,  Hamilton 
and  Schmitz,  1962;  Gaeuman,  Hoover  and  Ashe,  1962;  Young,  et  al.,  1963;  Hood 
and  Higgins,  1965;  Hood,  et  al.  1966).  Hood  and  Higgins,  1965,  found  decreased 
oxygen  consumption  at  6  and  10  Hz  and  0.3  g  acceleration  amplitude;  at  1.3  g 
it  increased  out  of  proportion  to  the  increase  in  cardiac  output.  Levels  of 
17-hydrocortosone  (17-0H-CS),  epinephrine,  and  total  catecholamines  were  ele¬ 
vated  in  supine  dogs  exposed  to  prolonged  periods  of  whole  body,  Z-axis  vibra¬ 
tion.  The  effects  decreased  with  increasing  displacement  and  acceleration  am¬ 
plitude  and  were  greater  in  awake  than  in  anesthetized  dogs  (Black,  1965). 

Plasma  and  urine  levels  of  17-0H-CS  decreased  in  humans  exposed  to  three  periods 
of  three  minutes  of  vibration  at  selected  frequencies  from  1  to  20  Hz  and  accele¬ 
ration  amplitudes  of  0.75  tc  3  G;  even  though  decreased,  values  remained  within 
normal  1  imita  (Hofflna;  Kanda,  and  Watanabe,  1971). 
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Although  the  trend  of  these  data  bears  a  superficial  resemblance  to  the 
changes  occurring  with  exercise,  the  responses  show  great  variation  and  differ 
in  important  ways  from  those  usually  produced  with  exercise.  Analysis  of  the 
effects  of  vibration  has  been  complicated  by  the  limited  combinations  of  fre¬ 
quency  and  amplitudes  of  displacement  and  acceleration  used  in  previous  studies. 
The  purpose  of  this  report  is  to  clarify  the  previous  observations  on  the  car¬ 
diovascular  effects  of  vibration  by  summarizing  results  from  an  extended  series 
of  tests  over  a  wide  range  of  sinusoidal  vibration  combinations  applied  to  both 
awake  and  anesthetized  dogs. 

B.  Experimental  Protocol 

Response  of  acutely  and  chronically  instrumented  dogs,  primates  and  pigs 
have  been  studied.  Special  emphasis  has  been  placed  on  the  chronically  implant¬ 
ed  preparation  which  permits  repeated  awake  and  anesthetized  tests  on  the  same 
animal.  Additional  advantages  of  this  preparation  include  the  closed  chest  pro¬ 
tocol  and  more  stable  transducer  characteristics,  especially  important  for  the 
vibration  environment. 

Animals  are  restrained  on  the  hydraulic  table  with  the  spine  vertical,  and 
vibrated  along  that  axis.  Short  term  tests  consist  of  sequential  exposures  to 
sinusoidal  vibration  of  2-12  Hz  and  acceleration  amplitudes  ranging  from  0.5  g 
to  3.0  g  for  30  sec.,  2  min.,  and  5  min.  intervals.  Long  term  exposures  In¬ 
clude  durations  varying  from  30  min  to  6  hrs .  at  a  selected  vibration  frequency 
and  constant  acceleration  amplitude.  The  highest  g-level  for  the  awake  animal 
is  1.5  g,  and  for  the  anesthetized  animal  is  3.0  g.  The  cardiovascular  variables 
include  heart  rate,  flow  velocity  from  the  aorta  and  peripheral  arteries,  arterial 
pressures  from  acceleration -insensitive  intraaortic  and  intraventricular  trans¬ 
ducers,  Oj  consumption  and  body  temp,  and  blood  chemistries  such  as  blood  gases, 
cortisol,  free  fatty  acids,  and  glucose,  measured  on  samples  withdrawn  from  in¬ 
dwelling  atrial  and  arterial  catheters.  Tne  design  of  the  vibration  exciters 


permits  recording  of  force  transmitted  between  subject  and  table,  and  also  the 
vibration  exciter  velocity.  From  these  variables  and  their  phase  relationship, 
complex  mechanical  velocity  impedance  can  be  determined.  In  addition,  accelero¬ 
meters  implanted  on  large  body  organs  provide  information  on  the  displacement 
of  these  organs  relativ  he  animal  torso  and  the  input  vibration.  A  Raytheon 

704  data  acquisition  syste  _  ing  programmed  for  off-line  and  on-line  analysis 
of  the  massive  amount  of  cardiovascular  and  mechanical  data. 

1.  Animal  Preparation:  Two  basic  surgical  procedures  have  been  imple¬ 
mented.  The  most  consistent  results  have  been  obtained  with  the  intra-thoraeic 
implantations  involving  the  heart  and  large  vessels.  Instrumentation  implanted 
may  include  aortic  flow  probe,  aortic  and  left  ventricular  pressure  transducers, 
right  and  left  atrial  and  subclavian-aortic  cannulae,  pacemaker  on  right  atrium, 
temperature  probe  within  the  chest,  and  an  accelerometer  (single  axis)  on  the 
heart  or  proximal  blood  vessel.  Abdominal  and  neck  entries  have  been  employed  to 
a  lesser  extent  for  the  measurement  of  peripheral  cardiovascular  parameters. 

The  choice  of  transducers  to  be  implanted  in  both  preparations  is  flexible 
although  experience  has  shown  that  standardizing  the  surgical  procedure  and 
setting  a  limit  on  the  number  of  devices  implanted  is  necessary.  As  part  of 
the  program  design,  the  chronic  preparations  have  allowed  a  number  of  differing 
experimental  protocols  to  be  followed.  Such  experiments  include  vertical  z-axis 
(and  more  recently,  horizontal  z-axis)  vibration,  hypovolemic  shock,  barorecep- 
tor  response,  and  animal  performance  tasks.  A  combined  chronic-closed  chest 
and  acute-peripheral  vessel  preparation  has  also  been  used  successfully. 

A  healthy  animal  and  reliable  implanted  instrumentation  make  a  successful 
animal  preparation  and  will  yield  usable  physiological  data.  Individual  ani¬ 
mals  are  tested  for  acceptance,  and  these  guidelines,  combined  with  previously 
published  normals,  provide  control  values  for  post-operative  care.  Along  with 


Che  hematological  variables  examined  in  the  implanted  animal,  metabolic 
and  hormonal  determinations  now  include  blood  gases,  cortisol,  free  fatty 
acids,  and  glucose. 

Thorough  evaluation,  regular  testing,  and  redesign  of  transducers  and 
materials  have  become  an  integral  part  of  our  implantation  schedule.  In 
particular,  pressure  transducers  and  flow  probes  (as  well  as  flowmeter  in¬ 
strumentation)  have  been  subjected  to  these  continuous  testa.  Experimental 
access  to  the  implanted  transducers  and  catheters  stored  in  a  fabric  pouch 
of  the  chronic  animal  is  quick  and  non- traumatic .  Details  of  the  animal  pre¬ 
paration  are  presented  in  Appendix  A. 

2.  Experimental  Procedure:  On  the  day  of  the  experiment  the  dog  is 
walked  to  the  laboratory  and  placed  in  r.  nylon  mesh  sling.  The  monkeys  are 
transported  to  the  vibration  room  in  an  open-back  restraint  chair.  While  the 
leads  are  removed  from  the  storage  pouch,  the  vibration  facilities  are  cali¬ 
brated.  If  all  the  transducers  and  cannulae  are  functioning,  the  animal  is 
then  placed  in  the  vibration  chair  and  firmly  restrained.  Acute  placement 
of  a  transducer,  if  necessary,  is  accomplished  at  this  time  using  a  local 
anesthetic  (Xylocaine)  and,  if  required,  morphine  I.M.  (1-3  mg/kg). 

Calibration  of  the  pressure  gauges  is  done  with  chronically  placed  cannu¬ 
lae  and  an  extravascular  pressure  gauge,  making  sure  that  measurements  do  noc 
include  an  additional  hydrostitic  pressure  head.  Some  of  the  analog  transducer 
signals  are  then  sent  to  the  computer  for  A/D  conversion. 

Following  transducer  calibration,  the  dog  is  placed  in  the  restraint 
chair  and  the  restraint  chair  is  then  attached  to  the  hydraulic  table.  Pre¬ 
vibration  control  and  post-vibration  recovery  periods  normally  last  for  10 
to  30  minutes.  Depending  on  the  vibration  protocol,  periodic  blood  sampling 
and  instrument  re-calibration  and  balancing  may  be  done  before,  during,  and 
after  the  experiment. 


29 


For  an  anesthetized  preparation,  the  animal  is  given  an  intramuscular 
injection  of  morphine  (3  mg/kg) ,  a  blood  sample  for  blood  gas  controls  is 
withdrawn,  and  after  30  minutes,  the  .  nimal  is  anesthetized  with  chloralose- 
urethane  given  I.V.  The  level  of  anesthesia  used  is  that  required  to  sup¬ 
press  the  voluntary  ejection  of  a  tracheotomy  tube.  The  animal  is  then  res- 
pirated  with  room  air  and  blood  gases  are  regularly  monitored  and  controlled 
with  the  respiratory  rate  and/or  volume. 

Following  the  experiment,  if  the  animal  is  to  be  used  for  further  study, 
the  lead  connectors  are  sealed,  cleaned,  and  returned  to  the  pouch  which  is 
closed  with  interrupted  sutures  or  steel  clips.  The  pouch  area  is  bandaged, 
a  nylon  jacket  is  placed  on  the  animal  (for  monkeys,  a  leather  jacket  is 
also  fitted),  and  the  animal  is  returned  to  the  animal  care  facility.  Further 
details  of  the  experimental  procedure  are  presented  in  Appendix  A. 

3.  Data  Acquisition  and  Analysis:  The  Raytheon  704  system  (Appendix  B) 
is  designed  to  fulfill  two  major  requirements  for  the  vibrational  6tress  ex¬ 
periments.  First,  a  real-time  output  of  various  physiological  variables  in 
terms  of  graphs  on  a  storage  oscilloscope  and  a  printed  record  on  a  teletype 
is  needed.  This  will  provide  feedback  to  the  experiment  controller  so  that 
the  protocol  may  be  varied  according  to  the  results  being  obtained.  Secondly, 
the  data  reduction  process  must  be  speeded  up  and  improved  in  accuracy.  The 

704  is  fast  enough  to  provide  several  real  time  outputs  and  at  the  same  time 

store  data  on  magnetic  tape  for  later  operator-controlled  analysis. 

The  two  program  systems  (Dye  Curve  and  Mean  Value)  which  have  been  imple¬ 
mented  to  date  are  described  in  detail  in  Appendix  B. 

The  more  ambitious  part  of  the  design  system  is  now  being  developed.  Out¬ 
puts  of  pulsatile  as  well  as  average  signals  will  be  provided  in  real  time. 

A  real  time  CRT  display  of  the  physiological  and  mechanical  variables  will  pro- 


vide  compression  of  the  data  and  feedback  to  the  experimenter.  Preliminary 
specifications  for  this  system  are  also  presented  in  Appendix  C. 

C.  Results 
1.  Experimental 

Physiological  Responses:  For  unanesthetized  dogs  receiving  short-term 
(iO  sec)  exposure  to  vertical,  whole-body  vibration,  the  major  finding  is  the 
overall  linear  relationship  between  mean  flow  in  the  ascending  aorta  (MAF, 
equal  to  cardiac  output  less  coronary  flow)  and  the  level  of  vibration  stress 
as  measured  by  the  ratio  of  peak  net  transmitted  force  to  body  weight  (PNF/BW) 
(Fig.  2).  For  10  experiments  on  6  dogs,  the  regression  equation  for  the  rela¬ 
tionship  between  percent  change  in  MAF  and  per  cent  change  in  PNF/BW  was  0.00 
+0.476  PNF/BW,  with  a  standard  error  of  the  estimate  of  ±26.71  and  a  R(correla- 
tion  coefficient)  of  0.70. 

At  a  PNF/BW  level  in  the  range  of  1.5  to  2.0  times  body  weight,  near  the 
maximum  values  applied  in  these  studies  and  in  the  region  of  resonance,  the 
mean  increase  in  MAF  was  1.6  times  the  control  value.  For  PNF/BW  levels  in  the 
neighborhood  of  0.8  times  body  weight,  a  few  of  the  MAF  values  were  unchanged 
or  decreased  from  control  levels.  Extreme  values  for  the  change  in  cardiac  out¬ 
put  were  +150  percent  and  -30  percent  of  the  control  value  taken  Immediately 
preceding  the  exposure. 

The  relationship  between  MAJ  and  PNF/BW  was  much  more  closely  correlated 
in  individual  dogs  (Figures  3  a  and  4a).  Also,  mean  aortic  blood  pressure  was 
generally  unaltered  (Figure  3a)  hence  total  peripheral  resistance  decreased 
(conductance  increased).  The  relative  contributions  of  heart  rate  and  stroke 
volume  to  the  increases  in  MAF  exhibited  three  distinct  classifications.  In 
three  of  the  10  experiments  (2  animals)  MAF  was  increased  by  increased  heart 
rate,  with  minimal  or  no  change  in  stroke  volume  (Category  1,  an  example  of  one 
animal  io  shown  in  Figures  3a,  b,  and  c) .  In  4  of  the  10  experiments  (3  animals) 
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MAF  was  altered  primarily  by  changing  stroke  volume,  with  little  or  no*thange  in 
heart  rate  (Category  2,  an  example  of  one  animal  is  shown  in  Figures  ^a,  b,  nnd  c) . 
Approximately  equal  contributions  from  changing  heart  rate  and  stroke  volume  were  pre¬ 
sent  in  3  of  the  10  experiments  (3  animals):  in  this  category,  both  variables  were 
linearly  related  to  PNF/BW  (Category  3,  an  example  is  shown  in  Figures  5a,  b,  and  c) . 

Domination  of  heart  rate  as  the  principal  mechanism  for  increasing  MAF  ap¬ 
peared  to  be  a  function  of  the  heart  rate  level  at  the  initiation  of  vibration 
exposure  of  Category  1  and  2  animals.  In  those  animals  with  initial  heart  rates 
under  150  beats/min  (BPM) ,  altered  heart  rate  was  the  major  response  (Category  1) ; 
for  those  animals  with  initial  heart  rates  above  150  BPM,  altered  stroke  volume 
was  dominant  (Category  2).  Initial  heart  rates  were  not  the  significant  factor 
for  those  animals  changing  both  heart  rate  and  stroke  volume  (Category  3).  Ir. 
this  category  initial  heart  rates  ranged  from  78  to  204  BPM. 

The  general  pattern  in  the  anesthetized  dogs  was  similar,  but  grouped  res- 
sponses  showed  much  more  scatter  (Figure  6).  The  regression  equation  for  the 
relationship  between  percent  change  in  MAF  and  percent  change  in  PNF/BW  was  0.40 
+0.15  PNF/BW,  with  a  standard  error  of  the  estimate  of  ±7.8  and  a  R  of  0.27.  In¬ 
dividual  correlations  between  MAF  and  PNF/BW  were  again  high. 

In  contrast  to  the  three  categories  of  MAF  alteration  in  the  awake  animalB , 
the  anesthetized  animals  increased  MAF  through  only  two  mechanisms,  some  exhibiting 
a  modest  increase  in  heart  rate  and  others  a  more  significant  increase  in  stroke 
volume . 

For  the  3  of  6  animals  altering  MaF  by  a  change  in  heart  rate,  the  mean  changes 
in  the  anesthetized  animals  were  95%  lower  than  those  in  the  unanesthetized  animals 
(Fig.  7a,  b  and  c) .  For  the  4  of  6  anesthetized  animals  altering  MAF  primarily 
through  changes  in  stroke  volume,  the  differences  were  not  statistically  significant 
from  those  of  the  awake  animals  (Fig.  8a, b,  and  c).  When  examined  individually, 
the  slope  of  MAF  vs  PNF/BW  varied  directly  with  the  initial  stroke  volume  of  the 


anesthetized  dogs. 


While  the  stress  level  as  measured  from  peak  net  force  is  of  major  impor¬ 
tance,  it  is  not  the  exclusive  determinant  of  the  response,  for  the  relation¬ 
ship  between  the  percent  change  in  MAT  and  the  log  of  the  ratio  of  mean  heart 
rate  to  vibration  frequency  is  also  linear  (Fig.  9). 

Mechanical  Responses:  A  portion  of  the  research  during  the  past  year  in¬ 
volved  a  coordinated  effort  between  the  AFOSR  and  the  AMRL  contracts.  Work 
performed  on  the  latter  was  primarily  directed  toward  investigating  the  effects 
of  vibration  upon  the  tracking  performance  of  Rhesus  monkeys.  Although  the  ulti¬ 
mate  goals  of  the  two  contracts  were  different,  each  involved  vibration  testing 
of  animals  during  which  mechanical  and  physiological  measurements  were  made.  A 
part  of  the  effort  on  each  contract  thus  complemented  the  other.  One  such  com¬ 
plementary  phase  involved  the  development  of  a  two  mass,  single-degree-of-free- 
dom  analytical  model  to  predict  the  mechanical  Impedance  response  of  the  sitting 
Rhesus  monkey  to  2  to  30  Hz  vibration  at  acceleration  amplitudes  from  0.5  to  l.Og, 
Figure  10  illustrates  the  sequential  development  of  the  model,  and  figure  11 
presents  the  final  model  with  linear  approximations  for  frequency  dependent  whole- 
body  coefficients  of  elasticity  and  damping  normalized  by  animal  body  weight.  Fig 
ure  12  contains  gr«ons  of  the  elastic  and  damping  coefficients  versus  frequency. 
The  values  were  obtained  for  each  test  at  each  frequency  by  requiring  the  Impedanc 
of  the  model  to  match  that  of  the  primate.  A  part  of  the  model  Is  an  empirical 
equation  defining,  as  a  function  of  vibration  frequency  and  intensity,  the  ratio 
of  non-reactive  to  reactive  animal  mass.  After  normalizing  for  body  mass  and 
averaging,  functions  for  the  elastic  and  damping  coefficients  were  developed  such 
that  the  impedance  response  of  the  model  closely  matched  the  experimental  results 
throughout  the  2  to  30  Hz  range  at  0.5  G  (see  Figure  13)  and  also  at  1.0  G  (see 
Figure  14).  A  detailed  description  of  this  study  can  be  found  in  the  appendix 
in  the  paper  "A  Model  to  Predict  the  Mechanical  Impedance  of  the  Sitting  primate 


During  Sinusoidal  Vibration"  by  R.  G.  Edwards  and  J.  F.  Laffert'  (Appendix  D) . 

The  amount  of  energy  dissipated  by  viscous  damping  in  a  vibrated  animal  can 
be  derived  from  the  experimentally  measured  impedance  data.  Figure  15,  for  a 
sitting  Rhesus  monkey,  contains  graphs  of  dissipated  energy  versus  time  of  vibra¬ 
tion  exposure  for  three  different  frequencies,  i.e.  6,  12,  and  20  Hz,  each  at 
0.5  g  acceleration  amplitude.  Although  the  data  from  the  12  and  20  Hz  exposures 
did  not  indicate  any  significant  temporal  trends,  the  6  Hz  exposure  did  indicate 
a  reduction  in  dissipated  energy  with  increasing  time  of  vibration  exposure.  Dur¬ 
ing  ten  1,5  g  tests  of  two  primates  the  acceleration  at  the  top  of  the  skull  was 
recorded.  This  was  accomplished  by  securely  taping  to  the  top  of  the  animal’s 
head  a  minute  piezo-resistive  accelerometer.  The  sensitive  axis  of  the  accelero¬ 
meter  was  approximately  along  the  axis  of  vibration.  This  measurement  enabled 
the  calculation  of  the  ratio  of  head  to  vibration  exciter  acceleration,  i.e.  the 
head  to  vibration  exciter  transmissibility .  The  purpose  here  was  to  establish, 
as  a  function  of  frequency,  how  much  of  the  input  acceleration  was  actually  trans¬ 
mitted  to  the  subject's  head.  Figure  16  is  a  graph  average  head  to  vibration 
exciter  transmissibility  versus  frequency.  At  3  Hz  the  1.5  g  applied  to  the  sub¬ 
ject  resulted  in  an  average  peak  acceleration  at  the  head  of  approximately  1.6 
times  greater  than  the  input  value,  i.e.  of  about  2.4  g.  From  3  to  30  Hz  the 
transmissibility  decreased  with  increasing  frequency.  The  acceleration  of  th> 
head  at  frequencies  less  than  12  Hz  was  amplified  relative  to  the  input  accelera¬ 
tion,  i.e.  the  transmissibility  was  greeter  than  1.0,  while  from  12  to  30  Hz  it 
was  attenuated,  i.e.  it  was  less  than  1.0.  At  30  Hz  only  20 %  of  the  input  accel¬ 
eration  was  transmitted  to  the  head. 

As  an  example  of  the  impedance  data  of  the  dog  tests  consider  the  impedance 
modulus  and  phase  angle  versus  vibration  frequency  graphs  contained  in  Figure  17. 
Two  curves  appear  on  each  graph;  one  for  a  1  g  test  at  frequencies  from  2  to  20 
ilz  with  the  animal  unaneaUici.i^ed ,  and  the  other  for  a  similar  vibration  protocol 
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but  with  the  subject  anesthetized  with  chloralose  -urethane.  Relative  to  the 
unanesthetized  data,  it  is  apparent  from  these  graphs  that  anesthetizing  the 
animal  resulted  in  a  2  Hz  decrease  in  the  primary  whole-body  impedance  resonant 
frequency,  i.e.  from  6  to  A  Hz,  and,  in  general,  a  more  "mass-like"  response 
(compare  to  the  inert  mass,  or  "mu"  impedance  response).  At  frequencies  from 
7  to  12  Hz  the  impedance  phase  angle  for  the  drugged  animal  is  approximately  20° 
greater  than  that  for  the  unanesthetized  condition.  Since  the  phase  angle  for  an 
inert  mass  is  9C P,  and  those  for  a  massless  damper  and  massless  spring  are,  respec¬ 
tively,  0°  and  -90°,  both  the  impedance  modulus  and  phase  angle  data  reflect  a  de¬ 
creased  muscle  tonus,  i.e.  a  decreased  whole-body  elastic  coefficient,  for  the 
anesthetized  (relative  to  the  unanesthetized)  dog.  This  kind  of  impedance  data 
interpretation  is  particularly  useful  where  applied  In  a  model  from  which  whole 
body  coefficients  of  elasticity  and  damping  can  be  extracted. 

Similar  to  the  aforementioned  transmissibility  measurement  on  the  Rhesus 
monkey,  it  was  desired  to  gain  information  on  the  transmissibility  of  certain 
internal  organs  and  vessels  for  the  vibrated  dog,  i.e.  the  amount  of  organ  and/or 
vessel  movement  relative  to  that  applied  by  the  vibration  exciter.  A  minute 
Konigsberg  piezo-resistive  accelerometer  (Model  A2)  was  securely  sutured  to  an  im¬ 
planted  aortic  flow  transducer  on  one  animal  for  the  purpose  of  gaining  informa¬ 
tion  as  to  the  amount  of  input  vibration  that  is  transmitted  to  the  aorta.  The 
upper  portion  of  Figure  18  la  a  graph  of  the  ratio  of  acceleration  at  th>  aortic 
transducer  legation  to  that  of  the  vibration  exciter,  i.e.  to  that  of  the  input 
vibration,  versus  vibration  frequency  with  0.5,  1.0,  and  1.5  g  as  parameters.  The 
lower  portion  of  Figure  18  is  a  corresponding  graph  of  the  phase  relationship 
between  the  aortic  transducer  acceleration  and  that  of  the  vibration  exciter. 

The  elastic  nature  of  the  thoracoabdominal  system  is  clearly  evident  from  these 
graphs.  Although  some  dif ference.  do  exist  as  a  function  of  the  input  accelera¬ 
tion  amplitude,  a  first  approximation  could  well  consider  the  response  linear  in 
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the  2  -  30  Hz  range  for  acceleration  amplitudes  from  0.3  to  1.5  g.  The  magni¬ 
tude  curves  clearly  show  a  resonance  at  4  Hz.  A  secondary  resonance  occurred 
at  9  Hz  during  the  0.5  and  1.0  g  tests,  but  was  not  present  in  the  1.5  g  test. 
During  vibration  at  the  4  Hz  primary  resonant  frequency  a  transmissibility  of 
from  2.3  to  2.4  was  recorded.  For  example,  during  5  Hz  vibration  at  1.5  g  an 
acceleration  of  approximately  3.6  g  was  recorded  at  the  aorta.  Beyond  approxi¬ 
mately  11  Hz  the  transmissibility  drops  below  1.0,  i.e.  theic  ia  less  accelera¬ 
tion  transmitted  to  the  aorta  than  is  applied  to  the  animal.  At  30  Hz  only 
about  15%  of  the  applied  acceleration  is  measured  at  the  aorta.  The  almost 
linear  increase  in  phase  angle  with  increasing  frequency  is  indicative  of  an 
elastic-type  response.  The  graphs  of  figure  18  illustrate  the  frequency  depen¬ 
dence  of  the  aortic  transmissibility,  i.e.  of  the  amount  of  acceleration  trans¬ 
mitted  to  the  aorta  as  compared  to  that  input  by  the  vibration  exciter.  These 
data  represent  the  transmissibility  of  the  aorta  with  electromagnetic  flow  trans¬ 
ducers  and  accelerometer  attached.  The  added  masB  of  these  transducers  result 
in  ttariuiulssibility  values  different  from  those  which  would  be  measured  from  a 
"mass-less"  accelerometer;  however,  these  reported  values  are  thought  to  be  rea¬ 
sonably  indicative.  These  graphs  illustrate  two  very  pertinent  facts  with 
respect  to  studies  such  as  the  current  one;  namely  (1)  the  magnitude  of  any  given 
applied  whole  body  vibration  transmitted  to  the  aorta  is  extremely  frequency 
dependent,  and  (2)  there  is  a  resonance  at  approximately  4  Hz,  at  which  frequency 
the  applied  vibration  amplitude  is  intensified  approximately  2.4  times  when 
measured  at  the  aorta. 

2.  Computer  Modeling 

Cardiovascular  modeling  has  progressed  in  two  directions.  The  first  of  these 
is  a  closed-loop  electrical  analog  model  of  the  hydraulic  and  mechanical  aspects 
of  the  CV  system,  which  is  an  extension  of  an  earlier  open-loop  model.  These 
effects  of  vibration  are  simulated  by  pressure  source  terms  generated  from  values 
of  the  vibration  g-level  or  resulting  transmitted  force.  These  pressure  distur- 
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bances  propagate  and  reflect,  adding  and  subtracting  with  the  pressure  generated 
by  the  action  of  the  heart,  and  as  a  result  produce  changes  in  flow.  The  closed- 
loop  model  contains  piecewise-linear  representations  of  venous  resistance  and 
compliance  which  were  included  to  accurately  account  for  blood-volume  shifts 
occurring  during  postural  changes  and  application  of  time-varying  whole-body 
acceleration.  The  possible  enhancement  of  cardiac  output  by  the  application  of 
whole  body  acceleration  synchronously  with  the  heart  cycle  was  investigated  by 
R.  L.  Starnes  in  an  M.S.  thesis  study  (see  Appendix  E) .  He  found  that  simulated 
cardiac  output  could  be  enhanced  by  only  a  few  percent  when  the  acceleratory  fre¬ 
quency  was  the  same  as  the  heart  frequency  (1.25  Hz).  This  study  has  recently 
been  extended  to  acceleratory  frequencies  at  whole  multiples  of  the  heart  fre¬ 
quency.  It  has  been  found  that  an  acceleratory  forcing  function  of  3.75  Hz  and 
4  G  peak  amplitude  enhances  cardiac  output  by  about  20  percent. 

A  preliminary  study  has  begun  in  which  experimental  data  of  HR,  SV  and  PNF 
are  used  as  input  parameters  to  the  model.  The  value  of  PNF/BW  vs  frequency  from 
the  animal  experiments  is  used  to  establish  the  values  of  the  pressure  source  terms 
for  the  different  simulated  vibration  frequencies.  Heart  rate  was  varied  in  the 
model  as  observed  in  the  animal  experiment.  The  computer  response  of  MAF  vs 
PNF/BW  was  similar  to  those  of  the  animal  experiments,  with  the  actual  percent 
change  values  close  to  those  found  for  the  group  of  anesthetized  dogs  changing 
MAF  with  minimal  change  in  heart  rate  and  little  change  in  stroke  volume,  e.g. 
maximum  percent  change  in  MAF  was  approximately  25%. 

The  second  direction  of  progress  is  the  development  of  a  digital  model  anal¬ 
ogous  to  the  analog  model.  Two  primary  reasons  for  moving  in  this  direction  are: 
(i)  the  digital  model  can  be  expanded  if  desired,  to  permit  more  detailed  regional 
investigations.  The  analog  model  has  been  expanded  to  the  point  where  present  com¬ 
puter  capacity  is  exhausted.  (2)  Nonlinear  baroreceptor  and  other  CNS  control  phe¬ 
nomena  can  be  more  easily  modeled  with  the  digital  approach.  The  digital  model  has 
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successfully  simulated  cardiovascular  performance  in  the  prone  position  (no 
acceleratory  stress)  and  is  presently  being  adapted  to  account  for  postural 
changes  and  time-varying  applied  accelerations. 

D.  Discussion 

For  the  range  of  vibration  stress  imposed  on  the  dogs  in  this  study  in 
the  vertical  (z-axis)  direction,  the  effect  of  vibration  on  mean  aortic  flow 
(MAF)  was  found  to  depend  primarily  on  the  net  force  absorbed  by  the  whole 
body  (peak  net  transmitted  force,  PNF).  Other  less  important,  but  significant, 
factors  are  initial  or  control  heart  rate  and  stroke  volume,  vibration  're- 
quency,  body  weight,  and  duration  of  exposure.  For  brief  exposures  at  PNF 
levels  quadruple  the  body  weight,  MAF  may  be  expected  to  reach  levels  nearly 
double  those  present  immediately  prior  to  the  exposure.  For  PNF  levels  less 
than  body  weight,  variation  in  response  may  include  a  few  instances  of  unchanged 
or  decreased  MAF. 

Detailed  comparison  of  these  findings  with  the  results  of  previous  investi¬ 
gations  is  difficult  since  frequency,  displacement,  and/or  acceleration  were  the 
only  vibration  parameters  available  in  earlier  data.  In  the  present  study,  test¬ 
ing  of  these  three  parameters  generally  resulted  in  much  greater  variability  with 
poorly  developed  trends.  Thus  the  major  determinant  of  the  response  appears  to 
be  the  net  effect  of  the  interaction  of  all  three  components  as  reflected  in  PNF. 

The  general  tendency  has  been  to  equate  the  effects  of  vibration  with  those 
occurring  in  mild  to  moderate  exercise.  Our  results  support  this  comparison  with 
regard  to  the  extent  of  the  change  in  MAF.  For  the  vibration  case,  however,  it 
is  unlikely  that  the  magnitude  of  the  increase  in  MAF,  and  the  interaction  between 
heart  rate  and  stroke  volume,  are  metabolically  dependent  in  the  manner  observed 
with  exercise.  For  instance,  heart  rate  and  cardiac  output  with  exercise  in  the 
dog  are  linearly  related  to  oxygen  consumption;  this  is  probably  not  the  case  for 
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vibration  stress  and  oxygen  consumption  are  probably  related  only  to  the  extent 
that  the  imposed  vibration  leads  to  muscle  contraction,  especially  for  postural 
adjustments  (Duffner,  Hamilton,  and  Schmitz,  1962). 

The  extent  of  muscle  contraction  is  somewhat  difficult  to  estimate.  The 
anesthetized  animals  of  Hood  and  Higgins,  1965,  responded  to  vibration  with  some 
increase  in  oxygen  consumption  which  was  reduced  by  curarization.  However,  in 
related  studies  in  our  laboratory  of  awake  Macaca  mulatta  exposed  to  force  levels 
similar  to  those  used  with  the  dogs,  oxygen  consumption  changed  with  varying  per¬ 
formance  demands  but  wa9  not  dependent  on  PNF/BW.  Results  from  selected  measure¬ 
ments  of  arterial  and  venous  pH,  pO^,  pCO^,  0^  and  CO^  content,  and  nonesterif ied 
fatty  acids  in  the  monkeys  and  in  the  dogs  in  the  present  study  do  not  show  any 
consistent  dependence  on  PNF  or  other  vibration  parameters  at  the  levels  used. 

Mediation  of  the  response  through  cyclic  muscle  contraction  induced  by 
variation  in  the  direction  of  the  I‘NF  is  unlikely,  since  cardiac  output  did  not 
follow  sinusoidal  exercise  above  0.2  Hz  (Ashkar,  1972):  a  rate  much  lower  than 
in  our  studies. 

Neural  mediation  appears  to  be  a  far  more  important  determinant  of  vibra¬ 
tion  responses  than  direct  metabolic  effects.  The  effect  on  MAF  was  greatly 
attenuated  in  the  group  of  anesthetized  dogs  with  limited  change  in  heart  rate. 
The  effect  in  the  unanesthetized  animals  i9  analogous  to  a  sustained  state  of 
alertness  or  alarm.  In  the  group  of  dogs  in  which  the  level  of  anesthesia 
limited  the  increase  in  heart  rate,  the  change  in  NaF  was  significantly  dif¬ 
ferent  from  that  produced  in  unanesthetized  dogs  Typically  with  exercise,  the 
rapid  onset  and  offset  of  the  heart  rate  response  is  dependent  on  neural  path¬ 
ways  (Rushmer,  Donald  and  Shepherd);  similar  effects  are  present  with  vibration. 
Mediation  through  the  carotid  sinus  is  unlikely  since  the  frequency  content 
added  to  the  phasic  pressure  waveform  is  relatively  high  (Edwards  et  al ,  1972) 
and  the  increased  heart  race  in  the  unane6thetized  animals  suggests  absence  of 


carotid  sinus  suppression. 

The  principal  neural  channel  which  could  lead  to  the  observed  effects  is 
afferent  sensory  information  from  peripheral  mechanoreceptors .  Small  myelinated 
and  unmyelinated  fibers  have  been  shown  to  be  primary  afferent  pathways  in  car¬ 
diovascular  and  respiratory  reflex  responses  (McCloskey,  Matthews,  and  Mitchell, 
1972).  Vibration-sensitive  receptors  are  included  in  the  proprioceptors  of 
these  afferents.  The  primary  afferents  of  muscle  spindles  are  not  significantly 
involved  in  these  effects;  high  frequency  (100  to  300  ilz)  vibration  of  the  triceps 
sura  muscles  of  both  hind  limbs  of  decerebrate  or  anesthetized  cats  had  no  appre¬ 
ciable  influence  on  systemic  arterial  blood  pressure,  heart  rate,  or  respiratory 
rate  or  depth  even  though  some  reflex  muscle  contraction  was  produced  (McCloskey, 
Matthews,  and  Mitchell,  1972). 

The  importance  of  neural  pathways  in  whole-limb  vibration  is  supported  by 
the  results  of  Liedtke  and  Schmid,  1969.  Peripheral  vasodilation  in  the  intact 
limb  was  much  greater  than  in  the  limb  following  denervation,  and  little  altera¬ 
tion  followed  section  of  the  carotid  sinus  and  vagal  nerves. 

The  initial  value  and  slope  dependence  of  the  variables  is  most  likely  a 
secondary  phenomenon  in  relation  to  the  effects  of  vibration,  but  is  of  consid¬ 
erable  interest.  The  dependence  of  the  magnitude  of  a  physiological  response 
cr.  the  preceding  level  is  to  some  degree  inherent  in  the  range  and  type  of  ad¬ 
justments  available.  At  high  heart  rates,  further  excitation  cannot  lead  to  a 
very  large  response  in  cardiac  output  through  further  increases  in  heart  rate; 
therefore  increases  in  stroke  volume  would  be  expected.  The  validity  of  the 
"law  of  initial  values"  (Wilder,  1957)  (Lacey,  1967)  for  non-extreme  levels  has 
beer,  difficult  to  substantiate,  but  for  heart  rate,  strong  presumptive  evidence 
of  its  presence  and  mediation  by  central  neural  activity  has  been  obtained  in 
monkeys  (Snapper,  Kadden,  and  Schoenfeld,  1971).  The  fact  that  many  of  the  un¬ 
anesthetized  dogs  exhibited  a  combined  heart  rate  and  stroke  volume  response  is 


unexplained  with  respect  to  vibration,  but  is  certainly  an  option  available  to 
the  system.  Dependence  of  MAF  on  the  ratio  of  mean  heart  rate  to  vibration 
frequency  suggests  the  presence  of  effects  derived  from  the  fluid -mechanical  and 
time-dependent  system  properties.  As  previously  discussed,  in  the  absence  of 
synchronization  between  vibration  and  cardiac  cyles,  this  effect  is  relatively 
minor. 

In  summary,  the  peak  net  transmitted  force  has  been  shown  to  be  a  major 
variable  determining  the  cardiovascular  response  of  dogs  exposed  to  brief,  whole 
body,  sinusoidal  vibration  applied  in  the  spinal  axis.  The  effect  is  markedly 
decreased  when  the  heart  rate  response  is  limited  by  anesthesia. 

Thus  vibration  appears  to  be  a  stress  which  can  be  graded  to  produce  a 
quantifiable  state  of  alertness  resembling  exercise,  but  separating  the  neural 
and  metabolic  components.  Analysis  of  system  properties  based  on  this  approach 
offers  great  promise  for  increasing  our  understanding  of  physiological  regula¬ 
tion. 

E.  Recommendation^ 

1.  Studies  should  be  continued  to  evaluate  further  the  role  of  the  five 
mechanisms  involved  in  producing  changes  in  pressures  and  flows  in  the  cardio¬ 
vascular  system  exposed  to  vibration.  Special  emphasis  should  be  placed  on  ex¬ 
periments  which  can  separate  and  elucidate  the  role  of  the  various  mechanorecep- 
tors  -  CNS  pathways  from  the  rest  of  the  mechanism.  This  is  especially  important 
since  these  mechanisms  are  thought  to  produce  about  70X  of  the  changes  observed 
in  the  awake  animal. 

2.  More  experiments  in  which  0^  consumption  is  measured  on  animals  under 
vibration  should  be  conducted  in  order  to  confirm  the  analogy  of  vibration  and 
exercise . 

3.  Experiments  should  be  conducted  in  which  the  vibration  frequency  is 
extended  below  2  Hz.  This  can  be  accomplished  by  using  a  variable  radius  centri- 


fuge  to  produce  slowly  varying  sinusoidal  acceleration. 

These  experiments  could  provide  some  of  the  missing  information  which  now 
exists  between  sustained  acceleration  (0  freq.)  and  vibration  (above  1  Hz),  i.e. 
the  slowly  varying  time  dependent  acceleration  range  (0-1  Hz) .  These  experiments 
would  provide  continuity  in  describing  the  response  of  the  cardiovascular  system 
from  sustained  acceleration  to  high  frequency  vibration. 

4.  The  mechanical  impedance  model  should  be  combined  with  the  analog  corn- 
outer  model  of  the  hydraulic  aspects  of  the  cardiovascular  system  in  order  to 
define  more  realistically  the  vibration  forcing  function,  PNF,  in  the  computer 
studies . 

5.  The  analog  computer  model  should  be  expanded  to  include  mechanorecep- 
tor  response  and  the  digital  model  should  be  expanded  to  include  more  detailed 
sections  of  the  arterial  and  vencus  system.  These  models  are  excellent  guides 
in  planning  experiments  and  evaluating  the  resulting  data. 
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Fig.  1  Block  diagram  of  the  interaction  of  major  physiological  systems  to 
produce  changes  in  pressure  and  flow  resulting  from  a  vibration 
input  into  the  whole  system. 
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Fig.  2  Precentage  change  in  mean  aortic  flow  as  a  function  of 
the  ratio  of  peak  net  force  to  body  weight  in  the  awake 
animals . 
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Fig,  3b  Percentage  change  in  mean  heart  rate  as  a  function  of  the 
ratio  of  peak  net  force  to  body  weight  in  animal  1600. 
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Fig.  4a  Percentage  change  m  mean  aortic  flow  as  a  tunc  cion  of 
the  ratio  of  peak  net  force  to  body  weight  in  an  awake 
animal  (5559)  whose  changes  were  mainly  due  to  stroke 
volume . 
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Fig.  4b  Percentage  change  in  stroke  volume  as  a  function  of  the 
ratio  of  peak  net  force  to  body  weight  in  animal  5559. 
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Fig.  5a  Percentage  change  in  mean  aortic  flow  as  a  function  of  the 
ratio  of  peak  net  force  to  body  weight  in  an  awake  animal 
(556)  with  fairly  equal  changes  in  mean  heart  rate  and 
stroke  volume. 
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Fi*  5c  Percentage  change  in  strode  volume  as  a  function  cf  the 
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Fig.  6  Percentage  change  in  mean  aortic  flow  as  a  function  of  the 
ratio  of  peak  net  force  to  body  weight  in  the  anesthetized 
animals . 
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Fig.  7a  Percentage  change  In  mean  aorcic  flow  ac  2  function  of  the 
ratio  of  peak  net  force  to  body  weight  for  the  group  of 
awake  animals  whose  changes  were  mainly  due  to  heart  rate. 
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Fig.  7c  Comparison  of  the  percentage  change  in  mean  aortic  flow 
as  a  function  of  the  ratio  of  peak  net  force  to  body 
weight  for  the  anesthetized  vs.  unai.esthetized  animals 
whose  changes  were  due  mainly  to  heart  rate. 
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Fig.  8b  Percentage  change  in  mean  aortic  flow  aa  a  function  of  the 
ratio  of  peak  net  force  to  body  weight  in  the  group  of 
anesthetized  animals  whose  changes  were  due  mainly  to 
stroke  volume. 
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Fig.  13  Impedance  phase  angle  and  normalized  modulus  versus  frequency 

for  0.5*g  vibration  of  four  primates  compared  to  that  predicted 
by  the  model  of  figure  11. 
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Fig.  14  Impedance  phase  angle  and  normalized  modulus  vs.  frequency 
for  1.0  g  vibration  of  four  primates  compared  to  that 
predicted  by  the  model  of  Fig.  11, 
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Fig.  16  Head  transmissibility  vs.  vibration  frequency  for  Rhesus  monkey 
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Fig.  1/  Impedance  moduluo  and  phase  angle  versus  vibration  frequency 
for  anesthetized  and  unanesthetized  dog  vibrated  at  1.0  g. 
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Aortic  transmi8sibility ,  magnitude,  and  phase  angle,  vs. 
vibration  frequency  for  a  dog  vibrated  at  0.5,  1.0,  and 


III.  Appendices 


Appendix  A 


Experimental  Techniques 

I.  Animal  Care 

A.  Pre-Operative  Animal  Care:  Surgery  for  the  chronically  implanted 
instrumentation  required  for  vibration  studies  has  been  performed  in  dogs  and 
and  monkeys.  Mongrel  dogs  (10-20  kg.)  supplied  through  the  University  of 
Kentucky  Medical  Center  have  been  generally  satisfac wt-iy .  An  attempt  to 
use  purebred  foxhounds  for  more  uniformity  was  not  successful  -  the  mon¬ 
grels  were  generally  hotter  behaved  for  awake  studies  and  were  much  hard¬ 
ier  for  thoracic  surgery.  Rhes  s  monkeys  (Macaca  mulatta),  weighing  5-10 
kg.  ,  have  been  generally  healthy  as  supplied  by  the  vivarium  at  Wright  Pat¬ 
terson  AFB.  As  an  alternate  to  the  Rhesus,  pigtail  monkeys  (M.  Nemistrina) 
of  a  similar  size  purchased  from  Primate  Imports  (Long  Island,  N.Y.)  have 
undergone  surgical  experimentation  although  none  have  been  exposed  to  vibra¬ 
tion  . 

Successful  chronic  implantation  procedures  and  techniques  and  surgical 
procedures  used  in  cue  chronically  instrumented  animal  must  begin  with  a 
health\  animal.  Success  is  directly  proportional  to  the  standards  set.  Be¬ 
cause  some  of  our  investigations  include  unane^thetized  preparations  in  ar. 
environment  wnich  tends  to  be  stressful  to  the  animal,  behavior  patterns, 
such  as  nervousness,  must  be  discerned  before  choosing  a  surgical  subject. 

The  physiological  guidelines  for  determination  of  normality  include 
physical  examination,  qualitative  estimation  from  behavior  patterns,  and 
results  of  hematological,  bacterial  culture,  and  x-ray  examinations.  Therapy 
for  canine  intestinal  parasites  is  repeated  until  stools  are  negative.  For 
the  dogs  a  hematocrit  vali ^  of  4 f  and  white  cell  count  of  10, COO  or  below  is 
expected.  The  hematocrit  ^or  monkeys  should  ha  40  to  42  with  a  white  count 
of  less  than  15,000.  In  addition  to  these  determinations,  values  for  many 


commonly  used  variables  in  dogs  appear  in  Table  A-l.  Multiple  similar 
variables  for  monkeys  are  listed  in  the  report  by  McCutcheon,  et  al.  (to 
be  published). 

B.  Surgery-General  Procedure:  Inclusion  of  each  of  the  following  criteria 
involving  the  surgical  procedure  has  been  found  necessary  to  produce  a  success¬ 
ful  chronic  animal  preparation: 

1.  sterile  preparation 

2.  surgically  convenient  and  properly  prepared  implants 

3.  minimal  trauma  and  resultant  minimal  fibrosis  within  surgical  area 

4.  isolation  of  Implants  from  external  environment 

The  above  criteria  complement  each  other  -  details  of  each  and  specific  aspects 
of  surgery  are  found  in  the  sections  below.  Following  immediately  is  a  general 
overview  of  thj  surgery  performed  on  animals  which  have  met  pre-implantation 
standards . 

Surgery  is  performed  with  complete  sterile  precautions  using  the  facilities 
in  the  Experimental  Surgery  Laboratory  of  the  University  of  Kentucky  Medical 
Center.  The  principles  of  laboratory  animal  care  as  outlined  by  the  National 
Society  for  Medical  Research  were  rigorously  observed.  The  monkey  is  preme¬ 
dicated  with  sernylan  and  atropine,  induced  with  thiopental,  intubated,  and 
enesthetized  with  halothane.  For  the  dog,  the  procedure  is  identical  with 
the  exception  that  sernylan  and  atropine  are  not  given.  The  approach  to  the 
chest  is  through  the  left  fourth  intercostal  space  with  the  various  trans¬ 
ducers  placed  as  required.  Each  lead  is  passed  through  different  spots  using 
intercostal  spaces  other  than  the  one  for  entry.  The  rib  cage  is  closed  in 
the  usual  manner  and  the  leads  are  tunneled  under  muscle  layers  to  emerge  be¬ 
tween  muscle  and  skin  in  the  center  of  the  back,  between  the  scapulae.  The 

B-  2 


> 


skin  ia  undermined  and  freed  up  from  the  original  incision  to  the  point  where 
the  leads  emerge.  The  leads  are  gathered  together  and  placed  in  a  fabric  ve¬ 
lour  bag  of  our  own  design  (McCutcheon,  et  al.,  1973a)  and  left  under  the  un¬ 
broken  skin.  The  muscle  layers  and  skin  of  the  original  incision  are  closed 
in  the  usual  manner.  A  chest  tube  connected  to  a  water  bottle  suction  is  left 
in  for  twenty-four  to  thirty-six  hours. 

The  abdominal  approach  Is  through  a  mid-line  incision  in  the  linea  alba. 
Females  are  much  more  suitable  for  the  abdominal  procedure.  The  probe  leads 
are  placed  on  the  back  in  a  manner  similar  to  that  used  for  the  thoracic  im¬ 
plantation. 

The  fabric  velour  pouch  used  to  store  the  transducer  connectors  under 
the  skin  is  sewed  from  nylon  velour  fabric  (fig.  A-l).  More  recently,  a  dacron 
velour  material  has  been  tested.  The  pouch  is  initially  gas  sterilized,  and 
after  its  implantation  the  r  t-*  is  left  undisturbed  to  allow  connective  tissue 
growth  into  the  fabric. 

Approximately  8  to  IT  days  after  surgery,  the  level  of  recovery  is  adequate 
to  open  the  pouch  ("windowing"  process)  and  obtain  access  to  the  leads.  Unacr 
tranquilization  and  local  anesthesia,  an  appropriate  length  of  akin  over  the 
pouch  is  incised.  The  pouch  i6  opened  with  scissors  and  the  leads  removed  and 
tested  (fig.  A-2).  The  edges  of  the  skin  are  retracted  and  sewn  to  the  outer 
layer  of  the  oouch.  The  inner  fabric  layer  is  closed  with  suture  or  steel 
clips  (fig.  A-3). 

After  the  bag  has  been  opened,  the  area  is  kept  covered  with  a  dressing  and 
a  nylon  mesh  Jacket.  The  jacket  has  prevented  access  by  the  dogs  to  any  exposed 
connectors  or  to  the  bag  itself.  Additional  caution  is  exercised  with  the  mon¬ 
key  -  a  custom-fitted  jacket  of  porous  girdle  material  is  the  first  covering,  and 
a  second  jacket  of  leather  is  fitted  over  this.  Both  Jacket.''  are  closed  by  sew- 
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C.  Surgery  -  Instrument  Implantation  Schedule:  The  standard  intrathoracic 
instrumentation  in  the  dog  includes  the  following: 

1.  Flow  probe  on  ascending  aorta 

2.  Pressure  gauge  in  apex  of  left  ventricle 

3.  Pressure  gauge  in  thoracic  aorta  just  distal  to  arch 

4.  Pacemaker  ECG  leads  on  right  atrium 

5.  Cannulae  in  right  atrium,  left  atrium,  and  through  subclavian 
artery  into  aorta. 

Other  implanteu  devices  which  have  been  substituted  for  some  of  the  above 
include  an  intrathoracic  temperature  probe  (thermistor  bead)  and  single  axis 
accelerometer . 

The  standard  intrathoracic  instrumentation  in  the  monkey  includes  the 
following: 

1.  Flow  probe  on  ascending  aorta 

2.  Pressure  gauge  in: 

a.  Descending  thoracic  aorta  with  subclavian  arterial  cannulae  or 

b.  Apex  of  left  ventricle 

3.  Cannulae  in  right  atrium  and  left  atrium  (unless  subclavian  cannula 
present) 

4.  Thermistor  temperature  probe 

/■  few  carotid  Implants  have  succeeded  but  recording  from  this  vessel  is 
not  y?t  a  standard  procedure.  In  order  to  acquire  move  ayperience  with  im¬ 
plantation  techniques  in  the  Rhesus  and  Pigtail  monkeys,  several  animals 
were  implanted  with  right  and  left  atrial  cannulae  and  intrathoracic  tempera¬ 
ture  probe. 

Instrument  placement  in  the  abdomen  (dogs  only)  includes  the  following: 

1.  Flow  probe  and  occluder  on  renal  artery 

2.  Flow  probe  and  occluder  on  mesenteric  artery 
Flew  probe  and  occluder  or.  iliac  artery  . 
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D.  Post-Operative  Animal  Care:  Aside  from  the  procedures  involved 
with  "windowing"  the  implanted  pouch,  routine  post-operative  management 


includes  antibiotic  coverage,  blood  hematology  and  chemistry  tests,  dress¬ 
ing  changes,  exercise  and  training  (if  necessary),  x-rays,  and  maintenance 
of  transducer  lead  integrity  and  catheter  patency.  The  fact  that  these 
animals  have  undergone  major  surgery  and  need  most  if  not  all  of  the  pre¬ 
ceding  diagnostic  tests  and  care  cannot  be  overemphasized. 

The  post-operative  antibiotic  schedule  is  coordinated  through  veter- 
narians  at  Wright-Patterson  AFB  and  the  University.  Cultures  of  suspect 
areas  are  taken  and  upon  receiving  sensitivity  results,  corrective  medica¬ 
tion  is  administered.  Simultaneous  blood  hematological  (i.e.  hematocrit, 
white  bT^ood  cell  count,  differential  leukocyte  count,  platelet  estimation) 
and  chemical  diagnostics  (total  protein,  glucose,  electrolytes)  aid  deter¬ 
mination  of  the  animal's  condition  and  also  resultant  trends  following  anti¬ 
biotic  dosages.  Since  many  antibiotics  adversely  affect  renal  function, 
urine  chemistries  (glucose,  proLein,  bilirubin)  and  sediment  analysis  are 
monitored  and  trends  noted. 

In  the  dog,  the  original  chest  incision  site  has  occasionally  failed  tc 
heal  completely,  and  over  an  extended  implantation  time  the  tissue  actually 
begins  to  break  down.  The  source  of  the  infection  has  been  both  external  and 
internal,  the  latter  arising  from  one  of  the  lead  tracts  exiting  from  the  win¬ 
dowing  pouch.  Adequate  care,  including  fluid  drainage  if  necessary,  and  bac¬ 
terial  cultures  have  been  successful  courses  of  action.  Of  the  cultures  taken, 
most  have  shown  Staphylococcus  aureus  (coagulate  positive  and  negative)  while 
others  have  grown  S taphylococcus  epidermis.  Other  than  topical  ointments  or 
powders,  post-surgical  treatment  is  accomplished  routinely  using  Combiotic  (some 
times  in  combination  with  i.oridine)  for  7  Lo  14  days.  Later  cultures  having 
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the  above  mentl  ">ned  growths  have  usually  been  sensitive  to  and  treated  with 
Erythromycin  or  Keflin. 

Intensive  antibiotic  treatment  in  the  monkey  has  been  the  only  route  to 
curb  a  frequently  occurring  infection  usually  beginning  in  the  outer  tissue 
surrounding  the  fabric  pouch.  Common  species  found  have  been  primarily  Pseu¬ 
domonas  aeruginosa  with  others  being  Proteus  mirabilis,  Enterobacter  aerogenes, 
and  Escherichia  coli.  As  with  the  dogs,  Combiotle  is  also  given  after  surgery 
followed  by  treatment  with  other  antibiotics  if  necessary.  This  latter  course 
has  included  Gentamycin  and  sometimes  Carbenicillin,  Chloramphenical ,  Loridine, 
Ample ill in,  and  Keflin. 

The  animal  is  weighed  periodically  and  quality  and  quantity  of  intake  and 
output  are  observed.  If  nutrition  seems  inadequate,  special  efforts  are  made 
to  tempt  him  with  favorite  foods,  and  hand-feeding  is  performed  if  necessary. 
Blood  status  is  monitored  as  indicated,  and  injections  of  iron  and  vitamins 
arc  giver.. 

The  r,.al  use  of  x-rays  over  the  period  of  implantation  provides  a  valu¬ 
able  diagnostic  tool  not  only  in  determining  the  health  course  of  the  animal 
but  also  the  status  of  the  implanted  transducers.  Lung  inflammation  and  other 
abnormalities  within  the  chest  region  can  be  detected.  Often  if  a  transducer 
does  fail,  an  x-ray  will  bring  out  a  broken  wire  connection  or  kink  in  a  can¬ 
nula  and  corrective  measures,  if  oossiblc,  are  taken.  A  representative  family 
of  x-raya  taken  from  a  dog  and  monkey  are  shown  respectively  in  figures  A-4 
and  A-5.  In  these  post-operative  pictures  the  absence  of  congestion  and  abnor¬ 
malities  with  concurrent  blood  analysis  and  observations  are  used  to  determine 
the  qualification  of  each  subject  for  experimental  investigation. 

One  further  aspect  of  post-operative  care  which  is  very  important  to  our 
particular  protocol  is  that  the  animal  receive  some  exercise  and  c  form  of  adap- 


Cive  training.  Both  dogs  and  monkeys  are  normally  caged,  and  the  routine 
of  daily  or  nearly  daily  walks  or  tasks  Is  necessary  to  examine  fitness  for 
vibration,  speed  recovery,  and  adapt  the  animal  to  the  personnel  involved  with 
the  experiment.  Experiance  has  shown  that  prior  -waraness  of  (not  necessarily 
adaptation  to)  the  vertical  restraint  chair  (especially  dogs)  tends  to  calm 
the  animal  and  result  in  a  more  meaningful  experiment.  Adaptation  to  over- 
the-head  helmets  for  respiratory  studies  has  also  been  benificial.  Further 
improvements  in  post-operative  behavioral  care  are  being  worked  upon  and  the 
preceding  examples  only  illustrate  this  complex,  yet  often  overlooked,  area 
of  animal  care. 

Since  the  chronically  instrumented  animal  is  usually  sacrificed  under  in¬ 
vestigator  control,  an  autopsy  can  be  performed.  Much  vital  information  is 
available  from  these  examinations  which  have  been  routine  upon  the  death  of  a 
chronic  animal.  Tissue  specimens  are  sent  to  a  pathology  laboratory  and  the 
results  of  the  microscopic  examinations  have  guided  ub  during  program  develop¬ 
ment,  particularly  In  three  areas:  1)  tissue  Ingrowth  patterns  of  the  nylon 
(and  dacron)  velour  materials;  2)  specific  diagnosis  of  any  pathology,  If  pre¬ 
sent,  in  the  major  organa  of  the  animal  (i.e,  heart  wall)  muscle,  kidney, 
aortic  vessel  wall);  3)  Improvements  in  surgical  technique  or  Implanted  trans¬ 
ducer  design.  Specific  references  to  such  results  are  found  in  succeeding 
sections  of  this  report. 
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II .  Chronically  Implanted  Instrumentation 


The  selection  of  instrumentation  and/or  materials  for  implantation  1b 
one  of  the  more  difficult  aspects  of  developing  a  chronic  animal  preparation. 
Experience  is  a  lar^e  asset  and  careful  evaluation  of  any  device  prior  to 
and  during  chronic  implantation  is  of  the  utmost  importance.  The  failure 
of  one  component  in  such  a  preparation  can  be  costly  ir  terms  of  money,  man¬ 
power,  and  time.  Therefore,  in  our  laboratory  we  have  devoted  a  great  deal 
of  effort  to  these  concerns.  The  following  sections  describe  the  implanted 
instrumentation  uf.ed  in  our  laboratory,  and  the  procedures  we  follow.  In¬ 
cluded  with  each  description  are  design  considerations,  specific  surgical 
techniques,  and  results  of  evaluation  teats.  The  techniques  are  covered  in 
considerable  detail  since  they  are  extremely  important  for  overall  program 
success . 

A.  Blood  Flow 

1.  Implantation:  Blood  flow  measurement  is  one  of  the  most  valuable 
components  of  the  techniques  used  for  evaluation  of  the  circulatory  effects 
of  vibration.  For  cur  standard  chronic  animal  preparation,  electromagnetic 
cuff-type  flow  probes  are  placed  around  the  base  of  the  ascending  aorta.  At 
this  location,  a  satisfactory  physiological  zero  flow  baseline  is  present  - 
vespel  occlusion  is  unnecessary,  in  contrast  to  the  situation  for  peripheral 
vessels.  Unfortunately,  however,  the  wall  of  the  blood  vessel  under  the  probe 
at  this  site  is  subjected  to  large  stresses  due  to  the  rhythmical  wall  pulsa¬ 
tions  and  limited  anatomical  space  for  probe  placement.  The  effects  of  vibra¬ 
tion  (see  results  section  concerning  accelerometer  placement  on  the  flow  probe) 
compound  this  trauma  to  the  vessel. 

To  strengthen  the  flow  probe-aortic  tissue  interface,  end  to  minimize  the 
occurrence  of  vessel  rupture,  son>e  of  the  probes  have  been  narrowed  to  decrease 
the  probe-vessel  contact  area  and  allow  the  probe  cuff  to  surround  a  shorter 


section  of  aorta.  A  common  curtain  material  (Sears  "Entree  Tailored  Panel,  100% 
polyester  Sheer  Marquisette")  around  the  vesael  under  the  probe  aeema  to  minimize 
deterioration  of  the  aortic  wall.  Dacron  mesh  or  silastic  used  in  a  simmilar  man¬ 
ner  was  not  as  successful.  Even  with  these  methods  coupled  to  careful  surigcal 
implantation,  rupture  of  the  aorta,  generally  between  the  heart  and  probe,  accounts 
for  60%  of  the  deaths  in  our  chronically-prepared  animals  (note:  improving  the 
procedures  for  implantation  of  the  aortic  pressure  gauge  has  nearly  doubled  this 
percentage) . 

A  representative  pathology  report  describing  the  effects  on  an  intact  aortic 
wall  under  the  probe  after  56  days  (  cause  of  death  in  Dog  1600  was  rupture  of  tho¬ 
racic  aorta  at  pressure  gauge  site;  flow  probe:  In  vivo  metric  SL-1B,  02108,  20 
mm  I.D.  probe,  curtain  material)  was  as  follows: 

Sections  of  the  aorta  from  the  region  of  the  curtain  reveals 
the  material  on  the  outer  surface  of  the  aorta  and  each  of  the 
individual  threads  is  surrounded  by  fibrous  tissue.  The  aorta 
progressively  becomes  normal  toward  the  intimal  side.  The  outer 
half  of  the  aorta  in  this  area  reveals  separation  of  the  elastic 
fibrils  with  fibrosis  between  them.  In  the  most  outer  portion 
there  is  some  fragmentation  of  elastica.  The  curtain  appears  to 
be  firmly  embedded  with  this  adventitial  tissue. 

The  changes  are  more  extensive  when  the  cause  of  death  is  aortic  rupture 
proximal  to  the  flow  probe.  The  pathology  report  typical  for  this  latter  situ¬ 
ation  (Dog  822,  Zepeda  #923,  18  mm  I.D.,  Dacron  mesh,  duration  of  implantation 
20  days)  was  as  follows: 

Section  of  the  aortic  wall  under  the  flow  probe  near  the  site 
of  rupture  reveals  a  portion  of  the  aorta  with  a  zone  of  fibro¬ 
sis  in  the  outer  two-thirds  of  the  aorta,  merging  with  fibrosis 
within  the  adventitia  and  the  meshwork.  At  one  point  at  the  very 
edge  the  wall  appears  completely  necrotic.  This  probably  repre¬ 
sents  the  point  of  rupture. 

Alterations  in  the  aortic  wall  in  monkeys  are  similar  to  those  of  dogs. 

The  intact  vessel  wall  under  the  flow  probe  (Pigtail  monkey  PT4:  Zepeda  AA- 
#1142,  9  mm  I.D..  curtain  mater ial,  duration  51  days)  was  the  subject  of  the 
following  report: 


Section  of  the  aorta  shows  an  intact  lntima  and  an  intact  media. 

There  is  evidence  of  the  foreign  material  of  the  curtain  in  the 
adventitia  associated  with  fibrosis  and  some  degeneration  in 
the  very  outer  portions  of  the  media.  There  is  also  associated 
intimal  fibrosis. 

The  above  results  support  the  conclusion  that  a  reinforcing  material 
under  the  probe  probably  provides  a  useful  additional  Interface  between  the 
vessel  wall  and  probe. 

2.  Flowmeter  and  Flovprobe  Evaluation:  For  a  flow  probe  located  on  the 
ascending  aorta,  a  satisfactory  physiological  zero  flow  baseline  reference  is 
present.  This  baseline  can  drift,  however,  depending  upon  the  flow  probe  and/ 
or  the  flowmeter  to  which  the  probe  is  connected.  A  number  of  specific  causes 
have  been  detailed  (Gordon,  1971;  Cappelen,  1968;  Hagnestad,  1961,  Dobson, 
et  al.,  1966;  Spencer  and  Denison,  1959;  Fryer  and  Sandler,  1971)  in  the  litera¬ 
ture. 

Electromagnetic  flowmeters  undergo  continuing  evaluation  in  our  laboratory. 

Those  currently  in  U9e  are  clearly  not  optimum  instruments  (Biotronex  laboratory, 
Maryland,  Model  610  Pulsed  Logic  Flowmeter,  and  Zepeda  Instruments,  Seattle, 
Washington  Model  SWF-2  Square  Wave  Flowmeter).  For  instance,  the  Biotronex  de¬ 
vice  has  no  "zero"  capability.  The  Zepeda  flowmeter,  vith  careful  bench  cali¬ 
bration,  is  reported  to  provide  a  "magnet  zero"  within  5 X  of  the  true  blood  zero. 
Preliminary  bench  teste  indicate  that  for  moderate  flows  (]  L/min«f lov-5  L/min), 
the  magnet  zero  is  no  more  than  5%  from  actual  zero  baseline  flow.  However,  no 
Internal  calibration  for  checking  flowmeter  and  recorder  amplifier  gain  is  avail¬ 
able  on  the  Zepeda  instrument.  More  long-term  drift  studies  need  to  be  done  uo 
permit  adequate  comparisons  of  different  probes  over  a  period  of  time  or  the  nam: 
probe  over  different  time  intervals.  Evaluation  of  flow  meters  measuring  smaller 
peripheral  flows  with  chronically  implanted  probes  requires  vascular  occlusive  de¬ 
vices  more  satisfactory  than  any  developed  to  date  in  our  laboratory.  Further 
acute  animal  studies  are  now  underway  to  add  more  information  and  direct  our  efforts 
to  solve  this  problem  (cee  also  section  B) . 


Three  different  manufacturer a '  probes  have  been  Implanted  in  our  chronic 
animals.  In  the  dogs,  we  have  used  flow  probes  made  by  Zepeda  Instruments, 
Blotronex  Laboratory,  and  In  Vivo  Metric;  in  monkeys,  we  have  implanted  Zepeda 
and  Blotronex  transducers.  When  connected  to  the  EL  610  flowmeter,  all  the 
larger  probes  exhibited  essentially  the  same  output  characteristics  in  all  chron¬ 
ic  animals.  The  IVM  probes  tended  to  be  the  lea9t  stable  during  bench  calibra¬ 
tions  and  implantation,  probably  as  a  result  of  their  narrower  width  and  smaller 
electrode  area.  The  Zepeda  probe,  now  the  standard  in  our  laboratory,  is  quite 
stable  although  its  greater  width  contributes  to  some  of  the  trauma  seen  within 
the  wall  of  the  ascending  aorta.  Both  the  Blotronex  and  iVM  probes  used  have 
been  more  succeptible  to  internal  or  connector  failures  while  implanted  and  dur¬ 
ing  recycling  checks  between  lioplants. 

The  outer  diameters  of  the  ascending  aorta  of  the  monkeys  implanted  to  date 
ranged  from  7  to  12  mm.  In  this  range  the  flow  probe  most  used  has  been  the  axial 
lead,  symmetrical  core  Zepeda  probe.  A  minor  drawback  of  this  probe,  the  narrow 
elot  key  opening,  which  makes  surgical  manipulation  slighly  more  difficult,  is 
more  than  offset  by  its  stable  performance  and  durable  construction. 

A  still  smaller  range  of  probe  sizes  has  been  briefly  studied  during  abdo¬ 
minal  implantation.  Three  to  five  millimeter  lumen  diameter  probes  from  the  three 
manufacturers  listed  above  have  been  used.  An  evaluation  of  the  performance  of 
these  probes  will  be  available  following  more  experimentation  and  bench  calibra¬ 
tion  now  in  progress. 

The  technique  of  flew  probe  connector  placement  within  a  subcutaneous 
;:abric  pouch  places  additional  restraints  on  the  probes  selected.  Requiring  a 
larger  pouch  volume  and  window  slit  area,  the  bulky  button  or  bullet  connector, 
standard  wir.h  many  flow  probe3,  is  unsuitable  for  our  particular  chronic  prepara¬ 
tion.  The  Zepeda  connectors  have  a  base  of  shrink-tight  tubing  and  male  pins 
sealed  with  uhrlnk-tight  during  storage.  These  have  proven  satisfactory  in  terms 
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of  storage  space  and  reliability.  It  has  yet  to  be  determined  whether  a  con¬ 
nector  system  of  two  small  diameter  (2.25  mm),  silastic  covered  leads  is  more 
compatible  with  our  particular  chronic  animal  preparation  than  the  single  lar¬ 
ger  (O.D.  2.5  mm)  silastic  insulated  lead  common  with  the  button  connector 
found  on  large  probe  sizes. 

3.  Bench  Calibration:  All  incoming  flow  probes  are  tested  and  bench 
calibrated  before  use.  Probes  in  use  are  tested  for  electrical  continuity 
between  implants  and  are  periodically  calibrated  using  saline  and/or  blood. 

The  purpose  of  the  bench  calibration  is  to  evaluate  the  linearity  of  the 
probe  with  age  (especially  for  low  flows),  and  to  provide  a  flow  calibration 
number  (CAL)  which  can  be  compared  with  in  vivo  dye  curves.  A  listing  of 
some  of  these  results  is  found  in  Table  A-2  which  also  illustrates  some  of 
the  variables  employed  for  bench  calibration. 

For  larger  probes  (18  to  22  mm  I.D.)  cellophane  dialysis  tubing  has  been 
used  as  a  substitute  for  the  blood  vessel  wall.  Soaked  for  approximately 
thirty  minutes  prior  to  testing  in  normal  saline,  the  tubing  is  connected  into 
a  gravity  flow  system  where  the  pressure  at  the  probe  level  is  held  constant 
at  125  mmHg.  Besides  being  convenient,  the  dialysis  tubing  forms  a  stable  seat 
for  the  cuff-type  probes  tested.  Multiple  probe  testing  can  be  accomplished 
with  this  system  providing  care  is  taken  to  minimize  possible  electrical  inter¬ 
ference  between  nearby  probes. 

Typical  bench  calibrations  using  normal  saline  and  blood  of  varying  hema¬ 
tocrits  (derived  from  packed  red  cells  mixed  with  normal  saline)  are  plotted 
for  two  Zepeda  flow  probes  in  Figures  A-6  and  A-7.  Representative  straight 
lines  fitted  to  the  points  (lea3t  squares  method)  indicate  that  in  both  plots 
(except  for  one  case),  the  probes  had  greater  sensitivity  with  saline.  Higher 
hematocrits  tended  to  result  in  lower  sensitivity.  For  other  probes,  the  di¬ 
rection  of  sensitivity  when  using  dialysis  tubing  has  been  consistently  the 
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same.  The  overlap  of  points  indicates  other  factors  that  may  influence  the 
curve  slope  such  as  daily  variations  in  methodology  (Fryer  and  Sandler,  1971; 
Cutherbertson  and  Gilfillon,  1971;  Cappelen,  1968).  No  definitive  pattern 
has  been  seen  in  relation  to  the  duration  of  implantation  and/or  age. 

For  the  smaller  monkey  probes,  a  variety  of  materials  have  been  used  to 
anchor  the  probe  for  calibration  trials.  Dialysis  tubing  has  only  been  avail¬ 
able  for  the  6mm  size.  Blood  vessels  (abdominal  aorta  of  small  pigs)  and 
latex  tubing  have  been  used  for  the  other  sizes.  The  probe  is  inserted  within 
this  latter  material  and  the  leads  are  brought  out  through  a  sealed  Joint. 

For  the  calibration  plot  of  a  Zepeda  (7mm  I.D.)  symmetrical  probe  appear¬ 
ing  in  Fig.  A-8,  the  probe  was  placed  inside  the  tubing.  No  specific  trend  in 
sensitivity  is  evident  in  the  plot.  Similar  probes  calibrated  from  a  position 
within  the  flow  stream  have  yielded  results  also  indicating  saline  produced  a 
greater  relative  output.  Bench  tests  using  saline  and  excised  pig  aortas  re¬ 
sulted  in  increased  sensitivity,  compared  to  the  system  employing  direct  blood 
electrode  contact. 

Flow  probes  having  a  lumen  diameter  less  than  6mm  and  intended  for  peri¬ 
pheral  vessel  use  have  usually  been  calibrated  on  an  intact  vessel  by  bleed- 
out  volume/time  determinations.  When  done  either  acutely  or  chronically,  this 
procedure  provides  an  adequate  calibration  factor  in  lieu  of  dye  curve  deter¬ 
minations  which  become  more  complex  for  peripheral  vessels.  Bench  calibrations 
on  small  prooes  have  been  accomplished  using  dialysis  tubing  or  latex  tubing 
(probe  in  flow),  but  such  techniques  are  used  only  for  linearity  evaluations  and 
CAL  signal  comparisons. 

4.  Dye  Curve  Technique;  A  standard  procedure  for  calibrating  chronically 
implanted  electromagnetic  flow  probes  by  the  dye-dilution  technique  has  been 
developed  in  the  laboratory.  Chronically  implanted  catheters  in  the  right  atrium, 
left  atrium,  and  subclavian  artery  make  possible  several  combinations  of  dye  in¬ 
jection  and  blood  withdrawal.  Normally,  a  bolus  of  indocyanine-green  dye  (approxi- 


mateiy  2.5  mg  for  dogs  and  1.0  mg  for  monkeys)  is  injected  through  the  dye- 
filled  right  atrial  line  without  flushing.  The  blood-dye  mixture  is  with¬ 
drawn  from  the  left  atrium  at  10  ml/min.  Both  right  and  left  atrial  injec¬ 
tions  combined  with  aortic  pickup  have  given  reasonable  results,  although  in 
several  cases  the  latter  situation  appeared  to  give  suspiciously  high  outputs, 
perhaps  due  to  inadequate  mixing.  The  concentrate  m  of  this  mixture  is  mea¬ 
sured  with  a  Waters  XP-300A  Densitometer  equipped  with  an  XC-302  cuvette  and 
recorded  on  a  Honeywell  Visicorder  (Fig.  B-3) .  The  output  of  the  electromag¬ 
netic  flow  probe  is  recorded  in  both  phasic  (to  establish  flow  zero)  and 
mean  modes  simultaneously  with  the  output  of  the  densitometer.  The  height  of 
the  mean  flow  trace  above  the  phasic  zero  is  measured  during  the  passage  of  dye 
through  the  cuvette  and  this  trace  is  analyzed  to  give  cardiac  output. 

The  analysis  of  the  dye  curve  is  by  the  Stewart-Hamilton  method  (semi-log 
replot  and  extrapolation  under  the  assumption  of  an  exponential  decay) .  An  on¬ 
line  computer  analysis  of  the  curves  is  operational  (see  Appendix  B)  and  allows 
immediate  interpretation. 

5.  Comparison  of  Flow  Probe  Calibration  Results:  When  the  bench  calibra¬ 
tions  are  compared  with  dye  curve  results  for  a  specific  probe  on  a  particular 
animal,  a  wide  range  of  differences  from  the  reference  CAL  signals  is  evident. 
The  dye  curves  have  generally  read  approximately  30%  higher  than  the  bench  Cal 
signal  assigned  the  Zepeda  probes.  This  is  approximately  the  result  expected 
from  theory.  However,  as  evident  in  Table  A-2  this  percentage  figure  is  not 
consistent;  in  fact,  some  probes  having  stable  bench  U'unbers  have  had  a  wide 
range  of  dye  curve  calibrations  assigned  to  them. 

The  varied  dye  curve  versus  CAL  results  are  particularly  noticeable  for 
the  Zepeda  probe  SN922.  There  seems  to  be  little  correlation  between  the  de¬ 
rived  CAL  and  hematocrit,  cardiac  output,  or  cardiac  output  normalized  to  body 
weight.  A  curious  and  unsolved  dilemma  concerns  Dog  /'1607  for  which  two  dye 


determinations  were  done.  This  particular  dog  eventually  died  from  an  aortic 
rupture  at  the  flow  prcbe  site,  indicating  that  aortic  wall  thickness  may  have 
been  decreasing  during  implantation.  Similarly,  Dog  L185  died  from  aortic  rup¬ 
ture  just  two  days  after  the  dye  curves,  and  the  dye  curve  CAL  was  low. 

Variations  in  probe  electrode  contact  and  wall  thickness  during  implanta¬ 
tion  may  account  for  the  range  of  differences  illustrated  in  Table  A-2.  In  the 
previous  literature,  there  are  many  conflicting  findings  and  uncertainties  re¬ 
garding  calibration  techniques  and  dependent  variables.  For  this  reason,  coupled 
with  our  findings,  we  rely  on  the  dye  curve  results  to  yield  a  CAL  signal  value 
since  the  dye  dilution  measurement  is  sufficiently  reproducible  and  is  recorded 
under  in  vivo  conditions, 

B.  Vascular  Occluder 


When  an  accurate  electronic  blood  flow  zero  is  not  available,  mechanical 
occlusion  of  a  vessel  to  obtain  a  zero  flow  baseline  is  essential.  Currently 
we  use  the  occluder  design  outlined  by  Khouri  and  Gregg  (1967).  A  dipping  pro¬ 
cess  yields  small  latex  cuffs  of  varying  sizes  dependent  upon  the  stainless 
steel  mold  used.  In  vitro  and  in  vivo  (dog)  testing  with  two  sizes  is  underway. 

A  smaller  size  (3  mm  diam)  la  planned  for  coronary  artery  implants,  and  a 
larger  size  (5  mm  diam)  for  peripheral  vessels  such  as  the  femoral  artery. 

After  removal  from  the  mold,  the  C-shapod  latex  body  is  cured  for  ten  min¬ 
utes  In  a  steam  autoclave  (15  psi).  An  umbilical  tape  backing  and  "Tygon"  tub¬ 
ing  (l.D.  •  .030",  Q.U.  -  .072",  Norton  Plastic  and  Synthetics  Div.,  Akron,  Ohio) 
are  attached  to  the  latex  with  DAB  cement  (R.M.  Hollingshead  Corp.,  Camden,  N.  J,), 
Stainless  ateel  wire  loops  are  glued  with  DAB  to  each  end  of  the  occluder  so  that 
the  ends  can  be  drawn  together  when  implanted  around  the  vessel. 

These  occluders,  when  tested  following  the  animal’s  post-operative  recovery 
period  of  approximately  two  weeks,  have  provided  only  a  few  successful  occlusive 
cycles  before  malfunction,  usually  the  result  of  a  blow-out.  Preliminary  results 


n  15 


J 


Indicate  the  addition  of  one  or  tvo  more  coats  of  latex  to  the  original  three 


layers  considerably  improves  its  durability  without  sacrificing  its  performance. 

C.  Blood  Pressure 

Recordings  from  the  implanted  pressure  gauges  have  been  reliable  and  stable. 
Our  signal  processing  system  consists  of  Honeywell  strain  gauge  modules  (Accu- 
cata  105)  and  associated  DC  amplifier  units  (Accudata  120  DC  Amplifier).  Vari¬ 
able  bridge  excitation  voltage  and  amplifier  gain  are  available  to  the  operator 
and  are  used  during  our  transducer  evaluation  tests. 

Standard  acceptance  evaluation  teats  are  run  on  all  pressure  gauges  that 
are  ised  in  our  laboratory  as  previously  described  (McCutcheon,  et  al.,  1972a; 
McCutcheon,  et  al.,  1973b).  The  transducers  are  examined  for  sensitivity,  fre¬ 
quency  response,  and  temperature  drift  -  schematics  of  the  teat  and  experimental 
apparatus  are  shown  in  Figure  A-9.  Between  the  transducer's  removal  and  subse¬ 
quent  implantation,  the  sensitivity  is  checked.  Further  testing  (e.g.  tempera¬ 
ture  drift)  is  done  periodically  to  insure  the  transducer  is  operating  properly. 

To  Illustrate  the  results  of  the  above  mentioned  test,  a  chronological 
aerieo  of  t?sta  and  recycling  checks  done  on  an  implantable  pressure  gauge  is 
listed  in  Taole  A-3.  The  transducer  chosen  is  a  Konigsberg  P-21  (SN  15)  although 
many  of  our  ocher  gauger  could  also  have  been  similarly  represented.  Including 
a  total  of  12  chronic  implantations  ranging  from  11  to  56  days,  Table  A-3  repre¬ 
sents  the  transducer's  history  from  its  evaluation  test  during  June,  7971  to 
the  most  recent  examination  in  September,  1973. 

The  column  labelled  "sensitivity"  needs  further  explanation,  for  in  our 
laboratory  we  set  the  excitation  voltage  to  match  a  constant  output  voltage  for 
a  given  input  pressure.  To  accomplish  some  standardization  with  previous  litera¬ 
ture,  the  output  voltages  were  normalized  to  the  excitation  voltage.  Percent 
changes  in  sensitivity  reflects  the  difference  from  implantation  to  subsequent 
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recycling  during  which  the  sensitivity  is  recalibrated.  The  change  in  sensiti¬ 
vity  varied  from  -1.9%  to  +0.8%  with  no  consistent  pattern  seen  as  the  transducer 
aged,  and  the  percent  change  in  sensitivity  remained  small. 

Evaluation  tests  on  the  P-21  015  showed  little  change  from  1971  to  1973.  The 
linearity  test  (including  hysteresis)  was  accomplished  by  subjecting  each  gauge  to 
25mmHg  pressure  increments  from  0  to  309  mmHg  followed  by  25  mmHg  decrements  to  zero 

O 

pressure  at  37  C.  Temperature  stability  teats  include  zero-pressure  drift  measur- 
raents  over  a  32°C  to  42°  C  range  in  a  water  bath,  water  (37°C)  to  ambient  air  to 

O 

water  (37  C)  transients,  and  a  long-term  three  to  four  hour  stability  test  in  a 
37°C  water  bath.  Further  test  specifications  and  methods  may  be  found  in  previous 
articles  (McCutcheon,  et  al.,  1972a.  McCutcheon,  et  al.,  1973b). 

Currently  in  use  as  implantable  pressure  transducers  are  Konigsberg  (Pasadena, 
California)  and  Bio-Tec  (Pasadena,  California)  models.  In  dogs  for  aortic  and  left 
ventricular  placement,  5.0mm  (P-19)  and  7.0mm  (P-21)  Konigsberg  gauges  as  well  as 
6.5mm  (BT-250T)  Blc-lec  transducers  have  been  used.  Smaller  4.0mm  (P-13)  and  4.5mm 
(P-12)  Konigsberg  transducers  have  been  implanted  in  the  monkey's  aorta  and  left 
ventricle.  The  above  gauges  have  been  satisfactory  in  terms  of  performance  and 
durability.  Occasional  transducers  have  failed  our  acceptance  evaluation  tests 
and  were  returned  to  the  manufacturer. 

After  the  pressure  gauge  has  been  tested,  e  silicone  rubber  (General  Electric 
RTV-112  Adhesive)  "washer"  is  molded  behind  the  sensing  element.  A  typical  pre¬ 
paration  is  shown  in  Figure  A- 10.  The  gauge  head  is  then  briefly  treated  with 
TDMAC  and  heparin  anticoagulants  (see  Cannulae  section  for  details)  before  gas 
sterilization. 

A  significant  improvement  in  aortic  wall  integrity  at  the  pressure  gauge 
site  has  occurred  since  placing  a  silicone  rubber  "washer"  between  the  sensing 
element  and  vessel  wall.  Since  the  initation  of  this  method,  16  dogs  have  been 
implanted  with  an  aortic  gauge  and  no  deaths  have  result  d  from  wall  rupture  at 


this  site.  Previously,  damage  at  the  aortic  gauge  site  accounted  for  approxi¬ 
mately  25?  of  the  deaths.  Apparently,  the  motions  resulting  from  the  heart's 
contractions  (and  probably  the  experimental  vibration)  caused  the  sharp  edges 
of  the  pressure  gauge  sensing  head  to  break  down  the  surrounding  tissue.  Pres¬ 
sure  transducers  in  the  left  ventricle  are  also  fitted  with  a  silicone  washer. 

At  the  aortic  implantation  site,  silastic  sponge  material  (Dow  Coming  #812)  is 
also  fitted  around  the  vessel  for  additional  support. 

Clot  formation  over  the  treated  portion  of  the  transducer  within  the  vas¬ 
cular  system  nas  generally  not  been  a  problem.  Proper  placement  within  a  vessel 
or  chamber  is  very  important.  Autopsy  results  have  shown  that  when  the  gauge  ii 
inserted  too  deeply  into  the  vessel  (aorta),  there  have  been  occasional,  yet  well 
developed  clot  formations  anchored  to  the  gauge's  wire  connector.  The  depth  of 
penetration  of  the  pressure  transducer  in  the  left  ventricle  appears  also  to  be 
a  critical  feature,  especially  in  the  monkey.  Clot  formation  in  this  location 
has  been  minimal,  although  when  the  gauge  is  pulled  closely  to  the  wall  at  the 
apex,  a  combination  of  tissue  growth  and/or  thrombus  formation  has  occurred.  With 
this  latter  placement,  or  when  placed  too  deeply  into  the  ventricle,  the  resulting 
ventricular  pressure  trace  has  exhibited  an  abnormal  waveform,  probably  as  a  re¬ 
sult  of  pressure  transients  arising  from  lateral  compression  of  the  transducer. 

None  of  our  dogs  have  shown  such  phenomena.  The  relatively  small  size  of  the  mon¬ 
key's  left  ventricular  chamber  coupled  with  possible  changes  in  contractile  pat¬ 
terns  associated  with  the.  Increased  cardiac  output  of  stress  may  also  account  for 
the  occasional  presence  of  these  sharp  artifacts. 

Both  the  Konlgsberg  and  Bio-Tec  transducer  connectors  have  been  generally 
reliable  and  durable.  When  not  opened  for  a  period  of  two  weeks  during  implanta¬ 
tion,  both  connectors  have  exhibited  some  fl vid  leakages  when  sealed  as  recommended. 
A  commercial  cleaner  (Freon  TF  Degreaser,  Miller  -  Stephenson)  has  worked  well 
for  cleaning  these  female  pin  connectors.  Being  relatively  small  (approximately 


0.25"  diam. ,  1.5''  long)  the  connectors  are  well  suited  for  the  fabric  pouch  used 
in  the  chronic  animals. 

E;:travascular  transducers  have  been  used  for  calibration  of  Implanted  gauges 
or  pressure  measurements  during  an  experiment.  These  transducers  undergo  much  the 
same  evaluation  tests  as  do  the  implantable  gauges.  In  vivo  calibration  of  the 
implanted  aortic  gauge  is  accomplished  by  using  a  gauge  (Ailtech  MS-10B,  City  of 
Industry,  Calif.)  connected  to  the  aortic  cannula.  The  ventricular  gauge  is  cali¬ 
brated  by  measurement  of  left  atrial  pressure  through  the  cannula  to  that  cham¬ 
ber.  A  manometer-tipped  catheter  (PC-350,  Millar  Instrument  Co.,  Houston,  Texas) 
has  been  used  to  measure  pressures  when  vessel  cutdowns  are  necessary  or  appropri¬ 
ate.  In  contrast  to  the  system  employing  a  fluid  filled  line,  the  catheter-tip 
gauge  exhibits  minimal  vibration  artifact. 

D.  Pacemaker  Probe 

Since  some  of  the  more  recent  investigations  require  control  of  heart  rate, 
a  pacing  probe  has  also  been  implanted  routinely.  Earlier  attempts  to  place  a 
pacing  catheter  acutely  (femoral  and  external  jugular  vein  entry)  prior  to  the  ex¬ 
periment  wei  .  not  always  successful,  especially  when  the  animal  was  subjected  to 
vibration.  Capture  of  the  cardiac  pacemaker  was  intermittent,  probably  due  to 
failure  of  the  pacemaker  electrode  to  seat  itself  firmly  in  the  right  ventricle 
during  vibration.  To  pace  the  animal  we  have  been  using  a  Grass  Stimulator  (Model 
S4G)  connected  to  the  pacemaker  lead  via  an  isolation  unit  (Grass  Model  SIU-4B) . 
The  bipolar  pacing  lead  has  also  allowed  relatively  artifact  free  monitoring  of  an 
ECG  signal  during  non-pacing  experiments.  No  pacemakers  have  been  implanted  in 
monkeys . 

During  implantation  the  pacemaker  wire  is  sutured  to  the  surface  of  the  right 
atrial  wall.  A  commercial  semi-floating  pacing  probe.  (Elecath  //561,  4  Fr.,  100cm) 
has  been  implanted  in  the  majority  of  dogs,  and  because  this  pacemaker  is  intended 
to  be  disposable,  and  is  in  fact  quite  fragile,  a  new  pacemaker  must  be  used  for 


each  dog.  Recently,  we  have  also  used  silicone  rubber  insulated  wires  (#29-51  x 
46  Stranded  Copper  Conductor,  Caltron  Industries,  Berkeley,  Cal.)  terminated  on 
the  atrial  end  by  a  soldered  ring  and  on  the  pouch  end  by  a  small  male  connector 
fitting.  This  latter  approach  again  employs  two  conductors  (bipolar  pacing)  and 
we  further  coat  the  wires  with  medical  grade  Silastic  (#861,  Dow  Corning). 

Further  studies  need  to  be  accomplished  to  determine  if  lead  corrosion  over  a 
nominal  implantation  time  of  30  days  adversely  affects  the  pacing  pulse.  To 
date,  qualitative  results  have  been  satisfactory. 

E.  Cannulae 

Placement  of  cannula  directly  into  the  right  and  left  atria  of  both  monkeys 
and  dogs  has  become  a  routine  feature  of  the  animal  surgery.  To  calibrate  an 
implanted  pressure  gauge  located  Just  below  the  arch  in  the  descending  aorta,  an 
additional  cannula  is  inserted  into  the  left  subclavian  artery  which  is  then  tied 
off.  With  these  cannulae,  mixed  venous  and  arterial  blood  samples  can  be  taken. 
When  an  extravascular  pressure  transducer  is  connected  to  the  left  atrial  line, 
a  pressure  equivalent  to  left  ventricular  end  diastolic  pressure  may  be  obtained. 

The  reliability  of  the  cannulae  for  blood  sampling  and.  pressure  calibration 
has  varied  depending  upon  the  animal,  the  care  taken  in  placement,  and  cannula 
design  and  construction.  The  tubing  used  in  all  our  cannulae  is  "Tygon"  plastic 
tubing  (Norton  Plastic  and  Synthetics  Div. ,  Akron,  Ohio)  having  an  I.D.  of  .038" 
and  O.D.  equal  to  .072".  Joints  between  the  tubing  and  fitted  materials  (e.g. 
polyethylene  y-yoke)  are  made  with  DAB  cement  (R.M.  Hollingshead  Corp . ,  Camden, 

N.  J.),  Resistance  to  kinking  has  been  excellent  and  iu  has  held  up  quite  well 
for  implantation  times  of  four  months.  The  anti-coagulant  treatment  (see  below) 
tends  to  reduce  the  tubing's  plasticity,  but  not  enough  to  cause  kinking  or  split¬ 
ting  after  implantation.  Tissue  growth  over  the  cannula  tips  in  the  atria  has 
caused  occasional  problems  by  acting  as  a  one-way  valve.  Different  shapes  for  the 
inlet  portion  of  the  cannula  are  still  being  investigated.  The  portions  of  the 
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cannulae  which  would  be  placed  in  the  atrial  chambers  are  illustrates  in  Fig¬ 
ure  A-ll.  Both  designs  are  being  used  presently,  the  straight  tip  in  the  left 
atrium  and  y-yoke  in  the  right  atrium.  Getting  a  consistently  reliable  right 
atrial  sample  has  posed  a  still  unsolved  dilemma  as  autopsies  have  regularly 
demonstrated  rampant  tissue  growth  over  any  cannula  tip  in  the  right  atrium 
while  the  left  atrial  tip  has  remained  free.  The  y-yoke  tip  (s)  present  a  large 
surface  area  for  sampling  -  in  fact,  their  use  has  been  prohibited  with  some  of 
the  smaller  Rhesus  monkeys  -  but  the  same  type  of  one-way  valves  have  formed.  At 
present,  our  solution  to  maintaining  right  atrial  sampling  has  been  to  check  cath¬ 
eter  patency  as  soon  as  possible  after  surgery,  and  to  flush  daily.  The  cannulae 
inserted  through  the  subclavian  artery  (Fig.  A-ll)  into  the  aorta  do  not  seem  to 
have  the  problem  of  tissue  dockage  and  have  been  very  reliable  for  sampling  and 
pressure  measurements. 

Another  method  for  obtaining  mixed  venous  blood  using  the  cannula  in  Fig.  A-ll 
has  been  in  the  testing  stage  to  date.  Entry  through  a  branch  of  the  vena  cava, 
such  as  the  azygous  vein,  into  the  atrium,  and  then  possibly  across  the  tricuspid 
valve  into  the  right  ventricle  has  shown  some  preliminary  success.  Apparently 
less  tissue  can  accumulate  on  a  "floating"  type  of  cannula  as  opposed  to  the  tubing 
rigidly  attached  to  the  atrial  wall. 

Anticoagulant 

A  coating  of  "TDMAC"  (Tridodecylmethylammonium  Chloride  -  Polysciences,  Inc.) 
and  heparin  prior  to  implantation  has  been  successful  in  preventing  clot  formation 
on  the  implanted  tubing.  A  two-step  process  is  currently  used  (Grode,  1972),  and 
the  tubing  is  allowed  to  dry  thoroughly  between  and  after  treatments.  Because  the 
TDMAC  solution  tends  to  remove  plastomers  from  the  "Tygon"  tubing  (Leinlnger,  1972), 
the  cannula  tip  is  briefly  dipped  into  and  flushed  with  the  TDMAC  solution.  After 
drying,  the  tubing  is  soaked  in  and  flushed  with  the  heparin  solution  for  approxi- 
rately  10-15  seconds,  then  allowed  to  dry  before  gas  sterilization. 


F.  Accelerometers 


To  enhance  the  study  of  vibration-induced  acceleration,  individual  organ 
accelerations  may  be  examined  with  an  implanted  commercial  accelerometer  (Koniga- 
berg  single  axis  accelerometer,  Models  A1  and  A2) .  Data  have  been  obtained  from 
the  ascending  aorta  (.rigidly  attached  to  the  flow  probe)  and  from  the  descending 
thoracic  aorta  below  the  arch.  To  facilitate  surgical  placement  and  allow  axis 
identification,  a  silastic  base  was  molded  to  the  accelerometer  at  a  cost  of 
adding  some  mass  to  the  accelerometer .  Further  development  of  surgical  technique 
and/or  smaller  silastic  anchoring  hinges  is  progressing  to  reduce  that  particular 
drawback.  Ocher  organ  acceleration  sites  are  planned,  but  these  should  present  nc 
implantation  problems. 

The  implanted  accelerometer  received  in  our  laboratory  is  subjected  to  eval¬ 
uation  and  calibration  tests  (some  of  which  are  similar  to  those  involving  im¬ 
plantable  pressure  transducers  described  previously).  Static  tests  Include  tem¬ 
perature  drift,  ±  1G  sensitivity  ratings,  and  cross-axis  responses.  Dynamic  sen¬ 
sitivity,  frequency  response,  and  cross-axis  response  are  tested  using  Che  vibra¬ 
tion  system  facilities  and  precalibrated  accelerometers  in  the  Wenner-Gren  Labora¬ 
tory.  The  Konigsberg  accelerometers  have  generally  been  satisfactory  (although 
one  accelerometer  developed  an  erratic  sensitivity  due  to  an  internal  malfunction 
shortly  after  receipt).  An  in-vivo  i  1G  calibration  for  an  experiment  is  accom¬ 
plished  by  rotating  the  animal  in  the  desired  axis. 

G.  Temperature 

To  record  an  accurate  internal  temperature,  a  precalibrated  thermistor  bead 
(Yellow  Springs,  Inc.)  has  been  implanted  in  the  thoracic  cavity  and  sutured  to 
the  muscle  layers.  Earlier  troubles  with  fluid  leakage  into  the  connector  wire 
and  wire  breakage  have  been  nearly  eliminated  by  encasing  the  lead  in  flexible 
silastic  tubing,  sealing  the  bead  end  with  silastic  adhesive,  and  closing  the 
connector  end  with  suture  tie  and  ahrink-tight  tubing.  In  cases  where  an  intra- 


thoracic  probe  was  not  implanted,  rectal  temperature  has  been  measured  via  an 
acutely  placed  probe  (Yellow  Springs,  Inc.)*  This  recording  has  not  been  as 
reliable  as  information  from  the  implanted  sensor  since  we  have  had  some  pro¬ 
blem  maintaining  a  standard  anatomical  (and  temperature)  reference  point  dur¬ 
ing  the  z~axis  vibration.  Either  a  Honeywell  Temperature  Module  or  Yellow 
Springs  Scanning  Telethermometer  serves  as  the  first  interface  between  the 
probe  and  paper  chart  recorder  and/or  computer. 

H.  Implant  Support  Materials 

Included  in  this  section  are  the  qualitative  evaluations  of  materials  asso¬ 
ciated  with  the  devices  Implanted  in  the  animals.  The  investigation  revolves 
primarily  around  the  use  of  nylon  and  dacron  velour  fabrics.  Determination  of 
the  success  of  the  implant  material  has  depended  upon  the  results  of  tissue  exam¬ 
inations  and  gross  observations  made  during  the  chronic  animal's  implantation 
period . 

In  the  majority  of  implanted  animals,  the  skin  tissue  -  velour  fabric  pouch 
margin  has  formed  a  tight  seal  upon  both  gross  and  microscopic  examinations  as 
summarized  in  Table  A-4.  Sections  generally  uhow  skin  with  an  unremarkable  epi¬ 
dermis  and  upper  dermis.  Further  microscopic  findings  indicate  the  pouch  fibers 
are  embedded  in  fibrous  tissue  arising  from  the  subcutaneous  layers. 

Excepting  instances  where  many  pouch  openings  can  be  expected  (i.e.,  most  of 
our  monkeys  have  required  almost  daily  attention  to  maintain  right  atrial  cannula 
patency) ,  the  double  layer  nylon  velour  pouch  has  become  the  standard  pouch  im¬ 
planted.  The  dacron  pouch  serves  as  an  alternative  as  it  is  a  stronger  material 
over  long  implant  times  than  nylon.  Different  tissue  ingrowth  patterns  have  not 
been  noted  between  these  two  velour  fabrics,  but  our  experience  with  dacron  has 
been  limited. 

Because  of  che  nylon  velour's  capacity  r.o  form  a  firm  fabric-tissue  inter¬ 
face,  this  material' 8  use  has  been  extended  to  other  areas  of  the  surgical  pre- 


paracion.  Autopsy  examination  of  the  velour  covered  leads  has  revealed  that 
an  effective  tissue  seal  has  formed  along  what  would  have  been  the  lead  track. 

The  extent  of  the  foreign  body  (tubing,  DAB  cement,  velour)  reaction  has  not 
been  sufficiently  defined  at  present.  Without  the  velour,  various  degrees  of 
fistulae  have  formed  along  the  lead  sites  regardless  of  the  composition  of  the 
implant  devices'  leads  we  currently  employ:  silastic  coated  wires,  polyvinyl¬ 
chloride  insulation,  "Tygon"  tubing,  and  various  other  plastics. 

The  addition  of  a  velour  layer  to  the  leads  does,  however,  create  more  lead 
volume  and  more  tissue  must  be  disturbed  when  routing  leads  from  the  chest  to 
the  pouch.  If  bacteria  invade  the  implanted  velour,  the  treatment  of  the  infec¬ 
tion  is  a  much  more  severe  problem  than  if  the  material  were  not  there,  for  the 
fabric  site  containing  the  bacteria  19  nearly  inaccessible.  Whether  the  addition 
of  velour  to  the  leads  is  beneficial  to  the  overall  chronic  animal  preparation 
has  not  yet  been  determined. 

III.  NON-TNVASIVE  MEASUREMENT'S 

A.  ECG  &  Temperature 

Depending  upon  the  particular  experimental  protocol,  skin  electrodes  (Beckman) 
are  also  available  for  external  measurement  of  ECG.  Skin  temperature  probes  (Yellow 
Springs,  Inc.)  are  applied  for  measurements  of  surface  temperature. 

B.  Oxygen  Consumption 

Oxygen  consumption  has  been  measured  on  both  unanes Lhetized  and  anesthetized 
animals.  A  flow-through  system  (Fig.  A-12)  is  employed  for  awake  dog9  and  monkeys; 
that  is  room  air  diawn  past  the  animal's  head  will  show  a  decrease  in  oxygen  content 
dependent  on  the  amount  of  consumed  by  the  animal.  In  an  anesthetized  prepara¬ 
tion,  oxygen  consumption  is  determined  from  the  known  respiratory  volume  and  rate 
and  by  sensing  the  partial  pressure  of  oxygen  in  the  output  port  of  the  respirator. 

Further  helmet-mask  development  for  the  dog  (Fig,  A-13)  is  occurring  with  em¬ 
phasis,  as  was  done  for  the  monkeys,  on  obtaining  a  uniform  room  Air  flow  rate 


past  the  exhalation  site. 


Not  all  dogs  have  accepted  the  plexiglass  helmet  upon  initial  trials  and 
it  has  been  necessary  to  acquaint  the  dogs  with  the  helmet.  The  helmet  in  its 
preliminary  design  is  secured  to  the  vibration  restraint  chair  with  tape  and 
padding.  Helmet  design  requires  allowance  for  space  tor  the  dog  to  pant  as  they 
do  during  a  vibration  protocol.  Our  brief  experience  with  a  latex  mask,  which 
did  not  have  enough  dead  space  to  allow  panting,  showed  that  a  deg's  internal 
body  temperature  rose  noticeably  in  a  short  time  period. 

The  monkey  helmet  (Fig.  A-14)  is  formed  from  lightweight,  transparent  plastic 
and  is  secured  to  the  vibration  restraint  chair  by  a  large  hose-clamp.  With  the 
present  design  for  monkeys,  we  use  two  different  size  helmats  depending  on  the 
size  of  the  monkey.  The  bottom  plate  of  the  helmet  has  holes  allowing  a  controlled 
room  air  inflow  and  a  rubber  seal  which  fits  loosely  around  the  monkey's  neck. 

Room  air  flow  rates  are  controlled  by  an  electric  blower  and  monitored  contin¬ 
uously  by  a  Fleisch  pneumotachograph  (Instrumentation  Associates,  Inc.)  connected 
to  a  differential  pressure  gauge  (Statham) .  This  combination,  after  calibration, 
results  in  a  hard-copy  recorder  printout.  For  the  monkeys,  an  air  flow  rate  of 
10  L/min  is  generally  used.  For  dogs,  this  is  increased  to  15  or  20  L/min.  Flow 
rates  below  these  values  tend  to  compromise  the  quality  of  the  mixing  of  exhaled 
and  room  air  plus  altering  the  normal  partial  pressures  of  inhaled  gases.  The 
sensitivity  of  the  measurement  decreases  with  increasing  flow  rate. 

The  mixed  ratio  of  exhaled  air  and  room  air  is  then  sampled  continuously  with 
a  vacuum  pump  at  approximately  400  ml/min.  To  eliminate  varying  temperature- 
humidity  effects  in  the  open  system,  before  entering  the  oxygen  sensor  cr'l  (Applied 
Electrochemistry  Inc.),  this  sampled  mixture  is  pulled  through  a  drying  tube.  A 
readout  voltage  proportional  to  percent  oxygen  content  is  then  available  for  the 
recorder  and  computer  (where  0^  consumption  over  a  selected  period  of  time  is  cal¬ 


culated)  . 


Calibration  of  the  air-flow  system  is  accompli,  hed  with  the  use  of  a  bubble- 
flowmeter  (precalibrated)  and  vacuum  'ine.  The  0  sensor  has  been  periodically 
checked  with  tank  oxygen  concentrations :  a  high  end  of  98%  0^  and  low  end  of  11% 
are  used.  Daily  standardization  is  accomplished  by  reading  "dry"  room  air. 

The  response  time  of  the  AEI  thermaL  cell  sensor  is  in  milliseconds,  but  the 
overall  response  times  of  the  flow-through  system  are  much  larger.  Depending  upon 
the  dead  space  volume  in  front  of  the  sampling  line,  the  delay  In  generally  2  to 
5  seconds.  The  drying  tube  adds  another  6  seconds  to  the  response  time  and  there¬ 
fore  the  final  breathing  pattern  seen  is  quite  damped.  To  determine  a  respiration 
rate  more  accurately  using  this  method,  the  air  in  the  system  can  be  heated  to  the 
internal  body  temperr.ture  and  drawn  through  the  sensor  whila  saturated  with  water 
vapor.  At  present,  this  latter  method  Is  being  developed. 


IV .  Summary  of  Chemistry-Hematology  Measurements 


Depending  upon  the  animal's  condition  and  the  experimental  protocol, 
the  following  tests  as  summarized  are  conducted  on  a  routine  basis : 

1.  pH,  pC>2»  pC(>2  (Radiometer  determination  -  Radiometer  A/S, 
Copenhagen,  Denmark). 

2.  and  C02  content  (Van  Slyke  technique  or  nomogram  calculations 
from  radiometer  and  hemoglobin  measurements). 

3.  Free  Fatty  Acids  (modified  technique  from  No^ak,  1965,  in 
Laffertv,  et  al.,  1973). 

A  Cortisol  (competitive  protein  binding  assay  of  Murphy,  1967 
and  Beamer,  et  al. ,  1972). 

5.  Glucose  (glucose  oxidize  -  peroxidase  method-GOD-Perid , 
Boehrlnger  Mannheim  Corp.  ,  N.  Y.,  N.  Y.). 

6.  Hemaglobln  (cyanmethemoglobln  method  -  Globintest,  Pfizer 
Diagnostics  Dlv. ,  N.  Y. ,  N.  Y.)  and  Hematocrit  (capillary 
tube) . 

7.  White  Blood  Count  (hemacytometer)  and  Differential  WBC 
.ount  (slide  using  Quik  Stain,  Chemstat,  Inc.,  Inglewood, 

Calif.) . 

8.  Other  hematology  testa  3uch  as  sedimentation  rate,  red  cell 
fragility,  and  reticulocye  count  ar;  performed  when  indicated. 


METRIC  1 


Fig.  A-l  Isolated  nylon  velour  pouch  with  inner  surface  exposed  to  show  construction 
details  of  double  layer  and  overall  shape. 


CaSbration  of  18mm  Zepeda  Flow  Probe  SN923 


Output  Voltage  of  Flowmeter'  (Volts) 

Least  squares  fitted  calibration  curves  for  18  ram.  Zepeda  probe  SN  923 


Output  Voltage  of  Flowmeter  (Volts) 


Fig.  A-10  Silicone  rubber  washers  are  used  as  an  interface  between  the  vessel 

wall  and  the  gauge  head  of  Bi.o-Tec  and  Konigsberg  pressure  transducers 
to  be  implanted. 
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Fig.  A  - 1 2 


Clock  diagram  of  open  system  to  measure  oxygen 
consumption  of  a  conscious,  vibrating  animal. 
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BE1JCH  CALIBRATIONS  USING  VARIOUS  FLOW  MEDIA  AND  DYE 
CURVE  CALIBRATION  RESULTS  ARE  COMPARED 
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TABLE  A -3 


CHRONOLOGICAL  HISTORY  OF  A  KONIGSBERG  P21  (SN15)  IMPLANTABLE 
PRESSURE  TRANSDUCER  CURRENTLY  USED  IN  THE  LABORATORY.  VERY 
SMALL  CHANGES  IN  SENSITIVITY  ARE  EVIDENT  AND  THE  LATER  TESTS 
DONE  ON  THE  GAUGE  SUGGEST  THERE  IS  MORE  TEMPERATURE  STABILITY 
WITH  LONGER  GAUGE  LIFE. 
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F.  S.  LINEARITY 

5k 0.47.  (6-10-71) 

32<,C-42°C  ■  -.25  nsnHg  (6-16-/1),  +2  B*nHg  (6-10-71)per  »c. 

TEMP.  STABILITY 

W.A.W.-+13.5  nmHg  (6-15-71) 

(V  ■  5V.) 

4  hre.-  ltanHg  (6-10-71),  0  nmHg  (11-10-71)  at  37*C 

SENSITIVITY 

101.1  yv/v/cmHg  at  1GV  Excitation  (6-10-71) 

FREQUENCY  RESPONSE- 

> 100  Hz  (6-10-71)  .  ' 
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•  <±  0.47.  (9-27-73) 

32  =  L-42‘t;  -  •+■  1.33  ramHg/" 

-  /  a  ne 

TEMP,  STABILITY 

-  W.A.W.-<1  tunHg  (9-28-73) 

(Vex-8.35V) 

SENSITIVITY 

3  hrs .  -  -0.5  cnHg  at  37” 

C  (9-28-73) 

-  100 . 6yv/v/cmHg  (9-27-73) 

at  8.35V  Excitation 

FREQUENCY  RESPONSE 

-  ">100  Hz  (9-27-73) 
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RAYTHEON  704  DATA  ACQUISITION  SYSTEM 


A .  704  Description 

The  Raytheon  704  Computer  System  is  a  general  purpose  digital  computer 
system  primarily  intended  as  the  central  element  for  data  acquisition  and 
control  applications.  Some  important  hardware  characteristics  are  16  -  bit 
word  length,  1.0  microsecond  cycle  time,  automatic  priority  interrupt,  direct 
input/output  to  control  processor,  direct  memory  access  and  a  full  line  of 
peripheral  options.  Only  assembly  language  is  supported.  The  configuration 
of  the  U.K.  system  is  shown  in  Figure  B-l. 

The  direct  memory  access  feature  allows  peripherals  to  communicate  directly 
with  core  memory  independently  and  simultaneously  with  any  program  which  may  be 
executing.  Slower  speed  peripherals  communicate  with  the  CPU  via  the  DIO  bu9  . 

The  16  priority  interrupts  may  be  connected  to  Raytheon  peripherals  or  user 
devices  such  as  the  'start'  and  'stop'  buttons  or  a  trigger  from  the  ECG  of  the 
animal  under  study.  Several  peripherals  such  as  a  time  code  generator  controller, 
experiment  status  switches,  and  a  digital  link  to  a  remote  computer  will  be 
connected  co  the  DIO  bus  and  the  priority  interrupts. 

B .  Mean  Value  Program 

The  programs  of  this  system  are  designed  to  give  a  real-time  print-out  of 
the  average  values  of  up  to  five  analog  input  channels  and  to  store  these  average 
on  digital  magnetic  tape.  The  voltages  on  these  channels  represent  slowly  chang¬ 
ing  parameters  such  as  body  temperature,  mean  blood  flow,  etc.  A  sampling  rate 
of  256  readings  per  10  seconds  (25.6  Hz)  was  selected  with  a  minimum  print-out 
interval  of  10  sec.  Sampling  at  a  rate  of  1  Hz  would  have  been  sufficient  but  we 
wanted  to  examine  the  9peed  of  the  computer.  The  print-out  interval  may  be  set 


at  any  integral  multiple  of  10  sec.  but  the  data  are  stored  on  magnetic  tape  every 


10  sec.  The  printed  values  always  represent  the  average  taken  over  the  period 
since  the  previous  print-out;  a  typical  print-out  appears  in  Fig.  B-2. 

The  programs  in  this  system  were  designed  to  be  very  interactive  and  almost 
completely  dependent  upon  interrupts  for  execution.  All  programs  are  in  core  at 
all  times  and  a  lower  priority  program  is  interruptable  without  damage,  by  one 
of  higher  priority.  The  704  has  sixteen  interrupt  levels  which  may  be  connected 
to  external  devices  such  as  the  'start'  and  'stop'  buttons  or  peripherals  such  as 
a  timer  or  the  analog  to  digital  converter.  Whenan  interrupt  occurs,  the  com¬ 
puter  begins  execution  at  the  specific  location  assigned  to  that  level,  provided 
no  higher  priority  level  is  in  service. 

Once  the  programs  are  loaded,  IN1TALZ  executes  automatically  and  performs 
various  initilization  chores,  enables  the  'start'  button,  starts  the  clock,  and 
exit;;  to  the  wait  loop  (WAIT).  In  this  wait  loop  the  sense  switches  are  checked 
and  the  type  buffer  is  examined  to  see  if  any  printout  is  requested.  When  no 
other  program  is  in  execution,  the  computer  runs  in  this  loop  continuously. 

Calibration  is  performed  by  reading  the  voltage  corresponding  to  at  least 
two  (maximum  four)  points  in  the  range  of  each  input  channel.  Since  the  inputs 
must  be  linear  (the  equation  y  *  Mx  +  B  is  used)  the  values  M  and  B  are  calculated 
using  a  least  squares  fit  over  the  points  which  were  read  and  then  stored  in  the 
calibration  tables.  The  print-out  will  be  in  the  appropriate  physiological 
units . 


When  the  'start'  button  is  pressed, a  program  (START)  starts  the  A/D  converter 
and  then  returns  to  WAIT.  When  the  A/D  converter  is  through  with  the  requested 
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brings  CALC  into  execution. 


15-2 


CALC  Is  che  program  which  performs  the  arithmetic  on  the  data.  The  volt¬ 
ages  just  read  are  added  to  their  summing  buffers.  If  it  is  time  to  print-out 
or  to  write  on  magnetic  tape,  these  sums  are  divided  to  get  the  average  and 
multiplied  by  the  calibration  factor  to  give  output  units  of  beats  per  minute, 
mmHg,  etc.  Then  they  are  stored  in  the  type’  or 'magnetic  tape’ buffers .  CALC 
also  starts  a  programmable  timer  whose  period  is  equal  to  the  desired  sampling 
period  minus  the  program  execution  and  A/D  conversion  times.  When  this  timer 
times  out  an  interrupt  causes  PIT  to  start  the  A/D  converter  again. 

This  cycle  is  repeated  (approximately  every  40  milliseconds  when  the  sample 
rate  is  25.6  Hz)  until  the  'atop'  button  is  pressed.  The  computer  then  loops 
in  WAIT  until  the  'start'  button  is  pressed  again. 

C.  Dye  Curve-  Program 

Measurement  of  cardiac  output  by  the  dye  dilution  technique  is  one  way  to 
calibrate  chronically  implanted  electromagnetic  flow  probes  (Rothe,  1973).  The 
classical  Stevart-Hamil ton  analysis  requires  a  plot  of  the  downslope  of  the  dye 
curve  on  semilog  paper  with  subjective  assessment  of  the  exponential  decay  con¬ 
stant.  According  to  this  theory,  the  shape  of  the  downslope  should  be  an  expon¬ 
ential  decay  until  recirculation  occurs.  Therefore,  fitting  an  exponential  to  the 
downslope  before  recirculation  and  extrapolating  to  the  baseline  should  give  an 
accurate  indication  of  dye  disappearance  with  time.  The  area  under  the  first  part 
of  the  curve  is  added  to  the  area  under  the  extrapolated  exponential  with  the  sum 
being  inversely  proportional  to  cardiac  output. 

Application  of  this  same  method  with  a  digital  computer  is  a  relatively 
straightforward  matter.  The  signals  from  the  densitometer  (Waters  Model  XP300 
with  Cuvette  Model  XC-302)  after  passing  through  analog  signal  conditioners  illus¬ 
trated  in  Fig.  B-3  are  sampled  by  the  A/D  converter  at  10  Hz  (Stenson,  Crouse, 


Harrison,  1972).  The  computer  must  have  at  least  one  calibration  where  a  quantity 
of  blood  with  a  known  dye  concentration  is  passed  through  the  densitometer  cuvette 
while  the  voltage  is  recorded.  For  the  monkeys  and  dogs  under  study,  concentrations 
of  2.5,  5.0,  and  7.5  mg  dye/liter  blood  are  normally  used.  The  calibration  signals 
are  monitored  on  an  oscilloscope  by  an  operator  who  presses  'start'  and  'stop' 
(interrupt  level)  buttons  to  allow  the  computer  to  read  only  the  flat  portion 
of  the  curve.  Each  calibration  signal  is  read  for  a  period  of  5  to  15  seconds. 

Due  to  some  drift  arising  from  the  densitometer  system,  operator  control  is 
necessary.  Following  the  calibration  levels,  a  second  pure  blood  zero  is  read 
which  allows  the  operator  to  determine  the  magnitude  of  the  shift,  if  any,  in 
the  blood  zero  baseline  and  gives  a  more  meaningful  calibration  factor  when 
averages  are  computed. 

The  dye  curve  is  read,  as  above,  under  operator  control  at  a  rate  of 
10  Hz  and  recorded  point  by  point  in  an  array  (300  point  maximum).  When  the 
1 8  top '  button  is  pressed,  the  computer  locates  the  peak  point  of  the  curve, 
and  then  by  a  process  of  iteration  the  blood  zero  baseline  prior  to  the  injection 
of  dye  is  calculated.  This  value  is  subtracted  from  each  point  in  the  curve 
to  eliminate  the  offset. 

The  area  under  the  dye  curve  is  then  calculated.  Starting  after  the  peak 
at  a  magnitude  of  0.9  peak  voltage,  the  curve  is  examined  in  3  second  overlapping 
segments  (30  points  each).  Each  segment  is  fitted  to  a  true  exponential  by  the 
least  squares  method  and  variance  and  exponential  decay  constant  (t)  are  calculated 
and  recorded.  After  each  segment  has  been  examined,  the  constant  t  which  had  the 
smallest  variance  is  selected  and  an  exponential  decay  curve  is  extrapolated  from 
that  point  down  to  the  baseline  (zero).  The  area  under  the  first  part  of  the  curve 


is  calculated  using  Simpson's  Rule  and  added  to  the  exponential  area.  The  cardiac 


output  is  then  calculated  with  the  following  equations: 


CAL 


mg  dye/llter  blood 
volts 


amount  of  dye  injected 
area  under  curve  x  CAL 


C  .0 .  =  — - - “ - 77 -  =  liter/minute 

volts  x  minutes  x  mg/ liter 

vol  ts 

The  cardiac  output  is  printed  out  along  with  several  intermediate  results  - 
see  sample  output  in  Fig.  B-4.  The  operator  may  now  choose  to  read  another  dye 
curve  using  the  same  calibration,  read  another  calibration,  or  exit  to  the  system 


monitor . 


Fig.  B-L  Block  diagram  of  Raytheon  704  system  currently  in  use  showing 
peripheral  communication  with  the  CPV. 
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Fig.  B-2  Printout  of  MEAN  VALUE  program  showing  five  mean  variables:  heart  rate  (BPM) , 

left  ventricular  pressure  (mmHg) ,  aortic  flow  (L/min) ,  oxygen  consumption  (ml/min) , 
internal  thoracic  temperature  (°C).  The  data  represents  control  and  vibration 
variables  for  a  Rhesus  monkey. 
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Fig.  B-3  Block  schematic  of  Instrumentation 
curves  on  both  dogs  and  monkeys. 


CV#  ID  #  DATE 

XX*  *XXXX  *  *XX /XX /XX 

01  1607  1 1  02  72 

SPAN*  185  FR EQ ■  >0  CAL*  1.8851 

BASE*  .2163  PEAK*  4.1085  AT  i»  96 


TAU 

VARIANCE 

I 

TO  I 

-.0268 

•0034 

106 

136 

-.0282 

•  0024 

1  1  1 

141 

-.0284 

•0022 

1  16 

146 

-.0281 

.0026 

121 

151 

-.0273 

.0026 

126 

156 

-.0271 

•0025 

131 

161 

-.0267 

.0026 

136 

i66 

-.0252 

•  0041 

141 

171 

-.0219 
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Fig.  B-4  Representative  printout  of  on-line  dye  curve  program 
illustrating  output  of  two  dye  cartes  done  on  Dog 
3607.  Intermediate  values  other  th.m  ctrdiac  output 
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PREFACE 


1.1  Ihe  purpose  of  this  document  is  to  present  an  overview  of  the  data  acquisi¬ 
tion  and  analysis  requirements  for  the  Vibration  Facility  at  the  Wenner-Gren  Labo¬ 
ratory  .  The  equipment  avrilable  is  a  combination  of  analog  signal  conditioners 
and  digital  data  processors  including  a  dedicated  Raytheon  704  digital  computer 
(Appendix  b)  »  supplemented  by  a  remote  data  processing  sy3tem. 

1.2  The  data  processing  protocol  presented  here  is  designed  to  fulfill  the  re¬ 
quirements  of  measurement  of  the  physiological  responses  to  vibration  (AFOSR 
Contract  No.  F44620-69-C-0127) ,  and  to  Externally  Applied  Synchronous  Energy 
for  Cardiac  Support  (EASE,  N1H  Contract  No.  N01-H6-3-2928 . ) . 

1.3  These  specifications  are  intended  to  provide  the  base  for  initiating  the 
design  and  implementation  of  processing  programs  which  will  serve  to  satisfy  the 
data  storage,  on-line  (real-time),  and  off-line  processing  requirements  and  will 
use  existing  equipment  at  the  Wenner-Gren  Research  Laboratory  Vibration  Facility, 

1.4  This  document  does  not  imply  any  commitment  on  the  part  of  the  investigators 
to  implement  any  or  all  of  the  specifications  outlined. 
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INTRODUCTION 


2 . 1  Evaluation  of  Physiological  Stres s 

Determination  of  the  effects  of  stress  imposed  on  an  organism  may  be  accom¬ 
plished  by  measuring  a  number  of  physiological  variables.  Heart  rate,  blood 
pressures  and  flows,  respiration  rate,  oxygen  consumption,  and  body  temperature 
are  examples  of  appropriate  physiological  measurements.  Evaluating  relative 
changes  in  these  variables  caused  by  a  given  stress  requires  correlation  with 
characteristics  of  the  imposed  function.  For  the  particular  stress  of  vibra¬ 
tion  with  which  the  current  studies  are  concerned,  such  parameters  are  table 
displacement,  and  internal  organ  acceleration. 

2.1.1  Dogs  and  monkeys  chronically  implanted  with  flow  and  pressure  trans¬ 
ducers  are  used  routinely  as  experimental  subjects  in  vibration  stress  studies 
at  the  Wenner-Gren  Research  Laboratory.  The  animal  is  strapped  to  a  table  and 
accelerated  (Section  3.1)  in  either  the  vertical  (vibration-stress  case)  or 
horizontal  (vibration-stress  and  EASE)  direction.  The  table  iB  actuated  exter¬ 
nally,  and  moves  in  a  sinusoidal  pattern  at  various  frequencies  and  levels  of 
acceleration.  Calibration  of  all  physiological  transducers  is  accomplished  im¬ 
mediately  prior  to  a  run.  After  a  control  period  of  stable  values,  the  experi¬ 
mental  series  of  appropriate  design  is  conducted.  A  recovery  period  follows. 

In  all  case,  both  hard  copy  records  on  chart-paper  and  a  multi-channel  analog 
tape  recording  are  produced  continuously  during  the  experimental  period.  Digi¬ 
tal  data  processing  based  on  the  Raytheon  704  Computer  is  used  to  supplement  the 
analog  data  acquisition  facilities  with  real-time  (on-line)  reports,  and  to  pro¬ 
duce  an  edited  and  condensed  record  of  experimental  variables  on  digital  magnetic 
tape  for  permanent  storage  and  off-line  analysis.  All  acquisition  equipment  re¬ 
cords  a  digital  time  code  in  HH:MM:SS  format  for  later  comparison  ami  synchroni¬ 
zation  of  data.  The  typical  experimental  arrangement  is  shown  in  Fig.  2.1. 

2.1.2  The  Digital  Data  Acquisition  System  (DAS),  described  elsewhere  in  this 
report,  requires  a  remote  operator  control  panel  with  lights,  switches  and  a  tele- 
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type  to  control  the  DAS,  as  well  as  a  bank  of  appropriate  indicator  lights  which 
inform  the  operator  of  the  DAS  measurement  conditions.  A  data  coordinator  is 
also  present  during  an  experiment  to  inspect  measurements  as  they  are  acquired. 
2.2  Definition  of  Measurement  Interval 


The  experimental  protocol  is  based  on  a  four-component  measurment  interval 
which  may  be  repeated  a  maximum  of  50  times  and  will  not  exceed  6  hours.  The 
procedure  is  outlined  in  Fig.  2.2.  and  shows  the  three  components: 

A.  Control  period,  of  duration  t  seconds.  -  Defined  as  the  period  of  no 
time-varying  acceleration  preceding  a  test  period. 

B.  Test  period,  of  duration  t^  seconds.  -  Defined  as  a  period  where  the 
animal  is  submitted  to  steady  whole-body  acceleration. 

C.  Recovery  period,  of  duration  t^  seconds.  -  Begins  upon  termination 
of  the  test  period  and  ends  with  an  operator  command;  no  time-varying 
acceleration  during  this  period. 

D.  Pause  period,  of  duration  t^  seconds.  Data  are  not  recorded  on  digi¬ 
tal  magnetic  tape  during  this  period,  but  listing  on  the  teletype  continues. 

2.2.1  Duration  of  Test:  The  total  measurement  period  within  a  given  experi¬ 
ment  for  a  certain  frequency  and  waveform  will  be  defined  as  ttota^  ■  t£  +  tt  +  t^ 

t  seconds.  The  maximum  total  duration  of  a  single  test  interval  (t  )  will 
p  total 

not  exceed  390  minutes  in  any  circumstance;  the  minimum  duration  of  a  single  test 
interval  will  not  be  less  than  150  seconds. 

2.2.2  The  range  of  the  control  (t  ),  test  (t  ),  recovery  (t  )  and  pause  (t  ) 

c  t  r  p 

periods  is  specified  as  follows: 


IDENTIFICATION 


MINIMUM  DURATION 


MAXIMUM  DURATION 


(seconds) 


(seconds)  (minutes) 


21,600 


21,600 


360 


The  digital  deta  acquisition  facility  will  halt  measurements  during  the  off-line 
period  (as  deffred  in  the  System  Operation  Section,  paragraph  7.1). 

2.2.3  Because  the  measurement  interval  and  its  components  may  vary  over  a 
wide  range,  a  procedure  for  period  identification  is  necessary.  The  procedure 
must  be  flexible  and  applicable  at  the  time  data  are  acquired.  The  onset  of  the 
Recovery  Period  (t  )  will  be  identified  by  data  acquisition  system  (DAS)  detec¬ 
tion  of  the  shake-table  oscillator  output.  The  other  periods  will  be  identified 
by  signals  from  the  operator  control  panel  teletype  (Section  7.0). 


DATA  ACQUISITION  CHARACTERISTICS 


3 . 1  Vibration 

Both  Horizontal  and  Vertical  Shake  Tables  are  in  use  at  the  Wenner-Gren  Lab¬ 
oratory.  During  any  test  using  animals  with  chronically  implanted  flow,  pressure, 
and  other  transducers,  several  parameters  are  measured  by  the  Raytheon  704  Data 
Acquisition  Processing  System  (DAS) .  The  horizontal  shaks  table  is  designed  for 
use  with  the  Externally  Applied  Synchronous  Energy  (EASE)  experiments.  The  verti¬ 
cal  table  is  reserved  for  use  in  measurement  of  animal  response  to  whole-body 
vibration,  especially  with  regard  to  physiological  stress  measurements. 

3.1.1  The  tables  themselves  are  actuated  hydraulically  on  cue  from  external 
control  equipment.  The  displacement  of  the  table  is  determined  by  a  voltage 
function  generator.  The  signals  supplied  to  the  table  are  specified  to  be  pure 
:::.:e  waves  with  frequencies  between  2  and  30  Hz.  The  frequency  of  the  sine  wave 
will  be  entered  by  the  operator  and  is  sampled  by  the  DAS,  with  the  direction  and 
zero-crossings  of  the  shake  table  oscillatory  output  noted  for  use  in  calculating 
phasic  relationships.  A  period  of  lost  data  not  exceeding  2.5  seconds  at  the  start 
of  a  test  period  is  available  for  frequency  analysis. 

3.1.2  A  digital  time  code  is  acquired  by  the  DAS  and  is  placed  on  the  digital 
magnetic  tape  for  subsequent  identification  of  test  periods.  The  time  code  is  re¬ 
corded  in  BCD  format  as  HH:MM:SS.  Hardcopy  records  of  the  time-of-day  and  various 
average  parameters  are  emitted  periodically. 

3.1.3  Calibration  of  all  flow  and  pressure  transducers  is  accomplished  imme¬ 
diately  prior  to  and  following  the  experimental  procedure.  Periodic  rechecks  may 
be  revested  during  the  experiment.  These  calibration  results  are  input  to  the 
DAS  both  prior  to  and  during  an  experiment.  The  DAS  calculates  calibration  factors 
from  the  input  data  and  uses  the  factors  tor  calculations  of  the  physiological 


C'8 


variables . 


3 . 2  Signal  Processing 

3.2.1  All  physiological  signals  are  conditioned  by  various  analog  circuits 
before  presentation  to  the  DAS.  Insofar  as  possib.le,  random  noise,  drifting 
baselines.  Interference  from  external  electrical  equipment  and  intermodulation 
distortion  are  maintained  at  the  lowest  possible  values.  However,  flawless  data 
acquisition  i3  not  possible,  and  a  certain  amount  of  error  detection  by  the  DAS 
is  required.  Channels  which  present  illegal  data,  such  as  zero,  constant  or 
out-of-range  values  should  be  omitted  from  calculations,  after  emitting  an  error- 
code  and  channel  identification  to  tna  operator. 

3.2.2  The  physiological  and  vibration  parameters  of  interest  measured  by  the 
DAS  may  be  classified  as  follows: 

1.  Physiological 

A.  Cardiovascular 

1.  Pressure,  Left  Ventricular 

2.  Pressure,  Aortic 

3.  Pressure,  Arterial 

A.  Flow,  Aortic 

5.  Flow,  Coronary 

6.  Heart  Rate 
E.  Other 

1.  Respiration  Rate 

2.  Oxygen  Consumption 

3.  Body  Temperature 

II.  Vibration 

A.  Acceleration 

1.  Table 


2.  Body  Organs 


B.  Other 


1.  Table  Displacement 

2.  Net  Force 

III.  Time  Relationships  (between  groups  I  and  II). 

3.2.3  Based  on  their  data  acquisition  requirements,  the  parameters  of  interest 
can  be  grouped  into  three  categories;  (1)  within-beat  processing,  (2)  10-second 
averages,  and  (3)  time-related  variables.  Along  with  these  groups,  a  24-bit  time 
code  (HR:MIN:SEC)  in  BCD  format,  and  periodic  experimental  comments  in  ASCII  are 
entered  on  the  digital  tape. 

3.2.4  The  following  analog  signals  are  sampled  during  the  experimental  proce¬ 
dure  (not  all  signals  are  sampled  continuously  during  the  measurement  interval)  : 


VARIABLE 

CODE 

SAMPLING  RATE 
(msec) 

NO.  OF  SAMPLES 
per  sec. 

Left  Ventricular  Press. 

LVP 

2 

500 

Aortic  Pressure 

AP 

5 

200 

Arterial  Pressure 

ARP 

3 

200 

Aortic  Flow 

AF 

2 

500 

Consumption 

02 

100 

10 

Acceleration 

(Internal  Organs) 

ACC  I 

3 

200 

Net  Force 

NF 

5 

200 

Displacement  of  Table 

TD 

5 

200 

Mean  Coronary  Flow 

CFM 

500 

2 

Mean  Heart  Rate 

HRM 

500 

2 

Mean  Aortic  Pressure 

APM 

500 

2 

Mean  Arterial  Pressure 

ARPM 

500 

2 

Body  Temperature 

T 

1000 

1 

Table  Acceleration 

acct 

5 

200 

for  possible  later  addition: 
Venous  Pressure 

VP 

5 

200 

Venous  Flow 

VF 

5 

200 

Left  Ventricular 

Dimension 

LVD 

5 

200 

Respiration  Rate 

RR 

10 

100 

Peripheral  Arterial 

Flow 

PF 

5 

200 

3-axis  acceleration  (Internal  Organ  ACCI^ ,  ACCI^,  ACCl^) 
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3 . 3  Anaiog-to-Digltal  Conversion 


The  Raytheon  704  DAS  incorporates  a  24-Channel  "Milliverter"  (Model  MADC  10-06) 
analog-to-digital  converter.  This  device  contains  a  24  single-ended  input  multi¬ 
plexer  with  sample-and-hold  amplifier.  Full-range  input  analog  signals  are  converted 
to  9  bits  plus  sign  at  a  50  kHz  rate.  Conversion  time  is  9.8  usee,  aperture  time 
is  100  nsec,  and  approximate  through-put  rate  (single-channel)  is  50  kHz  (max). 

3.3.1  In  order  to  maximize  the  precision  in  all  measurements,  the  range  of 
signals  measured  by  the  Analog-to-Digital  converter  (ADC)  in  the  Raytheon  704  Sys¬ 
tem  is  extended  to  full  9  bit  significance.  In  other  words,  the  voltage  range  of 
these  signals  is  extended  to  ±10.0  volts.  A  full  significance  reading  by  the  ADC 
is  expected  to  include  1  part  in  512  accuracy,  plus  the  sign  of  the  input  signal. 

Thi9  is  equivalent  to  1.95  mv  in  10  v.  For  a  blood  pressure  calibration  of  100  mmHg/10 
V,  1  mmHg  -  100  mv,  and  resolution  would  be  on  the  order  of  ±0.2  ramHg.  For  i.  flow 
calibration  of  10  L/min/10  V,  resolution  would  approximate  ±20  ml /min. 

3.3.2  During  Pause  periods,  any  A/D  channel  may  be  selected  for  output  on  a 
D/A  channel  to  a  monitor  CRT  for  direct  visual  comparison  in  order  to  verify  satis¬ 
factory  A/D  operation. 
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DATA  ANALYSIS  REQUIREMENTS 
4 . 1  Within-Beat  Signal  Processing 

A  beat-by-beat  analysis  of  several  parameters  is  required.  These  para¬ 
meters  are  illustrated  in  Figs.  4.1.1,  4.1.2,  4.1.3  and  include: 


VARIABLE  CODE 

Left  Ventricular  Pressure  LVP 

Aortic  Pressure  AP 

Arterial  Pressure  ARP 

Aortic  Flow  AF 


Upon  receiving  an  external  trigger  signal,  the  within-beat  interval  will  be 
identified  and  the  DAS  will  begin  searching  for  valid  data.  The  maximum  sam¬ 
pling  rate  shown  in  section  3.2.4  will  continue  until  completion  of  the  beat 
is  identified  by  the  DAS,  at  which  time  sampling  will  cease  until  the  next  trig¬ 
ger  is  received.  The  data  will  then  be  condensed  into  ten  second  intervals.  For 
those  variables  designated  as  maximum  and  minimum  values,  single  points  (maximum 
and  minimum)  for  that  ten  second  interval  are  to  be  identified  and  stored  on  tape. 

For  each  within-beat  parameter  (i.e.,  peak  aortic  flow  (Q  ),  peak  d (LVP) ,  peak 

S  dt 

d  (AF)  ,  etc.)  which  is  not  averaged  before  input  to  the  DAS,  values  for  all  beats 
dt 

within  the  ten-second  interval  are  averaged  and  this  mean  value  is  stored  on  mag¬ 
netic  tape. 

4.2  Processing  of  Averaged  or  Slowly-Changing  Values 

Preprocessed  averages  of  several  parameters  are  acquired.  The  variables  are 
to  be  sampled  aa  Indicated  and  averaged  over  Intervals  of  some  multiple  of  10  sec. 
An  option  is  to  be  available  for  reporting  these  10  sec.  averages  either  individ¬ 
ually  or  as  averaged  10  second  blocks  in  integer  multiples  up  to  30.  But  always 
data  recorded  on  digital  magnetic  tape  will  be  in  10  sec.  blocks. 
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Figure  4.1.2  Aortic  Pressure  (AP)  and  Arterial  Pressure  (ArP) 
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The  averaged  parameters  include: 


VARIABLE  CODE 

Mean  Heart  Rate  HPM 

Temperature  T 

Oxygen  Consumption  02 

Mean  Aortic  Pressure  APM 

Mean  Coronary  Flow  CFM 

Mean  Arterial  Pressure  ARPM 


4 . 3  Tine  Relationship  of  Analog  Signals 

An  analysis  of  the  time  relationship  between  the  physiological  variables 
(AP,  LVP,  AF)  and  the  shake  table  variables  (displacement,  organ  acceleration, 
table  acceleration,  transmitted  force)  will  be  required.  The  term  0^  will  be 
used  to  define  the  relationship  of  the  vibration  table  to  the  cardiac  cycle. 
Specifically,  0T  is  defined  as  the  time  from  the  null  position  of  the  table 
to  some  event  of  interest.  An  example  is  illustrated  in  Fig.  4.3.1.  For  ex¬ 
ample,  defines  the  elapsed  time  from  the  passage  cf  the  table  through 

zero  to  the  initiation  of  aortic  flow  for  the  next  heart  bear. 

4 . 4  Phasic  Variables 

Parameters  of  the  vibration  table  and  acceleration  of  internal  body  organs 
resulting  from  the  vibration  inputare  required.  These  parameters  include: 
VARIABLE  CODE 

Table  Acceleration  Gt 

Organ  Acceleration 

Net  Force  F^ 

Table  Displacement  W 

Examples  of  these  parameters  are  shown  in  Figs.  4.4.1,  4.4.2,  4.4.3,  and  4.4.4. 
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Figure  4.4.2  Organ  Acceleration  (G  ) 


Figure  4.4.3  Net  Force  (F  ) 


4.5  Other  Parameters  of  Interest 


As  time  permits,  the  following  variables  should  be  measured  by  the  DAS, 
and  recorded  on  digital  magnetic  tape.  Digital  processing  of  these  data  1b 
of  low  priority  and  should  be  accomplished  only  if  no  interference  with  the 
aforementioned  measurements  can  occur. 


VARIABLES  CODE 

Venous  Pressure  VP 

Venous  Flow  VF 

Left  Ventricular  Dimension  LVU 

Respiration  Rate  RR 

Peripheral  Arterial  Flow  PF 

4.6  Calibration 


A  calibration  table  will  be  constructed  by  reading  standard  signals  from  each 
channel  under  operator  control  at  the  beginning  of  the  experiment.  All  channels 
will  be  linear,  so  two  values  will  be  required  for  each  charnel,  i.e.  y  -  mx  +  b; 
m  and  b  will  be  calculated  by  the  calibration  program  and  stored  in  a  table.  It 
should  be  possible  to  re-calculate  the  calibration  for  a  channel  during  the  experi¬ 
ment  (during  a  pause  or  'stop*  period)  and  change  the  calibration  table. 
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DATA  STORAGE  REQUIREMENTS 


5.1  Magnetic  Tape  Storage 

Data  acquired  during  an  experimental  procedure  are  stored  on  both  analog 
tape  and  digital  tape.  All  raw  data  are  stored  in  unmodified  form  on  analog 
tape  along  with  an  IRIG-standard  time-code  superimposed  on  a  1  kHz  carrier. 

This  time  code  may  be  converted  into  HH:MM:SS  format  by  the  time  code  reader, 
permitting  the  off-line  retrieval  of  raw  data  a  a  later  time. 

5.1.1  Data  stored  on  digital  magnetic  ta;.  will  have  undergone  sampling 
and  considerable  data  compression.  An  identification  recc’-d  incorporating  the 
following  is  required  at  the  beginning  and  end  of  eve.-_  rest  tape  on  each  animal: 

1 .  Animal  number 

2.  Date:  MM/DD/YY 

3.  Comments 

The  header  for  each  data  block  will  include: 

1.  Time  of  Day;  HH:MM:SS 

2.  Vibration  Frequency,  (Hz) 

3.  G-Level 

4.  Period:  Control,  test,  or  recovery. 

Results  from  no  more  than  one  experimental  animal  are  to  be  stored  on  a  single 
digital  magnetic  tape.  Data  retrieval  at  a  later  time  may  be  accomplished  by 
entering  the  start  and  end  time  in  HH;MM:SS  whereupon  the  entire  contents  of 
the  tape  between  the  two  times  are  dumped  on  the  teletype. 

5.1.2  As  required  by  the  epecific  data  presentation  demands  explained  in 
Section  6.0  of  this  report,  CRT  display  of  the  physiological  data  may  be  accom¬ 
plished  off-line,  i.e.  during  pause  and  atop  periods. 

5.1.3  All  digital  magnetic  tape  files  for  a  specific  experiment  are  to  be 
kept  on  separate  tapes,  one  animal  per  tape.  Storage  is  accomplished  in  1BM- 
compatabie,  9  truck  SGO  bpi  format  with  HC  header  or  trailer  labels.  The  end 
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of  a  data  run  is  delinated  by  3  file  marks  on  the  tape.  Odd  parity  and  binary 
format  for  numbers  are  required.  Integer  values  for  physiological  data  are  used 
wherever  practical. 

5.2  Intermediate  Disk  Storage 

Buffer  area  on  the  fixed-head  disk  is  required  for  on-line  analysis  and 
data  presentation.  This  area  should  be  of  sufficient  size  to  allow  graphic  pre¬ 
sentation  of  all  the  variables  listed  in  Section  6.  These  displays  may  be  re¬ 
quested  by  the  system  operator  during  pause  or  stop  periods,  for  a  "display  up¬ 
date"  or  at  the  end  of  the  run  for  a  "summary  display".  The  formats  for  these 
display  patterns  are  explained  in  section  6  of  this  report. 


DATA  AND  RESULTS  PRESENTATION 


6.1  Teletype  Formats 

In  the  data  acquisition  protocol  described  here,  two  teletypes  (KSR33) 
are  required,  i.e.,  teletype  //I  for  a  chronolog  of  averaged  physiological 
data,  and  teletype  112  for  entering  operator  comment  statements  regarding 
the  progress  of  the  experiment.  Teletype  #1  output  is  required  for  an  on-line 
record  of  integer  multiples  of  10  second  average  values  of  the  physiological 
parameters.  This  record  is  used  to  assess  the  present  state  of  the  experiment 
and  to  review  data  from  the  beginning  of  the  experiment.  The  format  for  print¬ 
ing  these  data  is  pregented  in  Fig.  6.1.1.  Teletype  #2  is  used  to  type  and  pre¬ 
view  comments  before  they  are  entered  on  the  magnetic  tape  and  printed  output 
of  teletype  ill. 

6.2  CRT  Graphic  Display  Output 

A  Tektronix  611  storage  and  display  oscilloscope  is  available  for  use  by 
the  DAS.  The  CRT  display  may  be  requested  in  three  possible  formats.  Para¬ 
meters  of  interest  may  be  plotted  against  (1)  table  frequency,  (2)  test  inter¬ 
val,  or  (3)  time  relationship ,  An  example  of  each  with  the  parameters  of 

interest  is  shown  in  Fig.  6.2.1.  At  any  one  time  up  tc  four  dependent  variables 
may  be  plotted  against  the  same  independent  variable.  A  typical  display  is  illus¬ 
trated  in  Fig.  6.2.2.  An  ''update"  display  is  defined  as  a  graphic  plot  vhich  can 
be  viewed  as  the  data  points  are  added  at  the  completion  of  each  test  period. 

A  "summary  display"  is  defined  as  a  plot  with  all  data  points  for  the  experiment 
shown,  and  is  available  only  in  off-line  mode  at  the  end  of  a  complete  experiment. 
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Figure  6.2.2  Example  of  611  Update  and/or  Summary  Display 


OPERATION  OF  THE  SYSTEM 


.0  The  data  acquisition  system  will  require  a  control  console  for  interfacing 
between  1)  Data  Coordinator  and  2)  System  Operator  during  the  experiment. 

The  System  Operator  is  seated  at  the  main  console  and  is  responsible  for  con¬ 
trolling  and  editing  data  sent  to  the  DAS.  This  is  accomplished  by  input  from 
Teletype  No.  2  (TTY2).  The  Coordinator  also  uses  TTY2  for  calibration  inputs 
(/CaL),  comments  (/MSG),  and  judges  the  quality  of  data  and  results  by 
viewing  analog  signals.  Both  the  System  Operator  and  Coordinator  must 
communicate  in  order  to  define  the  various  measurement  periods.  A  layout 
of  the  control  console  is  presented  in  Fig.  7.0.1.  It  has  three  major  di¬ 
visions,  (1)  control  of  data  acquisition,  (2)  hardcopy  10  sec.  (or  an  in¬ 
teger  multiple  of  sec.)  averages  of  data  on  teletype  No.  1  and  hardcopy  ana¬ 
log  signals  of  physiological  data,  and  (3)  CRT  graphic  displays  for  data  in¬ 
spection.  The  layout  of  the  Control  Iniicator  light  panel  in  the  control 
console  is  shown  in  Fig.  7.0.2. 

7 • 1  Operational  Procedure 

A.  The  System  Operator  turns  on  the  Control  indicator  light  panel  by 
pressing  the  'ON'  panel  button.  He  also  enters  through  TTY 2  the  information 
for  the  identi rication  record  described  in  section  5.1.1. 

B.  The  Coordinator  inspects  the  incoming  analog  data.  When  these  are 
acceptable,  he  announces  readiness  to  begin  the  experiment. 

C.  The  System  Operator  then  enters  (ON  TTY2)  the  appropriate  frequency 
(/FRQ)  and  G-level,  "Begin  Control  Period"  (/CON),  which  initiates  data  acquisi¬ 
tion  for  transfer  to  digital  magnetic  tape  and  TTY1  output.  This  action  also 
initiates  a  timer  to  time  the  control  period.  (See  Table  7.1.1  for  TTY  commands). 

D.  After  consultation  with  the  Coordinator,  the  control  period  is  terminated 
by  entering  "End  Period"  (/ENO)  and  "Begin  Test  Period"  (/ TST)  code  on  the  TTY2, 

This  resets  the  timer  to  zero,  cues  the  shake  table  operator  to  begin  vibration,  and 
stops  transfer  of  aata  to  the  digital  magnetic  tape.  After  steady  state  values 
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Figure  7.0.1  Control  Con3ole 


Figure  7.0.2  Control  Indicator  Light  Panel  of  Control  Console 


TABLE  7.1.1 


TELETYPE  COMMANDS 

/CXL  Cancel  Present  Data  Period 

/CAL  Calibrate  Analog  Channels 

/CON  Begin  Control  Period 

/END  End  Period 

/RCV  Begin  Recovery  Period 

/RST  Begin  Test  Period 

/ EOT  Write  End  of  File  on  Tape 

/BVE  Terminate  Experimental  Procedure 

/MSG  Next  Line  Entered  Should  go  to  Comments 

Part  of  MAG  Tape  Record 
/'FRQ  Enter  Vibration  Frequency 

/GLV  Enter  G-Lovel 

*  Rubout  hey  cancels  command  just  entered 

*  Carriage  Return  and  Line  Feed  Initiates  Command. 


1 


4 
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of  vibration  have  been  reached,  the  Shake  Table  Operator  presses  his  Table 
Ready  button.  This  turns  on  the  Table  Ready  light  on  the  Control  Console  and 
again  resets  the  timer  to  zero  and  data  transfer  to  the  digital  magnetic  tape 
is  resumed. 

E.  The  Vibration  Test  period  is  terminated  by  typ'  ig  (/END)  and  "Begin 
Recovery  Period"  (/RCV)  by  the  System  Operator.  The  (/ RCV)  signal  1)  cues  the 
Shake  Table  Operator  to  stop  vibration;  2)resets  the  timer;  and  3)automatically 
initiates  the  recovery  sequence.  Sequence  (3)  consists  of  a  DAS  search  for 
zero  output  from  the  oscillator  driving  table.  When  the  output  reaches  zero, 
the  DAS  resets  the  panel  timer  and  activates  the  "Recovery"  panel  light. 

F.  Recovery  period  can  be  terminated  by  typing  (/END).  This  stops  trans¬ 
fer  of  data  to  digital  magnetic  tape  (TTY  1  output  continues).  CRT  611  is  now 
accessible  for  update.  Any  other  option  is  available  after  this  button  is  acti¬ 
vated.  For  example,  (/BYE)  Inactivates  all  digital  data  acquisition,  with  CRT 
611  remaining  accessible  for  update.  Or,  enter  f  (/FRQ) ,  G  (/CIV),  and  (/CON) 
to  repeat  above  sequence. 

G.  Typing  (/CAL)  Indicates  to  the  DAS  a  calibration  operation  is  to  be 
performed.  The  operation  will  be  carried  out  by  interaction  with  TTY2. 

H.  A  "Cancel  Present  Data  Record  Period"  (/CXL)  will  also  be  available, 
This  label  indicates  the  data  from  the  preceding  period  is  invalid  and  is  not 
to  be  used  in  the  standard  data  reduction  routine.  For  any  other  corrections, 
tne  TTY2  comments  will  be  used.  The  (/CXL)  command  also  enables  restart  of 
the  previous  period. 

I.  The  Coordinator  may  enter  comments  via  T1Y2  at  any  time  during  the 
experimental  period. 

7 . 2  Displays 

The  System  Operator  may  choose  the  desired  plots  of  data  defined  in  Section 
6.2.  The  611  display  is  initiated  by  (1)  pressing  "PLOT",  (2)  pressing  up  to 
four  variables  for  Y  and  one  for  X,  and  (3)  depressing  "ENTER". 
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This  paper  describes  e  novel  approach  to  quantifying  the  coefficients  of  a  bio- 
vibrational  model.  Unanesthetized  Rhesus  monkeys  were  exposed  to  vertical 
sinusoidal  vibrations  at  discrete  frequencies  from  2  to  JO  Hz  at  0.5  and  1.0  g  ac¬ 
celeration  amplitudes.  With  the  subjects  in  a  sitting  position  whole  body  mechan- 
leal  velocity  impedance  was  ob’alned  from  measurements  of  vibration  exciter 
velocity,  transmitted  force  and  the  phase  relationship  between  velocity  and  force. 
The  data  obtained  were  applied  to  a  two-moss,  single-degree  of-freedom  model 
to  determine  the  coefficients  of  elasticity  and  damping  as  a  function  of  frequency 
and  body  mass,  The  coefficients  were  obtained  for  each  test  it  each  frequency 
by  requiring  the  Impedance  of  the  model  to  match  that  of  the  primate.  A  part  of 
the  model  Is  on  empirical  equation  defining,  as  a  function  of  vibration  frequency 
and  Intensity,  the  ratio  of  nonreactlve  to  reactive  animal  mass.  After  normalizing 
for  body  mass  and  averaging,  function*,  for  the  elastic  and  damping  coefficients 
were  developed  such  that  the  impedance  response  of  the  model  closely  matched 
♦he  experimental  results  throughout  the  ?  to  30  Hz  range  at  0.5  g  and  also  at  1.0  g. 
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A  Model  to  Predict  the  Mechanical 
Impedance  of  the  Sitting  Primate  During 
Sinusoidal  Vibration 

R.  a  COWARDS  J.  F.  LAFFERTY 


High-speed  vehicles  frequently  transmit  dy¬ 
namic  forces  to  their  occupants.  Depending  upon 
the  intensity  and  duration  of  such  disturbances, 
serious  lmpalrmer.c  of  operator  or  passenger  func¬ 
tioning  may  occur.  This  problem  has  led  to  ex¬ 
tensive  research  directed  toward  defining  and 
understanding  the  dynamic  response  of  the  human 
body.  Laboratory  testing  of  man  at  l,w  accelera¬ 
tion  amplitudes  has  demonstrated  that  whole  body 
vibration  In  the  2-  to  30-HZ  frequency  region 
produces  the  most  pronounced  mechanical  effects 
for  iny  given  acceleration  amplitude.  Ziegen- 
ruecker  and  Magld  reported  the  Intensity 

tolerance  of  sitting  man  to  vertical  vibration 
to  be  lowest  In  the  4-  to  8-Hz  frequency  range. 
Coermann  (2)  recorded  the  force  transmitted  to 
both  sitting  and  standing  man  and  found  a  whole 
body  primary  resonant  frequency  of  4  to  6  Hz. 

He  also  demonstrated  a  remarkably  good  approxima¬ 
tion  to  this  experimental  data  by  use  of  a  simple 
linear  sprlng-mass-damper  model  of  the  human 
body.  Other  Investigators  ‘-av*  confirmed 

these  result*  and  contributed  add'  nonal  informa¬ 
tion  to  help  define  man's  vibration  resoonse. 
Although  high  acceleration  Intensities  are  at 
times  encountered  m  vehicles  and  machinery,  man 
cannot,  for  safety  reasons,  serve  an  a  laboratory 
subject  In  high  intensity  vibration  tests.  Xr- 
such  cases,  animals  have  typliiliy  been  used  as 
human  surrogates  ( 9-12 ) . 

Analytical  models  to  predict  the  dynamic 
response  of  man  and  animals  have  proven  useful  in 
many  Instances  (_2  and  J)  •  This  paper  describes 
ouen  a  model  to  accurately  predict  the  Impedance 
response  for  the  Bitting  Rhesus  morucey  when  vi¬ 
brated  vertically  (l.e.,  +  ge) • 

ANALYSIS 

For  the  system  shown  In  Fig.  1,  whole  tody 
mechanic*’,  veloolty  Impedance  (hereafter  termed 


1  Underlined  numbers  In  parenthesis  assign*!* 
Reference*  at  end  of  paper. 


simply  "impedance")  Is  defined  as  the  complex 
ratio  of  the  Interface  (transmitted)  force  to 
the  velocity  at  the  point  of  force  transmission 
(the  vibration  exciter  velocity  In  the  case  of 
whole  body  impedance). 


l  - 
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where 


Z 

V 

* 


■  complex  Impedance 
•  amplitude  of  the  transmitted  force 
-  vibration  exciter  velocity  amplitude 
»  phase  angle  relationship  of  the  transmitted 
force  relative  to  the  vibration  exciter 
ve loc lty . 


Coermann  (2)  has  shown  that  for  a  syBtem 
consisting  of  one  mass  (m)  atop  a  linear  spring 
(of  spring  constant,  or  coefficient  of  elasticity, 
k)  and  linear  damper  (of  damping  constant,  or 
coefficient  of  damping,  c)  In  parallel,  the  fol¬ 
lowing  equations  define  the  Impedances 


2 


k2  ♦  a2^  ♦  c2w^ 


2lu*»5 


1/2 
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(2) 


where  w  -  2rf  and  f  Is  the  vibration  frequency. 

He  illustrated  that  the  impedance  of  sitting  and 
standing  subjects  could  bo  approximated  by  that 
of  tne  model  defined  by  (1)  and  (2)  with  appropri¬ 
ate  values  of  whole  body  demplng  and  elastic 
coefficients.  It  Bhould  be  noted  that  this  model 
assumes  linearity  of  cc  ;ff  dents,  and  thus  the 
Impedance  is  independerc  of  vibration  exoiter  dis¬ 
placement  (end  therefore  acceleration  amplitude). 

Tha  model  we  shown  In  Fig.  2  evolved  from 
ills  ErsCtrtcn-Von  Olerk*  (9)  model  (Pig-  1)  and 
*rom  experimental  observation*  of  a  number  of 
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Pig.  1  Development  of  a  slngle-degree-of-fraedom 
model  for  impedance  response  of  the  restrained 
primate.  Prom  Broderson  and  Von  Olerke  (^) 


vibration  tests  with  Rhosus  monkeys.  The  Imped* 
ance  modulus  end  phase  angle  equations  for  this 
model  can  be  shown  to  be  (12): 


fl  ♦  r)*[fc*  ♦  [ho»,(l  ♦  t)  - 

ft2  ♦  «2  (e,2*2  ♦  c2  ■  lk»,) 

.1  I  (1  ♦  i)(iV  ♦  ft2)  ♦  a,«r2[etf«J  •  k(l  ♦ 

'  - , - 

I  <a<« 


1 1  n 
(3) 

(**) 


where  r  ■  the  ratio  of  "Inert"  to  "sprung" 

mast . 

The  development  of  equations  (3)  and  (4) 
require  that  k  and  c  be  Independent  of  amplitude, 
but  permit  a  nonlinear  frequency  dependence.  A 
nonlinear  variation  of  k  and  c  with  frequency 
would  Imply  that  system  strain  Is  sensitive  to 
the  rate  of  loading.  While  such  a  response  is 
not  uncommon  for  biological  materials,  a  frequency 
dependence  of  k  and  o  may  depend,  In  part,  on 
the  nature  of  the  model  Itself. 

The  function,  r,  Is  expected  to  vary  with 
both  frequency  and  acceleration  amplitude.  Body 
subsystems  may  well  be  amplitude  restricted  rela¬ 
tive  to  whole-body  motion,  l.e.,  the  displace¬ 
ment  of  some  subsystems  are  United  by  the  pres¬ 
ence  of  other  body  components.  Subsystems  which 
reach  their  relative  displacement  limits  will 
tend  to  Increase  the  Inert  mass  response  of  the 
total  system  as  deceleration  amplitude  is  further 
Increased;  however,  amplitudes  )n  excess  of  l.C 
g  may  well  result  In  a  more  c  implicated  amplitude 
dependence  of  r.  At  constant  acceleration  ampli¬ 
tude,  while-body  response  is  typically  character¬ 
ized  (2-B)  as  having  a  relatively  large  Inert 
mass  component  at  low  frequencies,  while  an  In¬ 
creased  damping  response  at  higher  frequencies 
Indicates  9  r‘>’"t,vn  decrease  In  the  Inert  mass 
component.  Tne  ratio,  r,  might  then  be  expected 
to  Increase  with  In.ivasl  acceleration  arr.pll- 
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Pig,  2  A  two-mass,  slngle-degree-of-freedom 
model  for  Impedance  response  of  the  sitting,  re¬ 
strained  primate 


tude  and  decreasing  frequency. 

By  requiring  that  k  and  c  be  Independent 
of  amplitude,  the  amplitude  dependence  of  r  can 
be  determined  from  Impedance  data  obtained  for 
carious  amplitudes  at  constant  frequency.  Suc¬ 
cessive  iterations  at  each  frequency  resulted 
In  an  expression  for  r  in  the  form 


where  A  is  tne  acceleration  amplitude  expressed 
In  units  of  g. 

Equations  (3),  (4),  end  (5)  define  the  vi¬ 
bration  response  of  the  two-mase  model.  The  ex¬ 
periments  and  results  dlscusse.  In  the  following 
sections  (a)  establish  the  experimental  impedance 
response  of  the  Rhesus  monkey  vibrated  from  2 
to  30  Hz  at  0.5  and  1.0  g  acceleration  amplitudes, 
(b)  establish,  based  upon  the  ..lodel  Just  presented, 
whole-body  elastic  and  damping  coefficients  as  a 
function  of  frequency,  and  (c)  compare  the  re¬ 
sponse  of  tne  model  to  that  experimentally  mea¬ 
sure  . 

METHODS 

The  We: ner-Oren  Research  Laboratory's 
vertical  electrohydraulic  vibration  exciter  ( 13 . 

14)  was  utilized  In  this  study.  A  force  cell 
arrangement  (12_,  14)  allowed  direct  recording  of 
the  force  transmitted  between  vibration  exciter 
and  subject.  Vibration  exciter  displacement  and 
velocity  were  recorded  from  the  output  of  recti¬ 
linear  potentiometers  attached  to  the  exciter. 

By  supplying  vibration  exciter  velocity  and  net 
transmitted  force  tb  a  phase  angle  meter,  an  analog 
signal  was  produced  which  was  proportional  to 
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Pig.  3  Primate  positioned  atop  vertical  vibration 
exciter 


Tabla  1  Summary  of  Animals  Tested  and 
Vibration  tarameters 


r — - 

ANIMAL 

(WItGKT) 

ACCtUJUTION 

AMFLIUjO* 

1—  - 

of  r.xpcjvui 

(MiwlCMi  Ord*r  >  to 

2.  J.  A,  5.  b.  1.  1.  *. 

10,  n.  1J.  U,  11.  30. 

)0  H« 

A  <J3  lb.' 

0.5  1 

3 

t  (Jft  lb.) 

0  5  « 

3 

C  (10. 3  lb.) 

0.5  , 

3 

D  Ml. 3  lb  ) 

0.5  , 

3 

C  (10. 3  kb.) 

1.0  • 

3 

0  (11. 9  lb.) 

1.0  • 

3 

the  phase  difference  of  the  two  Input  signals. 
Mechanical  velocity  Impedance  was  thus  obtained 
by  continuously  recording  (a)  net  transmitted 
fovea ,  (b)  vibration  exciter  velocity,  and  (o) 
tne  phase  relationship  between  equations  (1)  and 
(?).  pig.  3  1#  a  photograph  of  the  experimental 
setup. 
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Pig.  4  Complex  Impedance  verBuo  frequency  for 
0-5-g  vibration  of  primates  A,  B,  C,  and  D 

Tho  primates  were,  exposed  to  a  series  of 
frequencies  (see  Table  1),  e:  ch  at  the  same  ac¬ 
celeration  amplitude.  Pour  animate  wore  so  tested 
five  times  each  at  an  acceleration  amplitude  of 
0.5  g.  The  two  emaller  animals  (C  and  D)  were 
then  each  tested  five  times  at  an  acceleration 
amplitude  of  1,0  g.  An  exposure  duration  of  30 
sec  at  each  condition  was  specified  on  the  basis 
of  previous  tests  which  established  30  sec  to  be 
sufficient  to  obtain  reliable  results.  Extended 
exposures  also  introduce  temporal  effects.  Tests 
on  any  giver,  animal  were  separated  by  at  least 
one  day. 

The  monkeys  required  no  sedation  or  other 
medication  for  these  experiments  since  they  were 
well  adapted  to  the  experimental  environment. 

In  addition  to  dally  handling,  the  subjects  were 
previously  submitted  to  repeated  exposures  of 
extended  Duration  at  all  frequencies  including 
6-hr  exposures  at  resonance.  When  placed  In  the 
Impedance  chair  the  animals  were  securely  restrain¬ 
ed  by  a  lap  belt  and  by  belts  around  ea<  :>  ankle. 

Tne  subject  was  restrained  from  bending  ty  a 
rigiu,  mein!  i.cv;k  pises  which  physically  mntaoted 
the  animal  only  if  he  attempted  to  bend  out  of 
the  upright  position. 
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Pig,  5  Impedance  phase  angle  and  normalized  modu¬ 
lus  versus  frequency  for  0.5-g  vibration  of  pri¬ 
mates  A,  B,  C,  and  D 


RBSULTS 

The  oomplex  lmpedenae  for  the  0.5  g  experi¬ 
ments  Is  presented  In  Pig.  4  as  a  function  of  vi¬ 
bration  frequency.  The  data  from  Subjects  A  (25 
lb),  and  B  (26  lb),  are  grouped  together,  as  are 
those  from  Subjects  C  (10. 5  lb),  and  D  (11.5  lb), 
since  the  animals  In  each  group  differed  In 
weight  by  only  1  lb;  thus,  two  curves  are  shown 
for  both  Impedance  modulus  and  phase  angle. 

Since  eaoh  animal  was  tested  five  times  at  0.5  g, 
each  curve  In  this  figure  represents  the  mean  of 
10  tests.  The  limits  shown  correspond  to  the  mean 
value  plus  or  minus  one  standard  deviation.  The 
two  straight  lines  designated  as  "mu"  In  the  modu¬ 
lus  graph  correspond  to  tht  Impedance  magnitude 
of  an  "inert"  mass  of  weight  equal  to  that  of  the 
average  for  the  eubjects.  A  peak  Impedance  oc¬ 
curs  at  7  to  8  Hz  for  each  group.  At  lesser  fre¬ 
quencies  the  impedance  moduli  follow  closely 
those  of  an  Inert  mass  of  like  weight.  At  fre¬ 
quencies  immediately  greater  than  thos  1  of  the 
moduli  peaks  the  Impedance  magnitudes  decrease. 
Prom  12  to  30  Hz  the  moduli  remain  ab< ut  constant. 
The  phase  angle  curves  are  approximat >ly  90  deg 
(the  Impedance  phase  angle  of  an  Inert  mess)  at 
2  Hz,  and  then  decrease  rapidly  In  the  "Impedance 
moduli  peak"  region.  Beyond  about  U  Hz  the 
curves  are  In  the  10-  to  20-deg  ran/e.  This 
could  Indicate  a  predominant  damping  response 
since  the  Impedance  phase  angle  of  a  mass-less 
damper  la  0  deg. 

The  moduli  data  shown  In  PI/,.  4  were  nor¬ 
malized  by  body  weight  and  plotted  In  Pig.  5. 
versus  vibration  frequency.  Sln/e  the  phase 


Pig.  6  Normalized  whole  body  coefficients  of 
elasticity  and  damping  versus  frequency  for  the 
sitting,  restrained  primate 

angle  curves  as  shown  in  Pig.  4  did  not  appreci¬ 
ably  differ  between  the  different  size  animals, 
these  data,  for  each  frequency,  were  combined 
without  normalization.  This  figure  thus  repre¬ 
sents  the  complex  impedance  response  cf  any  size 
Rhesus  monkey,  In  the  sitting  position,  when  vi¬ 
brated  vertically  at  0.5  g  from  2  to  30  Hz. 

The  curves  or  Pig.  5  represent  mean  values 
as  determined  at  each  tost  frequency  by  twenty 
data  points.  For  each  frequency  20  complex  Im¬ 
pedance  values  were  used  In  equations  (3).  (4), 
and  (5)  to  solve  for  k  and  c.  Therefore,  for 
each  of  the  15  different  test  frequencies,  20 
sets  of  k  and  c  were  determined.  Each  set  repre¬ 
sents  the  magnitude  of  elastic  and  damping  coef¬ 
ficients  whim,  when  applied  to  the  model  of 
Pig.  2,  resulted  In  a  complex  Impedance  response 
which  matched  that  experimentally  measured  at 
that  frequency.  The  average  values  of  k  and  c 
thus  determined  ut  each  test  frequency  are  pre¬ 
sented  In  Fig.  6.  The  mean  elastic  coefficient 
varies  approximately  directly  with  vibration 
frequency.  On  the  other  hand,  the  mean  damping 
coefficients  exhibit  a  curve  which  looks  similar 
In  shape  to  the  impedance  modulus  plot  of  Fig. 

5;  however,  the  peak  of  the  damping  coefflclen: 
curve  Is  at  5  Hz  rather  than  at  7  to  8  Hz. 

The  curves  of  Fig.  6,  along  with  equation 
(5),  completely  quantify  the  model  presented  In 
Pig.  2.  For  computational  purposes,  values  cf 
k  and  c  from  Fig.  6  were  approxlnated  by  three 
linear  segments,  as  defined  In  Pig.  2,  for  the  2- 
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Fig.  7  Impedance  phase  angle  and  normalized  Fig.  6  Impedance  phase  angle  and  normalized 

modulus  versus  frequency  for  0.5-g  vibration  of  modulus  versus  frequency  for  l.O-g  vibration  of 
primates  A,  B,  C.  and  D  compared  to  that  pre-  primates  C  and  D  compared  to  that  predicted  by 

dieted  by  the  model  of  Pig.  2  the  model  of  Fig.  2 


to  5-Hz,  5*  to  10-Hz,  and  10-  to  30-Hz  regions. 

Complex  impeuance  obtained  from  tne  model 
is  compared  in  ?lg.  7  ulth  experimental  results 
for  an  acceleration  amplitude  of  0-5  g.  Inasmuch 
as  these  experimental  data  were  used  to  obtain 
the  frequency  dependence  of  k  and  c,  the  excellent 
fit  to  the  data  verifies  the  validity  of  the  com¬ 
putational  procedures. 

The  generality  of  the  model  is  demonstrated 
in  Pig.  8,  where  the  calculated  impedance  for 
1.0  g  acceleration  amplitude  If  compared  to  the 
experimental  results  from  Subjects  C  and  D.  The 
variation  of  the  modulus  as  predicted  by  the  model 
provides  a  very  good  approximation  to  the  experi¬ 
mental  results  in  both  magnitude  and  6nape.  The 
impedance  phase  angle  from  the  model  generally 
falls  within  tire  experimentally  measured  llmitB 
for  frequencies  up  to  14  Hz.  For  14  to  30  Hz 
tne  phase  angle  predicted  by  the  model  is  general¬ 
ly  about  10  deg  higher  than  the  upper  limit  of  the 
experimentally  determined  values. 

DISCUSSION 

The  object  of  this  investigation  was  to  es¬ 
tablish  experimentally  the  Impedance  response  of 
the  sitting  Rhesus  monkey  during  0.5  and  1.0  g 
vibration  at  2  to  30  Hz,  and  to  develop  a  model 
which  would  predict  an  Impedance  closely  matching 
that  experimentally  measured.  The  Impedance 
response  of  the  model  (based  upon  the  experimental¬ 
ly  derived  coefficients  presented  in  rig.  C)  vary 
closely  matches  that  experimentally  measured 
during  0.5  and  1.0  g  vertical  vibration  (PlgB. 

7J-* 


7  and  t) . 

The  model  of  Pig.  2  and  lte  manner  of  deriva¬ 
tion  are  somewhat  novel.  Although  the  model  is 
composed  of  linear  elements,  the  elastic  and 
damping  coefficients  are  defined  to  be  functions 
of  vibration  frequency.  In  addition,  the  ratio 
of  nonreactlve  (inert)  to  reactive  mass  in  the 
model  is  defined  to  vary  directly  with  input  ac¬ 
celeration  amplitude,  and  Inversely  ulth  the 
square  root  of  vibration  frequency.  Thus  for  any 
given  acceleration  intensity,  most  of  the  total 
mass  Is  inert  at  the  lower  frequencies,  and  be¬ 
comes  more  reactive,  or  "sprung,"  as  frequency  in¬ 
creases.  This  is  consistent  ulth  previous  experi¬ 
mental  measurements  which  show  the  subject  to 
exhibit  a  "mass-like"  response  at  lower  freauencle6 
and  a  predominant  damping  response  at  higher  fre¬ 
quencies.  Por  any  given  frequency.  Increasing 
the  acceleration  amplitude  increases  the  amount 
of  inert,  relative  to  total,  mass.  It  Is  postu¬ 
lated  that  increasing  the  input  intensity  could 
cause  "bottoming"  of  certain  suspended  Internal 
organs  (or  "parts")  and  that  such  a  system  would 
thus  respond  in  a  more  i.onreactive  manner  at  the 
higher  intensities. 

How  well  tne  model  responds  to  input  accel¬ 
erations  greater  than  1.0  g  remains  to  be  estab¬ 
lished.  It  is  interesting  to  compare  the  results 
of  Brcderson  and  Von  Qierke  (_2)  to  those  of  tho 
present  investigation.  A  whole  body  elastic  co¬ 
efficient  of  51.7  lb/in.  and  a  whole  body  damping 
coefficient  of  0.86  ib-sec/in.  were  reported  for 
a  l1-, 5  lb  Rhesus  monkey  vibrated  at  1.0  g.  These 
figures  were  extracted  by  matching  the  equation 


for  Impedance  modulus  from  tha  model  of  Pig.  1 
to  that  experimentally  established  during  5-  and 
12-Hz  vibration.  The  model  of  the  present  study, 
as  shown  In  Pig.  2,  yields  (for  a  14.5  lb  Rhesus) 
a  damping  coefficient  of  0.812  lb-sec/ln.  and 
an  elastic  coefficient  of  45.8  lb/in.  (using  a 
10-Hx  frequency).  A  favorable  comparison,  there¬ 
fore,  exists  between  these  two  models;  however, 
that  of  the  present  investigation,  being  derived 
from  both  impedance  modulus  and  phase  angle  data 
at  each  frequency  during  0.5  g  vibration,  yields 
an  exact  match  to  the  experimental  data  at  0.5  g, 
but  more  Importantly,  responds  to  changes  In  ac¬ 
celeration  In  a  manner  similar  to  those  experi¬ 
mentally  recordod. 

Future  experiments  are  now  planned  In  which 
Rhesus  monkeys  will  be  vibrated  at  higher  accelera¬ 
tion  amplitudes  for  short  periods  of  time  at  dis¬ 
crete  frequencies  from  2-  to  Temporal  ef¬ 

fects  due  to  extended  vibration  at  given  frequen¬ 
cies  will  also  be  determined,  and  hopefully  In¬ 
corporated  into  tha  model. 
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ABSTRACT  OF  THESIS 


ANALOG  COMPUTER  SIMULATION  OF  THE  RESPONSE  OF  THE 
CIRCULATORY  SYSTEM  TO  WHOLE-BODY  ACCELERATIONS 

The  closed-loop  response  of  the  human  circular::;.,  system  to 
whole-body  acceleration  patterns  of  both  constant  and  time-varying 
periodic  forms  was  simulated  in  this  study.  lilt  responses  of  the 
developed  model  were  first  determined.  In  a  head-up  tilt  cardiac  out¬ 
put  decreased,  and  in  a  head-down  tilt  it  increased.  men,  various 
types  of  time-varying  acceleration  patterns  with  zero  average  value 
were  applied  to  the  model  in  a  noosynchronous  manner  (acceleration 
cycle  was  independent  of  the  heart  cycle).  For  sinusoidal  acceleration 
patterns  the  model's  aortic  pressure  and  flow  responses  exhibited 
resonance  in  the  2  -  6  Hz  frequency  range,  as  has  been  observed  experi¬ 
mentally;  with  3  G  peak  acceleration  in  this  frequency  range,  peak 
left  ventricular  ejection  flow  rate  increased  75%,  aortic  puise  pres¬ 
sure  increased  330*4,  and  stroke  volume,  on  a  per-cycle  basis, 

Increased  over  40%,  However,  cardiac  output  was  not  changed  under  any 
applied  acceleration  pattern  with  zero  average  value,  although  it  was 
changed  with  asymmetrical  time-varying  patterns  (time-varying  patterns 
of  zero  average  value  combined  with  a  tilt).  Time-varying  acceleration 
patterns  were  also  applied  in  a  synchronous  manner.  It  was  found  that 
model  responses  effected  by  nonsynehronous  accelerations  producing 
long  beat  patterns  could  be  reproduced  synchronously. 
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CHAPTER  I 


INTRODUCTION 

It  has  become  evident  chat  acceleratory  stresses  Imposed  on  man 

in  many  of  his  modern  environments  can  have  substantial  effects  on 

cardiovascular  function.  Investigators  have  been  exploring  these 

effects  for  acme  time.  The  first  studies  in  cardiovascular  response  to 

whole-body  accelerations  dealt  primarily  with  arterial  responses  to 

constant  acceleration  patterns,  such  as  chose  produced  by  centrifuges. 

Later,  arterial  and  venous  responses  were  examined  for  prolonged 

conntant  acceleration  patterns,  as  well  as  for  constant  patterns  of 

shorter  duration.  Inue.  the  importance  of  certain  time  parameters  of 

whole-body  acceleration,  such  as  rate  of  onset  of  acceleration,  duration 

of  acceleration,  etc . .  became  apparent,  and  time-varying  accelerations 

of  various  forms  were  imposed  upon  the  cardiovascular  system*0.  More 

recently,  the  effects  on  the  cardiovascular  system  of  time-varying 
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accelerations  of  periodic  nature  have  been  investigated  ’  .  These 

studies  have  shown  that  certain  patterns  of  whole-body  acceleration  can 
produce  profound  changes  in  the  circulation.  For  example,  ; ome  patterns 
present  stresses  to  the  cardiovascular  system  which  can  drastically 
reduce  etrokc  volume,  whereas  other  patterns  may  be  used  to  assist  a 
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failing  heart  by  increasing  cardiac  output.  Primary  emphasis  in  these 


studies  has  been  given  to  arterial  responses.  However,  the  earlier 
constant-acceleration  studies  suggest  that  venous  responses  are  also 
an  important  part  of  the  total  cardiovascular  response  to  whole-body 
accelerations. 

The  purpose  of  this  study  was  to  analyze  human  cardiovascular 
response  to  whole-body  2-axis  accelerations  by  means  of  a  closed-loop 
analog  computer  simulation.  The  responses  to  both  constant  and  time- 
varying  periodic  acceleration  patterns  were  studied. 


CHAPTER  IT 


DEVELOPMENT  OF  THE  MODEL 

The  systemic  circulation  involves  the  movement  of  blood  through 
the  left  heart  and  high-pressure  arterial  system  to  the  microcirculation 
and  the  return  of  the  blood  from  the  microcirculation  into  the  right 
heart  through  the  low-pressure  venous  system.  The  pulmonary  circulation 
returns  the  blood  from  the  right  heart  to  the  left  heart.  The  arterial 
system  has  been  studied  extensively,  and  many  arterial  models,  some 
modified  to  account  for  whole-body  accelerations,  have  been 
developed^’ Although  some  venous  studies  have  been  mad<^7,  under¬ 
standing  of  the  role  of  the  venous  system  in  cardiovascular  response  tD 
whole-body  accelerations  is  not  nearly  as  advanced. 

The  divisions  into  which  the  circulatory  system  has  been  separ¬ 
ated  in  this  study  and  their  abbreviations,  also  used  as  subscripts, 
are  listed  below  (sec  Fig.  4): 

A.  Heart  and  Pulmonary  Section 

1.  Right  Ventricle  (RV) 

2.  Lungs  (LG) 

3.  Left  Ventricle  (LV) 

P.  Arterial  Sections 

1 ,  Aort  ic  Arch  (AA) 

2.  Upper  Arterial  System  (UAS) 
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3.  Lower  Arterial  System  (LAS) 

4.  Upper  Capillary  Bed  (UCB) 

5.  Lower  Capillary  Bed  (LCB) 

C.  Venous  Sections 

1.  Upper  Venous  System  (UVS) 

2.  Lower  Venous  System  (LVS) 

Right  Atrium-Venae  Cavae 
Junction  (RA' ) 

Heart  and  Pulmonary  Section 

The  equations  describing  flow  through  the  heart  and  lung  sec¬ 
tions  are  similar  to  those  used  in  an  earlier  open-loop  model6.  Con¬ 
tractions  of  the  ventricular  sections  were  achieved  by  means  of  time- 
varying  inverse  compliances  consisting  of  biased  rectified  sinusoids 
operating  at  1.25  Hz  (75  beats/minute).  The  pulmonary  section  was 
modeled  as  a  linearly  compliant  element;  its  resistance  was  included 
with  the  ventricular  valve  resistances. 

Arterial  Sections 

The  equations  describing  flow  through  the  upper  and  lowi  r 
arterial  sections  are  similar  to  those  used  in  an  earlier  model6.  The 
equations  for  flow  through  .ie  aortic  arch  section  have  the  same  form 
as  RLC  L-8ection  equations  with  the  compliance  parameter  divided 
between  the  ends  of  the  section. 
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The  equivalent  heights^  determined  by  Camill6  for  the  upper 
and  lower  arterial  systems  of  standing  man  were  modified  by  the  inclu¬ 
sion  of  the  ascending  aortic  arch  section.  This  section's  equivalent 
height  of  7.6  cm  reduced  the  eauivalent  height  of  the  upper  arterial 
system  to  29.7  cm  and  increased  the  equivalent  height  of  the  lower 
arterial  section  to  50.3  cm. 

The  capillary  beds  were  represented  as  resistive  elements  in  the 

model . 


Venous  Sections 

There  are  several  venous  system  characteristics  and  external 
mechanisms  which  affect  venous  system  function.  One  outstanding  feature 
is  the  nonlinear  pressure-volume  characteristic,  which  enables  the 
venous  system  to  contain  a  large  volume  of  blood  at  a  relatively  low 
pressure  without  being  excessively  compliant.  These  large  volumes 
become  especially  important  when  whole-body  accelerations  are  imposed 
since  accelerations  tend  to  produce  volume  shifts  in  the  circulation*®. 

An  empirical  pressure-area  relat lonship  for  a  segment  of  vein 
is  shown  in  Fig.  1A*^  .  Assuming  volume  is  proporcional  to  area,  the 


"height"  is  meant  the  difference  in  elevation  of  the  .wo  ends  of 
a  length  of  artery  in  an  applied  gravitational  field.  The  "equivalent 
height"  of  a  composite  section  representing  brand,  ing  flow-  raised  to 
different  heights  has  hcen  shown  to  be  the  average  height  if  the 
branches  weighted  as  their  flow  rates^'*®. 
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slops  of  tills  curve  represents  the  compliance  of  a  venous  segment.  In 


the  model  an  approximation  to  this  compliance  curve  was  used  to  represent 
the  composite  compliance  of  the  venules  and  small  veins  in  the  uppsr 
venous  system  as  shown  in  Fig.  IB  and  described  >y  the  following 

equation: 


puvs 


cuvs 


(11-1) 


where  VQ  represents  the  total  unstressed  volume  of  the  system,  or  the 
volume  at  which  the  compliance  changes  value  (see  Appendix  A).  The 
composite  compliance  of  the  lower  venous  system  wac  represented  likewise. 

A  second  very  important  characteristic  of  the  venous  system  is 
venous  valves.  When  accelerations  which  tend  to  reverse  flow  are 
applied,  the  valves  close  and  substantially  impede  backflow,  thus 
partially  preventing  venous  pooling  in  the  extremities.  Physically, 
these  valves  are  located  in  some  but  rot  all  of  the  venules,  small 
veins,  and  large  veins.  These  valves  were  modeled  by  incorporating 
ideal  valve  models  (ideal  diodes)  in  the  sections  representing 
equivalent  venous  flow  (flow  from  the  venules  and  small  veins  through 
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the  large  veins  to  the  region  of  the  right  atrium)  .  The  ideal  valve 
models  modify  the  flow  equations,  producing  unidirectional  flow  towards 
the  right  heart: 
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where  Q^,g  represents  bidirectional  flow,  which  could  occur  during 
accelerations  if  no  valves  were  present,  and  represents  modified 

or  physiological  flow.  The  analog  computer  representation  of  this 
plecewlse-llnear  equation  is  shown  in  Fig.  2,  The  position  of  the 
ideal  valve  in  the  circuit  is  at  the  end  of  the  equivalent  venous  flow 
section  distal  to  the  heart.  This  positioning  corresponds  to  the  obser¬ 
vation  that  no  anatomical  valves  exist  in  the  venae  cavae  near  the  heart. 

The  third  mechanism  which  affects  venous  function  is  the  ability 

of  some  of  the  large  veins  to  collapse  under  low  transmural  pressures. 

Physiologically,  this  collapse  phenomenon  manifests  itself  as  a  partial 

collapse3’4  occurring  locally  rather  than  throughout  an  entire  vessel. 

Examples  of  this  phenomenon  include  jugular  collapse  and  partial 

superior  vena  cava  collapse  upon  standing  and  collapse  of  abdominal  and 
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leg  veins  in  the  head-down  position  .  Therefore,  thi-  mechanism  tends 
to  prevent  or  delay  large  volume  shifts  due  to  acceleration  of  the  blood 
in  the  large  veins.  In  the  model,  changing  the  resistance  of  a  venous 
flow  section  without  changing  its  inertance,  since  both  vary  with  volume 
or  cross-sectional  area,  is  analogous  to  local  partial  collapse.  The 
resulting  flow  equation  is  of  the  same  form  as  Eq.  (II-2)  with  Ryyg 
nonlinear  as  shown  in  Figs.  2  and  3;  this  step  increase  in  resistance 
is  a  piecewise  approximation  to  the  continuous  collapse  depicted  in 
earlier  models1^’17.  Flow  in  the  lower  venous  system  was  modified 
slmilarl'-.  (In  Fig.  3  is  the  resistance  for  the  case  of  no  collapse, 
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or  "normal"  flow  resistance,  and  is  the  resistance  during  local 
partial  collapse,  which  is  produced  by  volume  shifts  due  to  accelera¬ 
tion.  The  volume  at  which  local  partial  collapse  is  assumed  to  occur 
is  Vc  and  is  also  show  in  Fig.  IB.) 

Although  the  external  mechanism  of  muscle  contractions  and  the 
presence  of  negative  intrachorac ic  pressure  and  its  variations  produced 

by  breathing  affect  venous  function,  these  were  not  included  in  the 
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model  for  the  sake  of  simplicity.'  '  However,  they  can  he  added  if 
desired. 


Since  the  heart,  lungs,  and  capillary  beds  were  assumed  to  have 
an  equivalent  height  of  zero,  the  equivalent  heights  determined  by  an 
earlier  investigator6  Cor  the  arterial  system  were  used  as  the  equi¬ 
valent  heights  of  the  venous  system;  therefore,  hyyS  =  37.3  cm  and 
W.VS  “  cm,  where  h  indicates  equivalent  height.  Accelerators- 

stresses  modified  flow  equations  in  those  sections  of  the  model  which 
had  non-zero  equivalent  height  (see  Appendix  Bl. 

The  electrical  equivalent  circuit  of  the  composite  model  is 
shown  in  Fig.  4,  and  the  corresponding  parameter  values  are  given  In 
Table  1.  The  polarities  of  the  constant  voltage  sources  in  Fig.  4 
represent  a  head-upward  tilt;  for  a  head-downward  tilt  t he  polarities 
arc  reversed.  Time-varying  accelerations  with  zero  average  value  are 
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3ure  4  Composite  Electrical  Equivalent  Diagram 


TABLE  1:  Parameters  of  the  M o del 


RESISTANCE 

INERTANCE 

COMPLIANCE 

dyne-sec 

dyne-sec2 

cm5  „ 

OTHER 

SECTION  cttiS 

cm^ 

.  X  10 

dyne 

PARAMETERS 

RV 

RtV  «  20; 

Rpv  -  10 

LG 

=  60,000 

LV 

*MV“  l°l 
rav  c  20 

AA 

raa  *  48-2 

Ha  -  2 

CLAA  =  600  •> 

cuaa  *  600 

UAS 

Has  "  241 

Has  “  33 

CUAS  *316 

LAS 

rlas  •  52-5 

llas  “  5-33 

CUAS  =  552 

UCb 

Hcb  •  4'000 

LOB 

rlcb  -  l**°° 

UVS 

Ru  *  510 

uvs 

Rr  *  4 ,100 

UVS 

Hvs  3  13 

cuvs  =  3>560^ 

V0  *  474  cc 

uvs 

vr  «  382  cc 

Hvs 

LVS 

\v,  ■  7,0  L^’35 

Rp  •  35,000 

''LVS 

[l] 

CLVS  =  13,000 

V_  *=  2,060  cc 

°LVS 

V_  “  2,000  cc 

Hvs 

ra' 

CRA'C  9,360 

1  1 

Vn  luc 

of  ( i ni  U-  cnmpl 

i  ;mcc 
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represented  by  the  sources  phVw  (see  Appendix  B).  The  analog  computer 
diagram  of  the  composite  model  is  shown  in  Figs.  5  and  6;  the  digital 
voltmeter  readings  of  the  coefficient  potentiometers  and  the  parameters 
represented  by  the  corresponding  potentiometer  setting*ampl if i?r  gain 
products  are  given  in  Table  2.  All  pressures  were  scaled  in  mm  Hg;  all 
other  variables  were  scaled  in  cgs  units.  The  heart  valves  and  their 
abbreviations  are  tricuspid  valve  (TV),  mitral  valve  (MV),  pulmonary 
valve  (PV),  and  aortic  valve  (AV) . 
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TABLE  2.  Analog  Potentiometer  Settings  ant 
Their  Corresponding  Parameter.. 


Pot  No. 


Pot  Setting 


Parameter  Represented 


l/Rpv 


1/cra' 


I/Clg 


1/cuaalaa 

1/cuaaluas 

1,cv,\al\as 

Acceleration^ 1  - 

^cra'luvs 

1/Icra'llvs 

Acceleration 

Acceleration 


'  AV 

Acceleration 

1/CLAARAV 

1/claaLaa 

RAA^'AA 
Accel  oration 


l"tu  occclcrat Inn  ..oUntl-ctce  HttlW  »>■<«"  1 

or  peak  accelerat Iona . 


V 


'V 


TABLE  2  Continued 


Pot  No. 

Pot  Setting 

24 

.015 

31 

.731 

32 

.950 

33 

.101 

35 

.084 

36 

.843 

37 

.256 

36 

.  392 

39 

.199 

40 

.190 

41 

.100 

42 

.315 

44 

.191 

46 

.985 

47 

.544 

48 

.625 

50 

.238 

51 

.184 

52 

.951 

53 

.203 

54 

.028 

Parameter  Represented 
Blood  Infusion^  - 

ruas/luas 


uas 


r  2 1 


l/C 

^LllAS 

1/Rucr> 
1^cuvs' 

1 /Lyyg 

rnuvs/luvs 

V°UVS 

V°LVS 

VC 

LVS 

RCUVS/LUVS 


'CUVS 
lLA 
l/C 


rlas/lla^ 


1/L 

1/R- 

l/C 

1/L 


US 

US 


LCB 

LVS 

LVC 

/I.x 


[21 


R«W  “LVb 

l.C. 


■potentiometer  24  was  used  to  produce  a  constant  signal  which  uas 
presented  to  a  volume  integrator  for  a  sufficient  period  of  time  to 
build  up  blood  volume  in  the  system.  Potentiometer  54  was  used  to  set 
volume  Integrator  Initial  conditions  so  tbit  overloads  would  not  occur 
before  the  system  reached  a  steady  state.  Together,  these  potentio¬ 
meters  established  a  total  system  volume  of  4,800  cc . 

|2l 

'Thin  potentiometer  represents  the  finite1  compl 1  ante  of  I  .10  corre¬ 
spond  Ing  venmi.i  system. 
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CHAPTER  III 


RESULTS  OF  SIMULATING  CIRCULATORY  RESPONSES  TO  ACCELERATIONS 

Constant  Accelerations 

The  model's  response  to  zero  acceleration,  simulating  man  in 
the  prone  or  supine  position  (all  the  blood  vessels  are  assumed  to  lie 
in  the  horizontal  plane),  is  shown  in  Fig.  7.  Only  those  variables  of 
the  model  which  have  a  physiological  counterpart  are  shown.  On  the 
arterial  aide  these  are  left  ventricular  pressure  (PLV),  ejection  flow 
out  of  the  left  ventricle  (.Qq^)  ,  pressure  in  the  lower  aortic  arch 

r  n 

(Pla^).  volume  of  the  left  ventricle  (V^,)  ',  flow  through  the  aortic 

arch  (Qy^).  and  pressure  in  the  upper  aortic  arch  (P^AA) •  The  simulated 

pressure  Plm  ^rresponds  to  aortic  pressure  Just  outside  the  left 

ventricle:  P,,,,  corresponds  to  aortic  pressure  at  the  junction  of  t Me 
UAA 

aorta  and  the  carotid  and  subclavian  arteries.  On  the  venous  side  the 
variables  shown  are  ptessure  outside  the  right  ventricle  or  composite 
pressure  of  the  right  atrium  and  large  veins  (P^),  flow  into  the 
right  ventricle  < >  and  pressure  In  the  right  ventricle  (P^y). 

Cl]ln  the  model  V^y  -  (QILV  '  ^0LV)dt  =  ’  °!LVdt  *  ■  <Uvdt  ■  where  °ILV 
i 9  flow  into  the  left  ventricle  and  the  last  integral  represents  stro' e 
volume  (SV);  since  QI(V  and  Q^y  always  occurred  in  different  portions 
of  the  heart  cycle,  tne  value  of  SV  Is  the  absolute  difference  between 
end-d lastol 1c  V^y  and  end-systolic  v^y. 
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Figure  7.  Response  of  Model  to  Zero  Acceleration: 

Typical  Simulated  Pressures,  Flows,  and 
Volume 
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Figure  7.  Continued 


Some  of  Che  phys iological  counterparts!  of  these  variables  are  shown  in 

Fig.  8. 

Model  responses  to  2ero  acceleration  which  are  not  shown  are 
pressures  and  flows  in  the  upper  and  lower  arterial  and  venous  systems 
and  in  the  lung  section.  The  upper  and  lower  arterial  systems' 
pressures  and  flows  wete  similar  in  detail  to  those  generated  by  an 
earlier,  open-loop  model^. 

The  volume  distribution  throughout  the  model  with  zero  accelera¬ 
tion  is  shown  in  the  following  table: 


TABLE  3.  Volume  Distribution  of  the  Model 


Venous  Sections 

Thoracic  Sections 

Arterial 

Sec  t ions 

UVS  580  cc 

RV 

50  cc 

LAA 

81  cc 

RA'  120  cc 

LG 

620  cc 

UAA 

8 1  cc 

LVS  3,000  cc 

LV 

144  c  c 

UAS 

x> 

o 

o 

n 

LAS 

7  l  cc 

The  total  system  volume  was  approximately  4,300  cc. 

The  model's  response  to  positive  1  G  acceleration,  simulating  a 
head-up  tilt,  is  shown  in  Fig.  9.  Response  of  the  model  to  negative 
1  G  acceleration,  simulating  a  head-down  tilt,  is  shown  in  Fig.  10. 

For  a  simulated  head-up  tilt,  lung,  ventricular,  and  upper 
venous  system  volumes  decreased, and  lower  venous  system  volume  increased 
The  volume  changes  were  quant  itativcly  similar  to  chose  observed  experi¬ 
mentally  in  man^.  These  effects  were  reversed  in  a  simulated  head-down 
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Acceleration  Left  Ventricular  Left  Ventricular  Lower  Aortic  Left  Ventricular 

(G‘s)  Volume  (cc)  Pressure  Arch  Pressure  Ejection  Flow 

(mm  Hg)  (mm  Hg  )  (cc/sec) 


Lower  Aortic  Arch  150 
Pressure  (mm  Hg)  gQ 

30 

0 

12 

Right  Atrial  9 

Pressure  (mm  Hg  ) 

3 

0 


Right  Ventricular 
Inflow  (cc/sec) 


350 

250 


50 

0 


S  DIVISION  r.Cr 

°f<iN  T  i  ()  i*.  ■ 


°n 


Acceleration  (G's) 


tilt.  In  the  head-up  tilt  collapse  of  the  upper  venous  segment  between 
UVS  and  RA'  did  occur  and  QyVS  was  altered  accordingly;  however,  in  the 
head-down  tilt  collapse  of  the  lower  venous  segment  between  LVS  and  RV 
did  not  occur. 


Wonsynchronoua  Time-varying  Accelerations 

Time-varying  accelerations  of  periodic  nature  were  applied  to 
the  modal  in  a  nonsynchronous  manner.  Nonsynchronous  accelerations  are 
Chose  which  are  applied  in  a  continuous  mode--the  acceleration  pattern 
is  not  related  to  tne  heart  cycle.  This  situation  of  two  independent 
oscillators  was  characterized  by  the  production  of  beat  patterns^1-  in 
Che  flows  and  pressures  in  the  model. 

A  typical  response  of  the  model  to  nonsynchronous  time-varying 
acceleration  producing  long  beat  patterns  is  shown  In  Fig.  HA,  where 
the  acceleration  pattern  is  sinusoidal  with  a  peak  amplitude  of  2  G 
and  a  frequency  of  nearly  1.25  Hz.  Figure  1 IB  shows  a  selected  portion 
of  FLg.  11A  recorded  at  a  faster  chart  speed. 

A  typical  response  of  the  model  to  nonsynchronous  time-varying 
acceleration  producing  short  beat  patterns  is  shown  in  Fig.  12a,  where 

^'The  period  of  a  beat  pattern  is  the  length  of  time  required  for  a 
nonsynchronous  periodic  acceleration  pattern  to  "drift  through"  the 
heart  cycle  so  that  its  occurrence,  or  timing,  in  a  particular  heart 
cycle  is  similar  to  that  in  an  earlier  cycle.  This  timing  can  be 
indicated  by  defining  a  parameter  tD  as  the  t imu  delay  in  milliseconds 
between  the  positive  zero  crossing  of  acceleration  and  the  initial  rise 
in  the  inverse  compliance  curve  or  pressure  of  the  left  ventricle. 


Figure  IIA.  Typical  Response  of  Mode!  to  Simulated  Non - 
synchronous  2G,l.25Hz  Sinusoidal  Acceleration 
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Figure  MB.  Selected  Portio? 
Figure  M  A . 
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the  acceleration  pattern  is  sinusoidal  with  peak  amplitude  of  3  G  and 
frequency  of  3  Hz.  An  experimental  response  to  this  same  acceleration 
pattern  la  presented  in  Fig.  12B9, 

Response  of  the  model  to  nonsynchronous  sinusoidal  accelerations 
of  different  amplitudes  and  frequencies  is  shewn  In  Figs.  13A,  13B,  and 
13C.  ihe  experimental  variables  are  defined  as  follows: 


(1) 

%AX: 

(2) 

<^MAX: 

(3) 

^MIN5 

(4) 

SW 

(3) 

SVM1N: 

maximum  difference  of  pressure  In  the  lower 
aortic  arch  (Pj^)  during  any  single-  excursion**  ’  ^ 
(see  Fig.  12B) . 

maximum  peak  left  ventricular  ejection  flow  (Q__  ,,) 

r  Q  vLV 

during  a  single  heart  cycleD,y  (see  Fig.  I2fl) . 

minimum  peak  left  ventricular  ejection  flow  (On,,,) 

AO  ULV 

during  a  single  heart  cycle0^  (see  Fig.  12B) . 

maximum  stroke  volume  of  a  single  heart  cycle 
(see  Fig.  12A) . 

minimum  stroke  volume  of  a  single  heart  cycle 
(see  rig.  12a). 


Response  of  the  model  to  other  nonsynchronous  time-varying 
accelerations  was  also  determined.  For  the  case  of  square-wave 

acceleration  patterns,  ^iP^,  QmaX>  fylN’  SVMAX>  Snd  SVMIN  curves  were 
obtained  which  were  similar  to  the  corresponding  curves  presented  in 
Figs.  13A,  13B,  and  13C  but  which  had  generally  greater  deviation  from 
control  than  the  sinusoidal  responses  for  any  given  G  level. 

A  nonsynchronous  asymmetric  rectangular  waveform  of  l  G  peak 


magnitude  and  a  frequency  of  approximately  1.25  Hz  was  applied  to  the 


Lower  Aortic  Arch 
Pressure  (mm  Hg) 


Figure  12  A  Typicoi  Respot .*-•  of  Model  to  Simulated  Nonsynch>onous  3G, 
3  HZ  Sinusoidal  Acceleration 
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Figure  I2B.  Examples  of  Recordings  of  A r t e » ia I 

Pressure  and  Aortic  F  iow  from  a  Dog, 
at  Rest  and  under  3(j  Vibration  oi  o  nz  . 
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Acceleration  Frequency / Hear t  Frequency 

Figure  13  B.  Simulated  QMAX  and  Q MIN  vs.  Normalized  Acceleration  Frequency  with 
Acceleration  Amplitude  as  a  Parameter. 


mcJel'-  Shown  In  Fig.  Id  is  the  model's  response  to  this  waveform. 


A  selected  portion  of  the  slow  chart  speed  is  shown  at  a  faster  chart 
speed.  Note  that  SV  is  above  the  normal  value  of  68  cc  for  all 
acceleration  timings.  For  a  2  G  rectangular  waveform  with  positive  and 
negative  portions  of  the  same  durations  as  for  the  1  G  case  above,  the 
model's  responses  were  similar  to,  but  greater  in  magnitude  than,  the 
1  G  responses. 

Synchronous  Time-varying  Accelerations 

A  time-varying  acceleration  pattern  is  considered  synchronous 
if  its  frequency  is  a  fixed  integral  multiple  of  the  heart  frequency. 

In  this  idealized  situation  no  beat  patterns  result;  l.e,  ,  each  cycle 
of  the  model's  response  Is  identical  to  the  previous  cycle,  once  the 
model  reaches  steady-state  (usually  after  only  a  few  cycles). 

Several  patterns  of  synchronous  time-varying  acceleration  were 
applied  to  the  model.  Shown  in  Fig.  15  is  the  response  of  the  model 
to  synchronous  sinusoidal  acceleration  of  1.25  Hz,  2  G  peak  amplitude, 
and  tD  *•  0  milliseconds. 

Presented  in  Fig.  16  is  the  simulated  stroke  volume  when 
synchronous  sinusoidal  accelerations  of  1.25  Hz  and  l  G  and  2  G  peak 

^In  the  model  an  acceleration  pattern  with  a  non-zero  average  value 
is  the  same  as  a  pattern  with  zero  average  value  combined  with  a  tilt. 
The  average  value  of  this  1  G  rectangular  waveform  is  about  -1/3  G, 
corresponding  to  an  approximate  19,r  head-down  tilt. 
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Volume  Pressure  Arch  Ventricular 
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Figure  15.  Response  of  Model  to  Synchr  onous  2G  , !.  2 5  Hz 
Sinusoidol  Acceleration  with  tp  =  0  Milliseconds 


amplitude  were  applied  to  the  model  at  various  values  of  tp.  The 
simulated  stroke  volume  obtained  with  square-wave  accelerations  of 
1.25  Hz  and  1  G  and  2  G  amplitude,  in  which  tD  was  again  the  parameter, 
was  similar  to  that  shown  for  sinusoidal  accelerations  in  Fig.  16  in 
the  sense  that  stroke  volume  was  never  substantially  different  from 
that  of  zero  acceleration. 

In  Fig.  17  the  model's  response  to  a  synchronous  asymmetric 
rectangular  acceleration  pattern  of  1  C  peak  amplitude  and  tp  =  760 
milliseconds  is  shown.  This  rectangular  waveform  has  the  same  positive 
and  negative  durations  as  the  nonsynchronous  rectangular  acceleration 
pattern  described  earlier  (p.28). 
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Figure  17.  Response  of  Model  to  o  Simulated  Synchronous 
Asymmetric  Acceleration  Pottern  with  tp*760 
Milliseconds. 
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CHA1TER  IV 


DISCUSSION  AND  CONCLUSIONS 


Conatant  Accelerations 

In  the  prone  position  the  arterial  waveforms  of  the  model  are 
similar  to  those  observed  physiologically.  However,  the  model's  peak 
left  ventricular  pressure  and  mean  left  ventricular  volume  are  higher 
than  tho<*e  observed  in  man,  whereas  stroke  volume  and  pulse  pressure 
are  lower.  Although  the  discrepancies  could  have  been  reduced  by 
reducing  both  the  frequency  and  magnitude  of  the  ventricular  compliance, 
a  normal  heart  rate  was  first  chosen  and  a  magnitude  of  the  compliance 
was  then  selected  which  compromised  stroke  volume  and  left  ventricular 
pressure  accuracies. 

The  venous  waveforms  of  the  model  do  not  match  physiological 
data  aa  accurately  as  the  arterial  waveforms  since  such  factors  as 
negative  Intrathoraclc  pressure,  breathing,  and  right  atrial  contrac¬ 
tions  were  not  taken  Into  account  in  the  model.  Thus,  right  atrial 
Junction  pressure  seems  rather  high,  but  this  pressure  is  a  composite 
pressure  due  to  the  extensive  lumping  of  the  model  and  its  exact 
physiological  identification  is  not  possible.  However,  if  is  somewhat 
similar  in  waveform  to  right  atrial  pressure  in  man.  The  right 
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ventricular  pressure  waveform  can  be  made  mure  like  its  physiological 


counterpart  by  decreasing  the  resting  ventricular  compliance,  but  this 
change  results  in  high  Initial  flow  Into  and  out  of  the  right  ventricle. 
Therefore,  more  steady  flows  were  achieved  at  the  expense  of  the  right 
ventricular  pressure  simulation, 

Because  the  values  of  arterial  compliance  needed  to  reproduce 
arterial  pressure  waveforms  accurately  in  the  model  were  much  smaller 
than  physiological  values,  volumes  were  very  low  in  the  upper  and  lower 
arterial  sections.  With  linear  compliances  it  is  impossible  to  achieve 
in  a  highly  lumped  arterial  model  the  degree  of  simultaneous  accuracy 
of  pressure  and  volume  possible  In  a  more  distributed  model.  If 
arterial  volumes  were  to  be  raised  by  Increasing  arterial  compliance 
In  the  lumped  model,  the  arterial  pressure  waveforms  would  lose  much 
of  their  characteristic  variation.  On  the  other  hand,  If  volumes  were 
to  be  raised  by  Increasing  mean  arterial  pressures,  serious  Inaccuracies 
In  pressure  and  flow  waveforms  throughout  the  model  would  result. 

Since  an  error  In  arterial  volume  is  a  small  percentage  of  the  total 
sybtem  volume,  arterial  volume  accuracy  was  sccrificcd  for  pressure 
accuracy. 

The  above  situation  was  circumvented  on  the  venous  side  of  the 
model  by  using  plecewise-1 inear  compliances,  which  provide  large 
residual  volumes  at  zero  pressures.  The  parallel  combination  of  many 
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venous  segments  which  have  the  pressure-volume  characteristic  depicted 
in  Fig.  1A  (p.  7 )  would  ideally  contain  a  rather  large  volume  at  zero 
pressure  but  would  begin  to  resist  volume  increases  once  all  the 
segments  have  been  filled  to  the  point  of  distention.  This  observation, 
together  with  the  fact  that  large  venous  volumes  do  exist  at  relatively 
low  pressures  although  the  venous  system  is  not  excessively  compliant, 
provided  the  basis  for  postulating  an  overall  p 'es sure  volume  relation¬ 
ship  such  as  that  of  Fig.  IB.  In  agreement  with  physiological  observa¬ 
tion,  the  upper  and  lower  venous  systems  of  the  model  contained  most  of 
the  venous  volume  .  The  thoracic  volume  of  the  model  was  low  because 
of  the  simplicity  of  the  lung  model.  However,  the  venous  and  thoracic 
volumes  were  similar  to  the  general  volume  distribution  obtained  in  a 
more  extensive  model17. 

Since,  in  a  head-up  tilt,  nervous  regulation  Is  involved  in 
adjusting  the  compliance  of  the  physiological  system  (venomotor  tone), 
the  slopes  of  the  pressure-volume  relationships  for  the  venous  sections 
of  the  model,  which  have  no  neural  control  simulations,  were  chosen  to 
result  in  volume  shifts  of  proper  magnitude^  3nd  therefore  to  repre¬ 
sent  compliance  values  averaged  over  all  conditions  of  venomotor  tone. 

When  a  man  stands  up,  aortic  and  left  ventricular  pressures 
and  stroke  volume  tend  to  decrease;  however,  nervous  control  changes 
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in  peripheral  resistance  and  venomotor  tone  occur,  tending  to  raise 
these  quantities'’11,1^1.  Still,  cardiac  output  will  be  down  by  as 
much  as  20  percent  because  the  decrease  in  stroke  volume  more  than 
offsets  the  increase  in  heart  rate^.  In  the  model  changes  in  peri¬ 
pheral  resistances  were  not  simulated;  therefore,  in  a  simulated  head- 
up  tilt  aortic  and  left  ventricular  pressures  and  stroke  volume 
decreased  to  about  60  to  70  percent  of  their  original  value.  These 
responses  are  similar  to  those  obtained  in  a  more  detailed  model 

1  g 

operating  with  only  one  nervous  control  loop  (heart-rate  control)  , 
but  the  percentage  decreases  in  the  detailed  model  were  greater  than 
those  In  the  model  used  in  this  study. 

A  majot  difference  between  simulated  and  physiological  tilt 
response  occurred  with  right  atrial  pressure  during  a  head-down  tilt. 
In  the  model  this  pressure  increased,  whereas  it  has  been  consistently 
found  to  decrease  with  a  head-down  tilt  in  man10,11.  This  apparent 
limitation  of  the  model  was  caused  in  part  by  the  extensive  lumping 
of  the  upper  venous  system  compliance,  much  of  which  was  placed  at 
heart  level  (C^y )  In  the  model  As  a  result,  the  shift  in  olood 
volume  away  from  heart  level  and  the  associated  shift  in  hydrostatic 
Indifference  point1*  induced  by  a  simulated  head-down  tilt  were  not 
as  large  in  tho  model  as  those  that  would  be  observed  experimentally. 
However,  till  limitation  was  an  advantage  in  other  ways.  Both  the 
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model  and  anesthetized  animals*  showed  increases  in  stroke  volume  and 
lower  venous  system  flow  under  conditions  of  head-down  tilta.  In  the 
model,  increased  stroke  volume  can  be  related  directly  to  increased 
,  which  should  cause  in  sequence  increases  in  right  ventricular, 
pulmonary,  and  left  ventricula:  pressures  and  volumes,  Theoe  increased 
pressures  and  volumes  should  bring  about  increased  stroke  volumes  sir.* 
the  inverse  compliance  curve  of  the  left  ventricular  model  is  invariant 
under  postural  changes.  In  the  animals,  neural ly-cor.trol led  changes 
in  peripheral  resistance  and  venomotor  tone,  heart  rate,  and  ventricular 
compliance  probably  occurred,  more  than  offsetting  the  decrease  in 
right  atrial  pressure  u  ,d  resulting  in  a  net  increase  in  stroke  volume. 

Excessively  large  volumo  shifts  related  to  large  acceleratory- 
induced  pressure  changes  on  the  venous  side  of  the  model  were  avoided 
by  using  piecewlse-linear  resistances  to  reduce  flowB  under  those 
conditions.  These  discontinuities  simulated  the  increase  in  resistance 
to  blood  flow  that  occurs  physiologically  when  local  partial  collapse  of 
large  veins  occurs  as,  for  example,  during  a  head-up  tilt. 

Elucidation  of  the  model's  responses  to  tilt  experiments  is 
provided  by  examination  of  the  heart-lung  section  of  the  model.  It 
appears  that  change  in  mean  flow  rate  through  the  simple  herrt-lung 
section  is  determined  by  the  change  in  mean  right  atrial  pressure 
(P  ,)  together  with  the  change  in  mean  pressure  outside  the  left 
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ventricle  (P^).  For  example,  if  P^,  is  held  constant  while  P^ 
is  decreased,  flow  rate  through  the  heart-lung  section  will  increase 
(as  has  been  observed  in  an  open-loop  model^  with  head-up  tilt);  or, 
if  Pj^'  is  held  constant  and  P]_^  increased,  flow  should  decrease. 
Conversely,  when  is  held  constant  and  PR4/  is  allowed  to  change, 

the  flow  should  decrease  if  P^'  is  decreased  and  increase  if  P  <  is 
increased.  Now,  consider  postural  changes  in  a  closed-loop  model 
where  both  and  change.  In  the  head-up  tilt,  where  Pg^/ 

decreased  367.  and  decreased  257.,  stroke  volume  decreased  387..  In 

the  head-down  tilt,  where  p  ,  increased  347.  and  Pjj^  increased  287., 
stroke  volume  Increased  327..  In  all  the  preceding  cases  it  is  seen 
chat  acroke  volume  of  the  model  varies  with  right  attial  pressure;  in 
addition,  if  a  linear  dependence  of  stroke  volume  on  Pp^/  and  P^A  is 
assumed,  it  is  found  that  the  Influence  of  Pp^'  on  stroke  volume  is 
about  twice  as  great  as  the  influence  on  P  .  Thus,  in  the  model 
constant  accelerations  affect  mean  right  atrial  or  filling  pressure, 
which  in  turn  affects  9troke  volume.  The  importance  of  the  effect  of 
postural  changes  on  right  atrial  pressure  hes  also  been  recognized  in 

7 

physiological  studies  . 

Time-varying  Accelerations 

Comparison  of  Figs.  12A  and  12B  shows  similarity  between  pres¬ 
sure  end  flow  responses  obtained  from  the  muuel  i u i  »imulaceo  sinusoidal 
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accelerations  and  responses  obtained  experimentally  with  whole-body 
sinusoidal  accelerations.  Also,  the  and  peak  Q  curves  of  Figs. 

13A  and  13B  seem  to  match  experimental  date  more  closely  than  these 
curves  produced  by  an  open- loop  model^  because  the  curves  of  Figs.  13A 
and  13B  have  a  wider  frequency  response.  Thus,  there  is  fairly  good 
qualitative  and  quantitative  agreement  between  the  model  investigated 
in  this  study  and  experimental  results.  However,  these  agreements  must 
be  viewed  in  light  of  the  fact  that  although  the  simulated  accelerations 
were  applied  directly  to  the  model,  che  experimental  accelerat lone  may 
be  modified  by  the  unknown  transfer  function  between  the  shake  table 
and  the  cardiovascular  system. 

The  response  of  the  simulated  stroke  volume  (Fig.  13C)  indicates 
the  fact  that  for  any  particular  sinusoidal  ecceleration  pattern,  atroke 
volume  was  enhanced  during  one  portion  of  the  Leat  pattern  and  diminished 
for  the  other  portion  of  the  beat  pattern  in  such  a  way  that  cardiac 
output  waa  about  the  same  as  for  zero  acceleration.  In  addition,  at 
multiples ,of  the  heart  frequency,  atroke  volume  was  near  control  through¬ 
out  the  length  of  the  beat  pattern. 

Although  the  square-wave  acceleration  pattern  effected  great 
changee  in  pressures  and  large  peak  flow  rates  in  the  model,  it  did  not 
Increase  cardiac  output  at  cither  1,  2,  or  3  G  peak  amplitudes.  As  In 
the  case  for  sinusoidal  acceleration,  simulated  stroke  volu.se  was  first 
enhanced  and  then  decreased  within  a  beat  pattern  in  such  a  way  that 
cardiac  output  did  not  increase. 
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The  rectangular  waveform  in  Fig.  14  provides  an  excellent 


example  of  simulated  stroke  volume  when  the  acceleration  pattern  was 
asynmetrlcal .  The  figure  shows  that  for  tinie-varying  acceleration 
patterns  cardiac  output  could  be  changed  using  an  asynmetrlcal  pattern. 
It  is  also  seen  that  this  pattern  resulted  in  increased  mean  right 
atrial  and  lower  aortic  arch  pressures,  which  is  similar  to  the  results 
obtained  for  constant  accelerations  in  a  head-downward  tilt.  Thus,  it 
nny  have  been  the  tilt  that  produced  the  increased  filling  pressure  and 
Increased  cardiac  output. 

Comparison  of  Fig.  1 1 B  to  Fig.  15  and  Fig.  14B  to  Fig.  17 
illustrates  a  result  that  was  consistently  found  in  this  study;  that  is, 
the  simulated  responses  of  the  model  occurring  at  a  particular  timing  of 
an  acceleratory  pattern  which  was  producing  long  beat  patterns  could  be 
reproduced  synchronously  by  using  the  same  pattern  and  timing.  Synchro¬ 
nous  accelerations  at  heart  frequency  and  multiples  of  heart  frequency 
were  used  and,  in  all  cases  examined,  reproduced  the  responses  obtained 
with  nonsynchronous  acceleration  of  corresponding  pattern  and  timing. 
These  results  indicate  that  any  transient  responses  resulting  from 
acceleratory  timing  changes  were  of  relatively  short  duration. 

Because  of  inertial  effects  the  responses  of  the  model  to 
simulated  time-varying  accelerations  were  more  difficult  to  analyte 
than  the  responses  to  simulated  constant  accelerations.  However,  some 
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observations  were  consistently  evident.  For  example,  when  the  direction 

of  acceleration  of  a  square  wave  pattern  shifted  to  headward,  the  initial 

response  in  P.  (lower  aortic  arch  pressure)  was  a  slight  rise,  indi- 

eating  a  shift  of  blood  toward  the  heart  in  the  ascending  aortic  arch 

section.  Since  this  section  had  a  smaller  inertance  and  was  closer  to 

the  left  ventricle  than  the  other  arterial  sections,  it  determined  the 

initial  response  of  P_  ,  .  The  length  of  this  section  (7.6  cm)  limited 

LAA 

the  magnitude  of  the  >-ire.  The  initial  rise  in  P,  ,  was  followed  bv  a 

LAA 

considerable  drop,  corresponding  to  a  shift  of  blood  footward  (away  from 

the  heart)  into  the  lower  arterial  system,  reducing  volume  and  therefore 

pressure  in  the  aortic  arch.  After  a  longer  interval,  P  increased  to 

b,AA 

near  original  values  as  the  blood  in  the  relatively  high- inertance  upper 
arterial  system  shifted  footward  (toward  the  heart).  This  aortic  refill¬ 
ing  may  have  been  assisted  by  backflow  from  the  lowe"  arterial  section^-1. 
These  effects  were  reversed  when  the  direction  of  acceleration  'hifted 


to  footward.  With  sinusoidal  accelerations  similar  cifects  were  present 
but  were  often  not  as  evident  because  this  acceleratory  pattern  tended  to 


smooth  the  inertial  effects  much  more  than  a  pattern  with  abrupt  changes. 


"The  oscillatory  frequency  of  the  step  response  of  the  lower  arterial 
section  was  about  twice  that  of  the  upper  section,  so  tha-,  depending 
on  the  initial  conditions  imposed  on  the  sections  by  previi.es  accelera¬ 
tion,  the  second  oscillatory  half-cycle  (flf-t  compliant  r< action)  of 
the  lower  section  may  have  occurred  during  the  first  oscil.atory  half¬ 
cycle  (primary  inertial  respon  e)  of  the  upper  section. 
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From  the  preceding  discussion  it  is  evident  that  if  left  ven¬ 
tricular  ejection  is  to  be  enhanced  by  modifying  aortic  pressure  in 
the  model,  an  acceleratory  pattern  must  be  chosen  which  optimizes  the 
inertial  effects  of  all  the  arterial  sections.  In  this  study  some 
symmetrical  acceleratory  patterns  increased  peak  left  ventricular 
ejection  flow  rate  by  nearly  1007,  for  a  portion  of  the  usually  long 
beat  pattern  (Fig.  13B).  However,  these  patterns  did  not  produce 
maximum  stroke  volumes;  rather,  maximum  stroke  volumes  occurred  at 
acceleratory  frequencies  producing  short  beat  patterns  (Fig.  13C). 

Heart  cyclss  with  large  stroke  volumes  were  preceded  or  followed  by 
cycles  with  low  stroke  volumes.  As  a  result,  over  a  beat  period 
cardiac  output  remained  about  the  same  as  for  zero  acceleration.  Thus, 
with  symmetrical  accelerations  large  peak  ejection  flow  rates  could  be 
produced  synchronously,  but  increased  cardiac  output  could  not  be  sus¬ 
tained.  With  an  asymmetrical  acceleratory  pattern  increased  cardiac 
output  was  sustained  although  peak  left  ventricular  ejection  was  not 
maximized  (Figs.  14  and  17).  Therefore,  optimization  of  the  inertial 
effects  of  the  arterial  system  may  not  result  in  Increased  cardiac 
output . 

The  data  obtained  in  this  study  do  not  indicate  that  zero-mean 
accelerations  caused  sustained  increases  in  right  atrial  pressure. 

This  observation  might  partially  explain  why  increased  cardiac  output 
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did  not  occur  with  symmetrical  acceleratory  patterns,  since  cardiac 
output  varied  with  mean  rlglt  atrial  pressure  under  applied  constant 
acceleration.  Although  time-varying  accelerations  combined  with  tllt9 
did  Increase  mean  right  atrial  pressure  and  cardiac  output,  It  Is  not 
evident  that  the  time-varying  patterns  were  responsible  for  these 
changes . 

Variation  of  Model  Parameters 

In  the  nonlinear  venous  sections  of  the  modal  three  parameter 
values  have  not  been  identified  physiologically;  these  are  V^,  V^,  and 
Rc  (Chapter  II),  The  value  of  Vq  of  a  section  was  chosen  in  conjunc¬ 
tion  with  the  value  of  finite  compliance  to  give  approximately  normal 
operating  volumes  and  pressures  in  the  prone  position.  The  value  of 
V  helped  to  determine  the  amount  of  volume  a  section  would  lose  due 
to  sustained  acceleration.  The  parameter  Rc  determined  how  the  flow 
was  curtailed  to  prevent  excessive  volume  Iosb  of  a  section.  With 
constant  accelerations  the  system  responded  in  an  overdamped  manner 
with  large  values  of  and  in  an  underdamped  manner  with  relatively 
small  values  of  R^.  Thus,  for  a  simulated  tilt,  R  could  be  adjusted 
so  that  the  steady-state  equivalent  flow  through  a  venous  section  was 
almost  constant,  as  with  zero  acceleration,  or  exhibited  flutter, 

1 2 

which  is  characteristic  of  flow  through  a  vertical  collapsible  tube 
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The  acceleratory  patterns  which  invoked Rc  and  the  infinite 

compliances  were  constant  accelerations  or  '  ^me-varying  accelerations 

of  low  frequency.  For  example,  with  sinusoidal  acceleratory  patterns 

of  1,  2,  and  3  G’s  peak  amplitude  and  frequencies  from  0.5  to  10  Hz, 

Vmrc^n  f°r  frequencies  below  1.5  Hz  and  R„  was  invoked  only 
uvs  °UVS  Cuvs  y 

by  the  1  G,  0.5  Hz  pattern.  For  the  sinusoidal  patterns  above  and  all 

other  acceleratory  patterns  in  this  study,  V,  >Vn  in  all  cases  and 

LVS  0LVS 

R^.  was  never  invoked. 

CLVS 

Variation  of  the  phase  between  the  simulated  right  and  left 
ventricular  contractions  (p.  4)  did  not  appear  to  change  the  model's 
responses  to  accelerations.  Pressure  and  flow  responses  to  simulated 
acceleratory  patterns  with  ventricular  contractions  in  phase  were 
identical  to  those  responses  obtained  with  ventricular  contractions 
180°  out  of  phase.  The  Inverse  compliances  simulating  ventricular 
contractions  used  in  this  study  are  presented  in  Fig.  18. 

Conclu3 ions 

Some  of  the  major  conclusions  drawn  from  this  study  are: 

1.  Simulated  tilts  produced  sustained  responses  in  cardiac  output:  in 
heab-up  tilts  cardiac  output  decreased,  in  head-down  tilts  it 
increased . 

2.  Simulated  zero-mean  time-varying  accelerations  did  not  change  cardiac 
output,  although  they  did  change  peak  left  ventricular  ejection  flow 
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Left  Ventricular 
Inverse  Compliance, 
l/CLv  (t)  (volts) 


Right  Ventricular 
Inverse  Compliance 
Times  4,4/CRV(t) 
(volts) 


Figure  18.  Simulations  of  Ventricular  Contractions  in 
the  Model. 


rate  by  757.  and  aortic  pulse  pressure  by  3307..  On  a  per-cycle 
basis,  stroke  volume  was  also  changed  by  407.  or  more. 

3.  A  simulated  head-down  tilt  combined  with  a  zero-mean  tine-varying 
acceleratory  pattern  did  produce  increased  cardiac  output. 

4.  The  model's  aortic  pressure  responses  to  zero-mean  sinusoidal 
accelerations  exhibited  resonance  in  the  2  -  5  Hz  frequency  range, 
and  the  peak  aortic  flow  responses  were  most  extreme  in  the  2-6 
!!z  range.  Resonance  of  these  quantities  in  these  frequency  ranges 
has  also  been  observed  experimentally9  and  In  an  open-loop  model^. 

5.  Model  responses  effected  by  nonsynchronoua  time-varying  accelera¬ 
tions  producing  long  beat  patterns  could  be  reproduced  synchronously. 

Recommendations  for  Improvement  of  the  Model 

1.  Active  right  atrium.  This  study  has  shown  that  the  effect  of 
accelerations  on  right  atrial  pressure  Is  important;  therefore, 
the  timing  of  accelerations  with  respect  to  venous  sytem  events 
may  be  an  Important  factor  in  changing  cardiac  output.  The  sensi¬ 
tivity  of  cardiac  output  to  right  atrial  pressure  In  the  model 
Indicates  that  the  timing  of  venour  acceleration  with  respect  to 
right  atrial  contraction  msy  modify  cardiac  output. 

2.  More  detailed  pulmonary  model.  This  addition  should  produce  better 
tilt  responses  in  the  model  since  the  lung  io  known  to  have  an 
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equivalent  heiglt^  and  might  produce  sustained  changes  in  cardiac 
output  with  zero-mean  time-varying  acceleration. 

If  the  inclusion  of  (2)  results  in  sustslned  changes  in  stroke 
volume  with  zero-mean  accelerations,  the  timing  of  left  ventricular 
filling  pressure  is  probably  important  and  an  active  left  atrium 


would  be  desirable. 


APPENDIX  A 


The  inclusion  of  a  piecewise- 1  inear  pressure-volume  character¬ 
istic,  postulated  for  the  upper  and  lower  venous  systems  of  the  model 
investigated  in  this  study,  could  conceivably  introduce  error  into  the 
performance  of  the  system  in  the  following  way.  Consider  the  plecewise- 
1  incur  pressure- volume  relationship 


V  -  PC 


(A- 1)  . 


The  derivative  of  Eq.  (A-l)  with  respect  to  time  is 


dV  =  CdP  PdC 
dt  dt  +  dt 


(A-2) 


where  —  represents  the  difference  between  flow  in  and  flow  cut  of  a 
dt 

section.  Since  C  may  be  discontinuous  in  time,  the  last  term  in  Eq. 
(A-2)  may  be  impulsive  in  nature  if  the  discontinuity  in  C  occurs  at 
P  ^  0.  If  volume  were  computed  by  direct  integration  of  Eq.  (A-2), 
step  changes  could  result  in  system  volume,  an  undersirable  condition. 
However,  volumes  were  computed  in  this  study  in  such  a  way  that  discon 
tlnuities  did  not  occur,  as  can  be  argued  by  considering  Fig.  5  (p  ). 
If  a  step  change  in  Vuvg  or  VLVS  is  assumed,  it  would  be  integrated 
and  produce  a  continuous  feedback  flow  signal;  thus,  the  assumed  step 
change  in  volume  must  he  a  consequence  of  impulsive  flow  from  the 
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preceding  section  of  the  model.  However,  flow  from  the  preceding  section 


cannot  be  impulsive  even  with  a  step  change  in  volume  in  that  section 

because  a  step  change  in  volume  would  produce  a  step  change  in  pressure 

and  therefore,  at  most,  a  finite  discontinuity  in  flow.  Applying  this 

dV 

reasoning  to  all  model  sections  shows  that  the  impulsive  ~  required 

dt 

to  produce  a  discontinuity  in  volume  was  not  generated  in  the  analog 
mode  1 . 
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In  order  to  determine  the  effect  of  external  forces  upon  a 
composite  section  of  flow  in  the  model,  consider  Che  axial  acceleration 
of  a  rigid  ibe  of  length  b  and  cross-sectional  area  A  in  which  a 
Newtonian  fluid  subjected  to  an  axial  component  of  gravitational 
acceleration  g  and  an  axial  pressure  gradient  (P^  -  P2)  i9  undergoing 
laminar  flow  (as  in  Fig.  19A).  A  force  balance  for  this  section  can 
be  written 


<Pt  -  P2)  A  -  W  -  F  •  m(V  +  Vw),  (B-l) 

where:  V  ■*  mean  velocity  of  blood  with  respect  to  wall, 

V,,  =  velocity  of  wall  with  respect  to  Newtonian  frame, 

w 

W  *>  weight  of  fluid  in  section,  and 
F  *  viscous  forces. 

Since  W  «  phAg,  where  p  is  fluid  density,  Eq.  (B-l)  becomes 

r , 

V  =  [(P1  -  P2)A  -  phAg  -  phA?w  -  Fjdt  (B-2) 

when  rearranged  and  Integrated  to  solve  for  V.  If  the  viscouc  forces, 
F,  are  assumed  proportional  to  the  velocity,  V,  then 


F  KV 
A  ”  A 


(B-3) 
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Area  =  A 


LI  .,T 

Figure  I3A.  Axial  Acceleration  of  a  Rigid  Tube 

Containing  a  Newtonian  Fluid  with  Direction 
of  Acceleration  the  Some  as  That  of  the 
Original  Fluid  Acceleration. 


Figure  I9B.  Axial  Acceleration  of  a  Rigid  Tube 

Containing  a  Newtonian  Fluid  with  Direction 
of  Acceleration  Opposite  That  of  the 

^  I  r  A  nnnlorfltiAh 

Ur  iuiiiui  r  iu*U  nvvv»»»wMVM. 


where  Q  =  VA  is  a  volumetric  flow  rate  and  R  =  — 7  is  interpreted  as  the 

A‘ 

fluid  resistance  of  the  section.  Making  these  assumptions  and  express¬ 
ing  L  *  as  the  inertance  of  the  section  produces  the  following  modl- 
A 

fication  of  Eq.  (B-2) : 


1  I 


Q  3  “  ■  [Pl  -  P2  -  Ph0v.  +  Vw)  -  ORjdt. 


(B-4) 


For  a  section  in  which  V  is  in  the  negative  2-direction  (Fig.  19B), 


]  r 

O  “  —  fP  y  '  P2  +  VW}  " 


(B-5) 


From  Eqs.  (B-4)  and  (B-5)  the  effect  of  external  forces  upon 
a  composite  section  of  flow  is  to  modify  the  flow  by  Instantaneously 
increasing  or  decreasing  the  pressure  gradient  depending  upon  the 
instantaneous  directions  of  the  external  forces  and  the  direction  of 
the  flow.  Since  accelerations  appear  as  pressure  modification  terms  in 
the  model,  they  we*re  simulated  by  series  voltage  sources  in  the  RLC 
vascular  sections  (Fig.  4).  The  polarities  of  all  sources  for  all 
acceleration  patterns  were  determined  from  Eqs.  (B-4)  and  (B-5). 

The  effect  of  the  earth's  gravitational  force  was  simulated  in 
this  study  by  the  conditions  of  55  0  and 

(a)  g  =  0  G  (prone  position) 

(b)  g  =  +1  G  (head-up  tilt) 

(c)  g  -  -1  G  (head-down  tilt) 

The  effect  of  time-varying  whole-body  accelerations  was  simulated  by 


* 


;r,d  V  a  time-varying  periodic  wave f on 
w 


iToaC  value 
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The  combined  effect  of  gravitational  and  time-varying  accelerations  was 
simulated  by  g  =  -  ~  G  (an  approximate  19°  head-down  tilt)  and  Vw  a 
time-varying  periodic  waveform. 
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